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O e-book: “Ciências exatas e da terra: teorias e princípios 2” é constituído 
por seis capítulos de livros que apresentam estudos de aplicação de modelos 
matemáticos, mudanças climáticas, importância das diferentes fontes de águas 
para fins potáveis e a presença de contaminantes que podem tanto afetar a 
qualidade da água, quanto se constituir em um meio no qual contaminantes com 
propriedades de modificar o sistema endócrino de animais e seres humanos. 

O primeiro capítulo apresenta o estudo de modelos matemáticos que 
regulariza dois sistemas de vigas termoelásticas definidas, conforme: o modelo 
de Timoshenko e a teoria de Green-Naghdiy do tipo Ill. Tais modelos se 
diferenciam por seus termos de acoplamento, como consequência da constitutiva 
leis consideradas no presente estudo. O capítulo dois se propôs a avaliar a 
estabilidade e regularidade do sistema constituído pela equação de vidração de 
terceira ordem de Moore-Gilson- Thompsom associada à equação clássica de 
calor com a lei de Fourier, demonstrou que tais modelos associados ao modelo 
MGT-Fourier são exponencialmente estáveis e analíticos. 

O terceiro capítulo investigou a caracterização e o comportamento de 
dados da Irradiância Global Horizontal (GHI) e Irradiância Difusa Horizontal (DHI) 
na cidade de Natal/RN, cujos resultados foram significativos ao se comparar o 
período de inverno e verão, com valores de 21,7% (GHI) e de 19,1% (DHI). 
Entretanto, para se identificar os fatores que ocasionaram este comportamento, 
faz-se necessário investigar a combinação das inúmeras variáveis que envolvem 
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em relação ao tema “Mudanças climáticas”, com ênfase na mitigação e 
adaptação que se constituem em um conjunto de ações a serem implementadas 
e desenvolvidas, com o intuito de reduzir a emissão de gases que contribuem 
para o efeito estufa e a minimização do uso de energia elétrica. 

O quinto capítulo avaliou a presença de diferentes classes de esteroides e 
suas respectivas capacidades de promover alterações nas atividades biológicas 
de organismos vivos, após entrar em diferentes ecossistemas aquáticos e 
promover inúmeras alterações nos diferentes organismos, sendo denominados 
de poluentes com propriedades de modificar o sistema endócrino. Por fim, o 
capítulo seis apresenta os conceitos de água subterrâneas e superficiais, bem 
como sua ocorrência, abundância e qualidade para fins potáveis. Além disso, 
realiza uma abordagem geral em relação aos principais aquíferos no Brasil. 

Nesta perspectiva, a Atena Editora vem trabalhando de forma a estimular 
e incentivar cada vez mais pesquisadores do Brasil e de outros países a 
publicarem seus trabalhos com garantia de qualidade e excelência em forma de 
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CAPÍTULO 1 


REGULARITY OF THE TIMOSHENKO'S SYSTEM 
WITH THERMOELASTICITY OF TYPE HI AND 
FRACTIONAL DAMPING 


Fredy Maglorio Sobrado Suárez 
Department of Mathematics, Federal 
University of Technology of Paraná, Brazil 


Lesly Daiana B. Sobrado 


Institute of Mathematics, Federal 
University of Rio of Janeiro, Brazil 


Filomena Barbosa Rodrigues Mendes 
Department of Electrical Engineering, 
Federal University of Technology of 
Paraná, Brazil. 


ABSTRACT: The article presents the study 
of the regularity of two thermoelastic beam 
systems defined by the Timoshenko beam 
model coupled with the heat conduction 
of Green-Naghdiy theory of type III, both 
mathematical models are differentiated 
by their coupling terms that arise as a 
consequence of the constitutive laws initially 
considered. The systems presented in this 
work have 3 fractional dampings:  (-A)" À, 
vA-My, and K(-A)*0, where &, y and O 
are transverse displacement, rotation angle 
and empirical temperature of the bean 
respectively and the parameters (T, 0, é) E 
[0, 1]º. It is noted that for values O and 1 of 
the parameter 7, the so-called frictional or 
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viscous damping will be faced, respectively. 
The main contribution of this article is to 
show that the corresponding semigroup 
S() = e&, with | = 1, 2, is of Gevrey class 
> E ET vio BE 
l+r'l+o'1+€ 

Ree= (0, 1º. It is also showed that S () = 
eºt is analytic in the region Ra=(T; O, 8) 
E [5, 1F) and St) = et is analytic in the 
region R, =((T,0, ) Ela, 1P /7= E). 

KEYWORDS: Gevrey class, Analyticity, 
Fractional damping, Semigroup theory. 


11 INTRODUCTION 


In this article, it is investigated the 
regularity of the semigroup associated with 
the thermoelastic beam system where the 
transversal vibrations are given by Timo- 
shenko's model (See Timoshenko [24]) 
and the balance of the energy is described 
by the Green-Naghdi theory, known as 
thermo- elasticity of type Ill(See Green 
and Naghdi [7]). The equation of motion 
and energy balance for this thermoelastic 
Timoshenko system is given by 


pAd—-S.=0 = 0, & 
plbu—-M:+S = 0, 


(0,L), te R*, 
ve(0,L),teR*, 


Capítulo 1 


pau tt +HIba = 0, zve(0,L), te R+. 


The constant p denotes the density, the cross-sectional area, and /the area moment 
of inertia. By S one denotes the shear force, Mis the bending moment and q is the heat flux. 
The function à is the transverse displacement, y is the rotation angle of a filament of the 
beam and uis the temperature difference. Here, tis the time variable and xis the coordinate 
space along the beam. 

In this investigation, two systems that differ only in their coupling terms are analyzed, 
which are obtained by starting from similar constitutive laws. 

The first best-known mathematical model in the literature, and studied for example in 
[5], [14] and [15], starts from the following constitutive laws: 


S = KAG(d, + UV), M = Ely, + Bu and q=-da; — Kas, (1) 


where a is the so-called thermal displacement whose time derivative is the empirical 
temperature u, i.e. a, = u, Eand Gare elastic constants, kthe shear coefficient for measuring 
the stiffness of materials (k < 1), ô and Kdenote the thermal conductivity, B the coefficient of 
linear thermal expansion and y a coupling constant. 

To simplify the notation let us denote by p, = pA, p, = pl, K=kAGand b= El. Under 
these conditions, the system can be written as 


próu—r(dr+Hy)a = 0, zve(0,L), teR”, (2) 
pode bbrr + (de +) + Bus = 0, ve(O,L), teR”, (3) 
paum — dUrr +HYbta — Kuro = 0, ve(0,L), te R”. (4) 


In order to exhibit the dissipative natural of system (2)-(4), it is convenient to 
introduce a new variable (see [25]): 


nt 
9(t,a) = / u(s, a)ds + px(e) (5) 
Jo 


where x E H;(0, L) solves the following Cauchy problem 


( Xez = paty — Kugra + VV in (0,1), 


x(2) =0, 2 =0,L. (6) 


Then, considering 3 fractional damping /;, (-A)" P, L1(-A)ºy, and K(-A)*6, and using (5) 
and (6) the starting system (2)-(4) is transformed to 


poe m(de +) +m(-A)Vo = 0, ze(0,L), teR”, (7) 
pote — bÃVy + (de + W) + BO +uo(-A) Ud = 0, ve(0,L), tEeR”, (8) 
p30u — 0A0 +Yybrr + K(-A)SO, = 0, zeE(O,L),tEeR*, (9) 


where the parameters 7, o and E take values in the range [0, 1] and yu ,, , are positive 
real constants. 

The second mathematical model addressed in this research that was proposed in 
2017 by Santos et al. [17], which consists of adjusting the couplings of the first system 
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couplings to satisfy type Ill heat conduction that acts mainly on the shear force, and as a 
result, the transverse shear force S = x(Q, + 1) is obtained from the transverse distribution 
of shear stresses acting on each cross-section of the beam. For this purpose, itis necessary 
to start with the used constitutive laws 


S = n(dr +) — uu, M = by, and q=-dar— Kas, (10) 


where a is the so-called thermal displacement, the time derivative of which is the 
empirical tem- perature u, i.e., a, = u. From these constitutive laws, by making appropriate 
changes, one arrives at the system to follow (for more details see [17]): 


pó rd ++ ud +m(-A) O = 0, re(0,L),teR*, (11) 
pode — bAp+n(dr +) uu +Uo(-A) db, = 0, ze(0,L),teR*, (12) 
p304 —0h0 + ud HU) +Y(-ASO, = 0, zEe(O,L),tER*. (13) 


To make the calculations in the results of the test more practical, the two systems are 
abstractly rewritten by using the operator: A: D(A) c LH0, L) —> LO, L), where 


A=-A=— e D(A)= H2(0,L)N HE(O,L). (14) 


tis known that the operator given in (14) is self-adjoint, positive, and has inverse 
compact on a complex Hilbert space D(4º) = LZ0, L). Therefore, the operator A” is self- 
adjoint positive for all v E R, bounded by v < 0, and the embedding 


D(A“) — D(A"2), 


is continuous for v, > v,. Here, the norm in D(A") is given by |ull5ay = || A“ul, u 
E D(AY, where ||.|| denotes the norm in the Hilbert space LO, L) = D(Aº). Some of these 
spaces are: D(A'?) = HI(O, L), D(4º) = LO, L) and D(A'2) = H(0, L). 

Studying the regularity (Gevrey class and/or Analyticity) of systems is relevant to 
science. In applied sciences, obtaining information about the solutions of the mathematical 
model under study such as regularity (smoothness) has the same importance as knowing 
the asymptotic behavior of the model's solution. From a mathematical point of view, when 
one talks about the regularity of solutions, concepts like differentiable, Gevrey class, and 
analyticity come to mind. It is already known that semigroups S(t?) = e*' of the Gevrey class 
have better regular properties than a differentiable semigroup, but are less regular than an 
analytic semigroup. The Gevrey rate s > 1 'measures' the degree of divergence of its power 
series. It should be noted that the Gevrey class or analyticity of the modelin particular implies 
three very important properties: The first one is the property of the smoothing effect on initial 
data, i. e., no matter how irregular the initial data is, the solutions of the models become very 
smooth in positive finite-time. The second property is that systems are exponentially stable. 
Finally, these systems benefit from the property of linear stability, which means that the type 
of the semigroup is equal to the spectral limit of its infinitesimal operator. 

During the last decades, various authors have studied some physical phenomena 
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for the Tim- oshenko beam system formulated in different mathematical models. Most of 
them focused on studying the asymptotic behavior, always trying to obtain the best decay 
rate and using dissipa- tions in some of the system equations. In the following paragraphs, 
some of these investigations are mentioned. 

In 1987, Kim and Renardy [9] studied the asymptotic behavior of the Timoshenko 
beam consid- ering two boundary dissipations. They demonstrated the exponential decay 
of energy associated with the model using the multiplier technique and also determined 
numerical estimates of the eigenvalues of the operator associated with this system. Later 
in 2008, Messaoudin and Said- Houari [14] also studied the asymptotic behavior of the 
Type Ill thermoelastic Timoshenko system with mixed boundary conditions (Dirichlet- 
Dirichlet-Neuman), they demonstrated using the en- ergy method that if the velocities of 
the waves associated with the hyperbolic part of the system are equal, then the system 
decays exponentially. The complement of the study of the asymptotic behavior of this same 
system due to different mixed conditions (Neumann-Dirichlet-Neumann) was studied again 
in 2013 by Messaoudin and Fareh [15], in this new research they show that if the speeds of 
the waves associated with the hyperbolic part of the system are different, then the system 
decays polynomially. In the following year, the work of Fatori et al. [5] who also studied 
this same system, but with two types of boundary conditions (Dirichlet-Dirichlet-Dirichlet) 
and (Dirichlet- Neumann-Dirichlet), using semigroup technique they showed that the 
corresponding semigroup is exponentially stable If and only if the velocities associated with 
the hyperbolic part of the system are equal, in the absence of exponential decay they show 
that the corresponding semigroup is polynomially stable and the determined rate is optimal. 

More recently, in 2017, Santos et al.[17] studied the type Ill thermoelastic Timoshenko 
beam system given by: 


pride — (de + Ur +uO a = 0, zeE(O,L), te R”, 


pote —bbra + rd +U)—uO, = 0, ve(0,L),t 
P30 — d0 pa = (Da + w) — YOtax = "0, TE (0, L), te R”, 


mM 
x 


and the boundary conditions (Dirichlet-Dirichlet-Dirichlet) given by 
(0,8) = G(L,0) = y(0,t) = W(L,0) = 0(0,t) = 0(L,t) = 0, vt>o0, 
or with the boundary conditions(Dirichlet- Neumann-Neumann) given by 
A(O, t) = 9(L,0) = Ya(0,t) = Úr(L,0) = 0:(0,t) = 0:(L,t) = 0, vt>0, 
note that this system differs slightly from the terms of coupling. In this research, the 
authors also applied the semigroup technique and showed that the associated semigroup 
can be exponen- tially or polynomially stable depending on the relationships between the 


wave propagation velocity coefficients, specifically in the case of polynomial decay it was 
proved that the rate found is optimal. 
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About investigations the regularity of semigroups associated with various 
mathematical mod- els, one could mention the paper/article of Fatori et al.[4], whose work 
the authors studied the differentiability and analyticity, in addition to studying the asymptotic 
behavior via semigroups. For analyticity, they used a theorem that can be found in the 
Liu-Zheng book[12] or the work of Hao et al.[8]. Other more recently published works have 
explored the regularity of solutions using the Gevrey class introduced in Taylor's thesis [19] 
(1989). In the same direction, one could mention [1, 2, 6]. 

Among recent research that has studied the asymptotic behavior and/or regularity of 
models with fractional damping, one could mention the works of Sare et al.[18], in that paper, 
the authors investigated coupled systems thermoelastic type, where they address two cases 
with Fourier's heat law and the other with Cattaneo considering in both cases the rotational 
inertial term. Further- more, they, study the exponential stability, possible regions of loss of 
exponential stability and polynomial stability, and, more recently, the work of Keyantuo et 
al.[20](2020) to be published. In this latter work, the authors studied the thermoelastic plate 
model with a fractional Laplacian between the Euler-Bernoulli and Kirchhoff model with 
two types of boundary conditions. In ad- dition to studying the asymptotic and analytical 
behavior, the authors show that the underlying semigroups are of Gevrey class s for every 
s > io for both the clamped and hinged bound-ary conditions when the parameter 6 
lies in the interval (0, 1/2). Moreover, one could cite the investigations [3, 10, 11, 16, 22, 23]. 

One motivation for deciding to study the regularity (determination of the Gevrey 
classes and analyticity) of Timoshenko's systems was in the direction of complementing the 
work published in 2005 by Raposo et al. [21], in this work the authors study the asymptotic 
behavior (exponential decay of the Timoshenko system), the studied system considers two 
weak (frictional) dampings given by the speed of the transversal and rotational movements. 
The investigated model is given by 


prumo k(us— Ur +tu=0 im (0,L)x (0,00), 
pode bre + ku) +y=0, in (0,L) x (0,00), 
u(0,t) = u(L,t) = MO,t) = Ú(L,t) = 0 t>0. 


The authors use the frequency domain technique (spectral characterization) to study 
the asymp- totic behavior via semigroups. 

This article is organized as follows. In section 2, using the theory of semigroup, 
the well- posedness of both systems is studied. Section 3 is dedicated to the study of the 
regularity of both models, this subsection is subdivided into two parts: in subsection 3.1 
Gevrey class of the semigroup associated with the first system S(t) = e!” in the region Rj. 


= (0, 1]º is determined and it is shown that S(t) = eº' is analytic in the region R,,. Finally, 
in subsection 3.2 we show that the Gevrey class of S(t) = eStis S > TS SR Ee for 


the parameters (T, 0, €) within the region AR, and this subsection is finished by proving 
that S(t) = eS% is analytic when the 3 parameters T, o, and E take values in the closed 
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interval [3, 1] such that T = E. The investigation ends with an observation regarding the 
asymptotic behavior of S(t) = eº' in which it is illustrated that the necessary estimates given 
in Propositions(5) and (6) imply that the semigroup associated with the system (15)-(17) 
and (21)-(22) and the system (18)-(22) respectivamently are exponentially stable for (T, o, 


Selo, 17º. 


21 WELL-POSEDNESS: SEMIGROUP APPROACH 
Using the A operator defined in (14) the system (7)-(9) will be given by 


pride (de HU) HMA O = 0, zeE(O,L),tEeR*, (15) 
pobe +bAb + R(O + UV) +80 +uÃ db = 0, xe(0,L),teR”, (16) 
p30 + 040 + YU + KASO, = 0, ze(O,L),tEeR?*. (17) 


And again using the Aoperator on (11)-(13), leads to 


proud HW + +mA O = 0, ze(O,L),teR”*, (18) 
pobe +bAv+Rr(d +HY)—-uO+uÃ Cd, = 0, zxe(O,L), te R”*, (19) 
p30u + 040 + ud, HU) +YASO, = 0, ze(O,L),tEeR*. (20) 


Additionally, the appropriate initial conditions for both systems are considered 


a(o, -) = 00. valo, -) = qd, ú(O, -) = Vo, ú(O, -) = a, 90, -) = do and B(0, -) = 61, (21) 


and boundary conditions for both systems of type Dirichlet-Dirichlet-Dirichlet 
LD =ALLD)=0, EO) =UA(tL)=0, 9(t0)=0(L.0)=0 t>0. (22) 
Taking = 4, V=y, O=0, considering U=(d, À, y, 4, 0,0)andU, =(d, PD, YU 
W, 9, 8,), both systems, can be written in the following abstract framework 


SU() =BU(t), U(0) =Uo, (23) 


[0 


where for the operator B, for the first system is defined by 


p 
| pride + Wa — mi 47] 
y 
2 TT um 9 
BU: | ba + (és 4.9) +90, + 4º] (24) 
[0] 
| —1 (540 + pb, + KA£O] | 


and the operator B, for the second system is defined by 
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pb 
o lk(dx + W)r — UO — 11, 47 O] 


v 95 
BU := a I-bAY — n(de HU) +UO — po AV] | (25) 
[9] 
os [0.40 — u(D, + W) — y 450) 
forU=(d, db, y !W, 6, O). Both operators are defined in the energy space 
Hi=H2=H:= [D(A?) x D(Aº). 
H, is a Hilbert space with the inner product given by 
(U,U dm = pnBD, DS) + pory(P, VS) + Bryldr + V, 6; + UV) + beylbe, Vi) 
+808(07,07) + BrK(O,0%, 
forU=(4, 4 4, 0,0,0), U'=(4, 4, y,W', 0,0) E H, and induced norm 
JU, = paByje + pay]? + Balde + VI] + dry] AZ U]]? + Bôs|A?0]2 + Brkjo!. 


In these conditions, he domain of B, is defined as 


D(B,) = f UcH/B(U) e H+ boundary conditions (22) | 
= D(A)x D(A)ND(A?) x D(A) x D(AZ)N D(A?) x D(A) x D(A?) N D(AS). (26) 

And H, is a Hilbert space with the inner product given by 
(DU = plBD) + poa + nlde +, OLA UVD + blbe, vT) + 5(0,,07) + pa(0,0%), 
for U = (6,94, W,0,0),U* = (98, Dt, W*,0*,0%) E H, and induced norm 

[Ul = pullBIP+ polly]? + nllós + W]2 + b] 424] + 6] 429]? + pal! 

In these conditions, the domain of B, is defined as 
D(Bo) := D(A)x D(AT)ND(A”) x D(A) x D(AZ)N D(A?) x D(A) x D(A?) N D(AS). (27) 

To show that both operators B, i = 1, 2 are generators of a C, semigroup we invoke 
a result from Liu-Zheng' book[12]. 

Theorem 1 (see Theorem 1.2.4in [12]) Let B, be a linear operator with domain D(B) 
dense in a Hilbert space H. IfB, is dissipative and O E p(B), the resolvent set ofB, then B, 
is the generator of a C,- semigroup of contractions on H. 

Next, we apply the Theorem 1 only for the first system, the proof of the second 
system is completely similar. Let us see that the operator B, defined in (24) satisfies the 
conditions of this theorem. Clearly, D(B ) is dense in H. Effecting the internal product of B,U 
with U, lead to 


BnkK? 


Re(BiU,U) = —Byn||AZ6]2 — sypo|| A? W]|? — I4ÍO|P, vUCD(B), (28) 


P3 


that is, the operator B, is dissipative. 
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To complete the conditions of the above theorem, it remains to show that O E p(B,). 
Let F= (f", f2, fº, Fº, fº, f) EH, let us see that the stationary problem B,U = F has a solution 
U=(4, & wy, IP, 0, O). From the definition of the operator B, given in (24), this system can 
be written as 


p=, (de +) = nf + mA fl, (29) 
vo fS, bAy + eldr +) = —paf! — no 4º fº — 8º, (30) 
O=º, 040 = —yf3 — KAS fº — ps fº. (31) 


Therefore, it is not difficult to see that there exists only one solution À, y and O of 


the system 
(de +) = pf + mA f € D(Aº), 
DAY + n(de +) = —pof! — no A fê — Bfô e D(Aº), 
040 = —yf3 — KAS fº — pa fê e D(Aº). (32) 
Therefore: 
JUlly < ClFlla, 


which in particular implies that ||B,'F ||, = C||F ||, so we have that O belongs to 
the resolvent set p(B,). Consequentily, from Theorem 1 we have B, as the generator of a 
contractions semigroup. 

As a consequence of the above Theorem(1), it follows 

Theorem 2 Given U, E H there exists a unique weak solution U to the problem (23) 
satisfying 


U e C([0,+00), H). 
Futhermore, ifU, E D(B9, KEN, then the solution U of (23) satisfies 


K 
U e (| cti(lo, +00), D(B)). 


j=0 


31 REGULARIZATION RESULTS 


This section is divided into two subsections, in each of them, the Gevrey class and 
analytics of the corresponding systems are studied. Using the characterization results 
presented in [20] (adapted from [19], Theorem 4, p. 153]) it is shown that the corresponding 


£ 
é 


semigroups S(t) = eS! are of Gevrey class s > é for r=min Liz: tis» Tre) E (O, 1), with (7, 0, 


8) E Rog := (0, 1)º. And for the study of analyticity the main tool used is the characterization 
of analytical semigroups due to Liu and Zheng[12] (See the book by Liu-Zheng - Theorem 
1.3.3). Ithas proved that S,(?) = eº"'is analytic in the region R, ;= [1/2, 1º and S() = e*tis 
analytic in the region R,:=((T, 0, DE [1/2,1]%/T = E). 

In what follows: C, C, and C, will denote positive constants that assume different 
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values in different places. 


3.1 Regularity of the first system 


Next, two lemma are presented where two estimates are tested which are fundamental 
for the determination of Gevrey class and the analytics of the associated semigroup S(?) = 
e!! t 

Lemma 3 Let S(i) = e*'! be a C,-semigroup on contractions on Hilbert space H=H,, 
the solutions of the system (15) (17) and (21) (22) satisfy the inequality 


lim sup ||(1X — By) Hen) < 00. (33) 


JAl>0e 
Proof: To show the (33) inequality, it suffices to show that, given O > 0 there exists 
a constant C, >0 such that the solutions of the system (15)-(17) and (21)-(22) for AI > & 
satisfy the inequality 
[Ux 
HPllx 
HAERand F=(f!,f2,f3,fº,fs, fe Hthenthe solution U= (4, À y, 4,0,0) E 

D(B,) the resolvent equation (iAl- B)U = F can be written in the form 


< Cs > |Ulu = (GA — By)" Fly < CollFlly- (34) 


Nó-D = fl in D(A?), (35) 
Do - (do +tuy) HAS = f2 in D(A?), (36) 
pa p1 
vv = Sin D(A?), (37) 
b a? B 9 
nv+Ay+ Edo +W)+ LO, +AV = fin D(A?, (38) 
P2 P2 P2 P2 
iNN-O = f im D(A?), (39) 
; Ô Y Ke 6 0 
NO + —40+—V, + — 4£0 fº in D(A?). (40) 
P3 P3 P3 


Using the fact that the operator B, is dissipative, result in 


z 9 2 2 o) 2 Ê 2 
Brum AF O]? + nyna] AS |? + So jAdO!? 
E) 


= Re((iA — Bi)U,U) = Re(F,U) < |F|uUlu. (41) 
On the other hand, performing the duality product of (36) for Byp, ), and remembering 


that the operators A” for all ve R are self-adjoint, resulting in 


Bed + W), dr) = Bypil|BI||2 + Bypi(D, f!) 
— im BY] AZ 6]? + a By(AZ FI, ATO) + Bypilf2, O), 


now performing the duality product of (38) for Byp,y, and remembering that the 
operators A“for all v E R are self-adjoint, leading to 
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Br Alda +), 0) + bBY]LABUIÊ = Bypa ]IVIÊ + Bypa(B, 1º) 
+ BO, 14) — iNByuo|| AF]? + Byuo( AF FS, AT) + Bypal FE). 


Adding the last 2 equations, result in 


Bryllóz + WI]? + b37 APV]? = 8yp1||B]|2 + Bypa | PJ? — By fm || AF ul]? 
+ o] AF UP + Bypo(D, Fº) + Bypo(W, fº) 
+ By (AF FI, ATO) + Byuo( AS f3, ATO) 
+ By, 4) + Bypalf!,b) + BO, va). (42) 
Taking real part, using norm ||F||, and || U||, and applying Cauchy-Schwarz and 
Young Inequa- lities, for e > 0 exists C, >0 which does not depend on À such that 


Bryllôe +] + bBy] AZ VI]? < BypiB|2 + Bypa |V]? + el 4? VI]? + Ce |O]? + Coll FIeU le. 


From estimative (41) and the fact 0 <Z, 0 <% and O < 3, the continuous embedding 
D(Aº ) — D(Aº!), 6, >6,, lead to 


Bnyllóe + WI] + bByAZWIP < CoFIalUn- (43) 


On the other hand, performing the duality product of (40) for Bxp,8, and remembering 
that the operators A for all v E R are self-adjoint and from (39), result in 


Brô||A20]2 = ô8x||O]2 + Brps(O, fº) + Bry(T, 04) 
— inBKAI|A20]2 + K8K(A? fº, 420) + Bxpa(1º,0) 


Taking real part and by using the inequalities Cauchy-Schwarz and Young and norms 
Fl), and |U], from IA >1, for e >0 exists C, >0 which does not depend on À such that 


Bróll4to|? < cre? + Wells + 4? FAO] + IfSIoIIy + CeIvIP + ejj4Zo|. 


Then, from estimative (41) and the factO <Z and 0 <5 < > the continuous embedding 
D(Aº2) >D(4º?), 8, > 6, lead to 


Bnô|A20]2 < Coll FIullU II. (44) 


Finally, ftom estimates (41), (43) and (44), the proof of this lemma is completed. 
Lemma 4 Let ô >0. There exisis a constant C, > 0 such that the solutions of the 
system (15) (17) and (21) (22) forlAl > ô satisfy the inequalities 


() MBlxllóe + 4]? + dIAZ UI] < [8y/pi BI] + pa |PI2] + Cs [Fa |U ||u- (45) 
(ii) BxoJaIAZO]? < BrpslAJOIÊ + Cs|lFIlu Un. (46) 


Proof: Item (i): Performing the duality product of (36) for Byp, AQ, and remembering 
that the operators A” for all ve R are self-adjunct, result in 
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BryM(de +), dr) = BypID|P + Byo(AB, Fº) — im By] AF DI 
in BATO, 1) + iByor($2,D) + iByo( PF"). 


Now performing the duality product of (38) for Byp,AW, and remembering that the 
operators A“for all v E R are self-adjunct leads to 


ABylik(óe +), 10) + bI AP VIP] = 8703 API + Bypa (NV, Fº) + 18240, va) 
+iB"y(O, f2) — ino By] AF VIP — iuaBy(ATY, fº) 
+iByoaf" DV) + iBypalf!, Fº). 


Adding the last 2 equations, result in 


AB/[sllôs + VIP + dA W]2] = 284/91 2]]2 + pa] VIP) — àBy(ma || AF SI] 
+ po AFI + Bypi(AD, 1) + Bypa(NE, PS) — By (ATO, f1) 
+ pa(ATE, PS + iBypl( 2, D) + iBypalf E) + iBypul ff!) 
+ iBypaf, 1º) + BO, a) + IB y(O, 8). (47) 


On the other hand, from (36) and (38), result in 


Bylor(aB, 1!) + pa(AP, Fº) = iByln((do + 14), F3) 
+u( ATO, 1) — pl, + (AS, AR PS) + (lda + 0), Fº) 
— B(O, fz) + uol ATY, Fº) — pol ft, fº)). (48) 
Using the identity (48) in the (47) equation and simplifying, lead to 


Aylkllós + WI]? + dIAZ W]P] = A8ylor |]P + pal PP] 
— (mi || AF OP + pa] AF UIP + iByn(de, f2) + iByn(y, f2) 
+ by (A? Ú, A? PYiByn(d, FS) +iBya(y, Fº) 
+iByo($, 6) + iBypalf! VD) +iBAO, ba). (49) 
As for E >0 exists C >such that lip?y O, pl =|| A2y||2+C,|0]|*, taking the real part, 
and applying the Cauchy-Schwarz and Young inequalities, estimative (34) of Lemma(3) and 
using the definitions of the F and Unorm in (49) the proof of item (1) this lemma is completed. 


Onthe other hand, performing the duality product of (40) for Bxp,AO, and remembering 
that the operators 4º for all O E R are self-adjunct, result in 


BrôoMIAZO]? = Brps(AO,iAO) — Bry(b,,A0) + BrK(AEO,ÃO) + Brps(f8,A9). (50) 
As 
pa(MO,iMO) = paO]2 + pa(ÃO, Fº) 
= paJOI)? + i5(420, 42 95) — i(W, 15) + iK(ASO, 1º) — ipa(fS, 5). (51) 
— Ve; Ab) = —iyttr; O) = iylba, 1º). (52) 
K(ASO,M0) = -iKI||AÍO]2 — iK(ASO, fº). (53) 
pa(f0,A0) = ipalf8,0) +ips(f8, Fº). (54) 


Using (51)-(54) in (50), lead to 
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8nôX||A20]2 = 8xp3A|O|2 + i2x6(A20, A? Fº) — iBRy(D, fº) — iBRy(ys, O) 
—iBrylbe, 1º) — iBRKI|AZO]2 + iBxps(1º,0). (55) 
Taking real part of the equation (55) and using the inequalities Cauchy-Schwarz and 
Young and norms ||F|, and ||U|),, from IAl >1, fore >0 exists C, >0 such that 


Brô|M|I420]2 < BxpsAIJO|]? + CL AZ0]| A? FS] + |U SS!) 


+ Ilha |]º + JO]? + [lb 111 $>]] + 11 ]]9]]). 


Then, ffom Lemma(3) and as ||y,||? = || 22 ||?, the proof of item (ij) of this lemma is 
finished. Q 


3.1.1 - Gevrey class of the first system 


Definition 5 Lett, 0 be a real number. A strongly continuous semigroup S(?), defined 
on a Banach space H, is of Gevrey class s > 1 fort>t, if S(t) is infinitely differentiable for 
t>t, and for every compact set Kc (t,) and each | >0, there exists a constant C = C(u, 
K) >0 such that 


[SC (Ile <Cu"(nl)*, foral te Kn=0,1,2.. (56) 


Theorem 6 ([19]) Let S(?) be a strongly continuous and bounded semigroup on a 
Hilbert space . Suppose that the infinitesimal generator of the semigroup S(t) satisfies the 
following estimate, for some O <r<1: 


lim sup |AVI|(GGAI — B)"lcça) < 00. (57) 
|A|- 00 


Then S(t) is of Gevrey class s fort>0, forevery s >. 

The main result of this subsection is as follows: 

Theorem 7 The semigroup S(t) = e”! associated to the system (15) (17) and (21)- 
(22) is of Gevrey class s for every s > : for r= minti&- 15: x) E (0, 1), with (T, 0, E) E 
Rea:=(0, 1º. 

Proof: From the resolvent equation F = (iAl-B)JU for A ER, U=(iAl-BJF. 
Furthermore, show (57) of theorem(6) itis enough to show, that for e >0 exists C, >0, such 
that: 


AFI(iA — By) 1F 
Apa Fio Es <Cs AF IUIa <CsllFlla, FeH and G>0. 
| H 


Equivalent to 
PIU <CotlFIulU|u) for re(0,1). (58) 


where r= min (ri- TÊs» da | O<tT<1, O<osiandO<Es<1. 
Next, |AT= |[d]| + |A[7= nv + |N = |O] will be estimated. 
Let's start by estimating the term | ||$]|: It is assume that IA > 1, some ideas 
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could be borrowed from [13]. Set 9 = À, + &,, where À E D(A) and &, E D(Aº), with 


Ad, + Ad, = f2, nd, = Lad + Ep, — PLATD 4 AB. (59) 
pa p1 pa 


Firstly, applying in the product duality the first equation in (59) by Q,, then by AQ, and 
recalling that the operator Ais self-adjoint, resulting in 


AB? + 1142042 = (92,01), (60) 
taking first the imaginary part of (60) and in the sequence, the real part in addition to 
applyingCauchy-Schwarz inequality, results in 
DB < CF and 142%] < CHF. (61) 


1 


Applying the 4-2 operator on the second equation of (59), result in 


DATtO, = DAto+ DA ty, — PLA IO + AjO, 


Pa Pa Pa 


then, as ||A-Iya||2 = (-A*W,z, A 3Y) = (At, 4-3) = ly]? and T—3 < Staking 
into account the continuous embedding D(Aº2) — D(Aº'), 0, >6, and using (60), result in 
IMIATE dal] < C(IIAZ 6] + 1436] + by + 114204], (62) 
using (34), (41) and second estimative the equation (61), lead to 
43 6a]? < CIP IaalIU a + IL FIZO: (63) 
On the other hand, from &,= P — &,, (41) and as 5 <> the inequality of (60), resultin 
46]P < C(AFO|P + |AFO|P) <C(FIa Ux + FIGO (64) 
By Lions' interpolations inequality O E [=4 , 3 | resultin 
[da ]2 < O(A 3 05|2)7F [AF 6a ]P)T+. (65) 
Then, using (63) and (64) in (65), for IAI > 1, resultin 


12] < CIT (Pa lU e + [Fo (66) 


Therefore, as | ||2< || &,||º + || &,|| from firstinequality of (61), (66) and |A-? < | 


, result in 


ASI <CFIa + MIPP<CIFIFIR for 0<T<L (67) 


On the other hand, let's now estimate the missing term || ||): It is assumed 
that >1. SetW=W +W, where lW E D(A) and W E D(Aº), with 


Ab + Ab =f! and Nba Ay — Eds — Ep — DO, — PAP 4 AP). (68) 
P2 P2 P2 P2 P2 


Firstly, applying the product duality the first equation in (68) by UW,, and as the 
operators A” forall v E R are self-adjoint, result in 
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APPA + AP Wa|2 = (14,04). (69) 


Taking first the imaginary part of (69) and in the sequence, the real part in addition to 
applyingCauchy-Schwarz inequality, results in 


AP | < CHF Im |4ºWi|| < CHF. (70) 


As follows from the second equation in (68) that 


b K B E ; 
AA ID, =—— Aly — — A Id, — — A Ip— 420, — ATI 4 AIG, 
p2 P2 P2 p2 2 
then, as [|AZ6,|? =(-A20,1, 4726) = |ló|2, |4720,]2 = |O]? ando — 1 < staking into 
account the continuous embedding D(Aº2) —> D(Aº'), 0, >0,, lead to 
IA Wal] < CAIA 6]] + |hól+ |O] + AF WI) + 4284] (71) 


Using estimative (34) and estimative (68), for IAI > 1, resultin 


A Wa|P < CINE IaU Ia + FIGO (72) 
On the other hand, from W, = W — W (41) andas 5 < > the inequality (70), result in 


I4fwalP < C(Afy|?+ AFP) <C(] 


Pla UI + IPO) (73) 


Now, by Lions” interpolations inequality O E [- 5, % ], lead to 


[Wa]? < C(||ATÉWo||2) = (|| AF Wa) T+. (74) 


Then, using (72) and (73) in (74), for IAI > 1, resultin 


Pal|É < CIATHS (FIO e + 11FIIoo)- (75) 


Therefore, as ||W||2 = |W||º + ||tv||º from inequality (70), (75) and IA? <IA 15, 


result in 


AT |EI<ClFIu + MIyê<S|FI for 0<o<l (76) 


Finally, let's now estimate the missing term jo]: Now we assume IAI > 1. 


SetO=0,+0, where O, E D(A4) and O, E D(4º), with 
ha) 4 K 


iÃO, + AO, = fº and NO, = AG- AO + AO. (77) 
P3 P3 P3 


Firstly, applying the product duality the first equation in (77) by O, and as the operators 
A“forall v E R are self-adjoint, lead to 


iX91]2 + | 4201]? = (48,09). (78) 


Taking first the imaginary part of (78) and in the sequence, the real part and applying 
Cauchy-Schwarz inequality, result in 
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MIO<ClFly and ||4/0;/]|<CI]Flk. (79) 


As follows from the second equation in (77) that 
e) ô a Y 4-l Kes 1 
AA 20, =—— AZ) — — A 2W,— — AS 204 420, 
P3 P3 P3 
then, as |A-> W,|2= ||W|?and =; <0, E-3 < 5 taking into account the continuous 


embedding D(4º2) > D(4º!), 0, >6,, resultin 


MIII4-265] <CH420] + | 420] + || + [4063] (80) 


Using estimatites (34), (41) and second estimative of equation (79), for IAI > 1, lead 
to 


4728,]" < CH 


[FlaUIa + Flo (81) 


On the other hand, ffom O, = O — O,, (41) and as 
(79), resultin 


toym 


< > the second inequality of 


sie) < cil4te|?+ AO) <C(PIUly + ILFIZO: (82) 


£ 
5 


, 3], resultin 


Now, by Lions” interpolations inequality O E [- 
182]? < C(I4-26,]2) 7 (142 Oo |P)r. (83) 
Then, using (81) and (82) in (83), for IAI > 1, resultin 


1922 < CAIS (PIU Ia + [FI (84) 


Therefore, as |0||2< |o|]º+ |0])> from first inequality of (79), estimative (84) and 
IN? < INT, resultin 


Ao <CFlu <> |oB<CHNFÉ]FIR for O<E<1 (85) 
Using the estimates (67) and (76) in the inequality item (7) of Lemma(4), result in 
Bysllóo + VIP + ByblAZV|P < Co(NTE + [NEM] + CoNHFIallUla. (86) 
Furthermore using estimative (85) in estimate (46) item (ii) of the Lemma 4, we have 
Bnôjl Ato]? < Brps|JO]2 + Co THE IeUlu < CONTEI]? + CONFIO. (87) 


Finally summing the estimates (67),(76), (85),(86) and (87), we have 


[UT < Co (AIEA + NÉ + CO INCH 
as for (1, 0, E) E [0, 1, we have -1 = 2 max(—-E=,— Tê» Tie) = -2 min 
T € 
(ir Tio» THE), then 
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JU, < cristo | FlaU la + IFIB, 
Equivalently 
pqmintrE-otês Tie JUlln < ColFly for (7,0,€) e [0,1]. 


As r= minto Tg) Ee (0, 1)for(t,o, &) E (0, 1]º the proof of this theorem is 
finished. 


3.1.2 Analyticity of S(t) = e*! for(t,o, Dela, 17º 
Next, we announce the Liu-Zheng Book Theorem that will be applied in the proof of 
the analyticity of S(t) = eº! for i= 1, 2 in their respective regions. 


Theorem 8 (see [12]) Let S(t) = e”! be C,-semigroup of contractions on a Hilbert 
space H. Suppose that 


B)D (i)N/AER)=iR (88) 
Then S(?) is analytic if and only if 


lim sup ||A(iA — B)"!le(m) < 00 (89) 


|A|—+00 
holds. 


Before proving the main result of this section, the following two lemmas will be 
proved. 


Lemma 9 The C,-semigroup of contractions S(t) = eº"' satisfies the condition 
p(Bi)D (iX/AER)=iR. (90) 


Proof: Let's prove iR c p(B,) by contradiction, it is supposed that R é o(B,). As 0 E 
P(B,) and p(B,) is open, the highest positive number À, is considered such that the interval 
(-iA,, iA,) c P(B,), then iA, or -iA, is an element of the spectrum o(B,). Itis supposed that iÀ, 
E o(B,) (if -iA, E o(B,) the proceeding is similar). Then, for 0 <v <A, there exist a sequence 
of real numbers (A ), with 0 <v <A, <A, À, — À, and a vector sequence U =(b, À, W, 
W, 0,0) €D(B,) with unitary norms, such that 
(in = Bi)l nllH = — Falla > 0, 


asn—> 0, From (43)for0<tT<1andO<o<1, resultin 


« À, ) Y - 
Br9 den + Un |? +b8y]A2Ym|" < ColBalalUnla- 
From (44), for 0 < & <1, resultin 
BrôA20n|2 < Cs] Falha llUnlly- (91) 


In addition to the estimative (41), for0 <T< 1 and0<o=<1, resultin 
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BK 


Brym | AÉ On? Aê On JP < Co] FnllyllUnlla- 


n 


Consequently, || U,|2— O. Therefore, lead to |U. ||, — O but this is absurd, since 
|u,|,=1foralneN.Thus, iR c p(B). This completes the proof. 


Lemma 10 Letô >0. There exisis a constant C, >0 such that the solutions of (18)- 


(20) and (21) (22) forlAl > satisfy the inequality 


. ç vuiida . . 1 

() IMISIÊ<SCoFIaUla for S<r<1 (92) 
(ii) MIÊ<CFIalUlu for S<o<i (93) 
a 2 . 1 
(iii) |MJOJP<CsFIulUlu for S<E<1 (94) 


Proof: Item (/): of this lemma will be initially shown, performing the duality product of 
pm 
(36) for ;, AT AQ, and recalling that the operator A“ is self-adjoint for all v E R, results in 
aIójl= É(M-A6+W,), 4776) + PL ÃO ATO) ii NA EO]? 
Ha ta 
= iÉMA6+Af!),ATTO) — E, +f3,ATO) 
nm 


ta 


ess do de Rê 
+— (final g — RÃ ba + ind — ipi AT” fo) — PINA TO]? 
Ha ta 


=iE IA o]? ti mt At FLA TO) +il(V,ATO,) —i(f3, ATO) 
Ha ta ta 
(12,0) Ro aa ALMA, va) + Acd 2 
1 
A EO], 
Noting that: |ATo,|?=(470,,474,) = (“A 4, 476) =| Ao], taking real part and 
considering that 3 <T < 1 using estimative (41) and using Cauchy-Schwarz and Young 
inequalities, for e >0 exists C,, such that 


MIBI? < ColFlaUly + cllATT Oo] 


As O < 5, then from estimative (41) |W|2 <C,|F||,||U|),. From the continuous 
embedding for IAI > 1, the proof of item (i) of this lemma is finished. 

Item (i): Again similarly, performing the duality product of (38) for E ASAUW, using 
(37), and recalling the self-adjointness of Ay, v E R, leads to 


c t v n) 2.49 E) c 
AIPI2 = (AN, 4-0) — E(Mór +14), ATT) + PE (f4 MAO) — PRA EV 
E diz tz uz 


= AE? + EA? 3 f3, AO) Hi (ACD, W) + io (11, ACO) 
[E 


E) 


2 b 
na Top+riÊ(, AD) if AL) (FIA Po) 
[E [AD P2 btz 


.K =g .B -—e - .P2 -g 2 .P2 4,92 EE; 5 
“(A U) = (IA O) = 4", 8) + if A EP IA Ty? 


ua 
Noting that: ||4=" 9]? = |A-7%,|2, taking real part, considering 3 < o < 1, using 
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Cauchy-Schwarz and Young inequalities, for e >0 exists C, such that 


1-29 
p) 


API? < ColFlulUly +C.|A7 


6I/* + elpyI2. 


As 2º <0 ++ | <o<1. Considering IA > 1, using estimative (34), norms ||F 


2 


|, and || U|),, the proof of item (ij) of this lemma is finished. 
Item (iii): Finally, let's prove the (iii) item of this lemma, performing the duality product 
of (40) for e AO, using (37), and recalling the self-adjointness of AY, v E R, resultin 


2 9314 -£and O nal£and o. .Ô, 4165 4d 0,6 41I& 
AO!" = INC] AT20]|* —I ATZ O]|* —(AZfº AISO) — i— (Fº, Al-EQ 
je] eb jA-doI ri SIA OP + ip (Ats ) = ip (ro, ATO) 


+L(yplv, 4-te,) 12. (48, 4-t4,) — (19,0) +i PEA S foj2. 
K [A K K 


MM 


Taking the real part, applying Cauchy-Schwarz and Young inequalities, for = > 0 
exists C, > O such that 


MIS < CEIA Te]? + A? PII A?-S6]) 


+ FILA! Sol + If I ASP | + CoAIPE|? + eJAj4“SO-|P. (95) 
1 1-2€ 
Considering 2<S6C<S1L 5º <0 and remembering the | identities 


JA SW, |? = Av? |A-SO, ||? = | AO]? in addition to using the item (ij) of this 
lemma the proof of the item (iii) of this lemma is finished. 

The main result of this subsection is the following theorem 

Theorem 11 The semigroup S(?) = eº!' associated to the system (15) (17) and (21)- 
(22) is analytic when the three parameters T, o and E vary in the interval [É 51. 


1 
2 


FIG. 01: Region R,, of Analyticity de S(t) = e*!! 


Proof: We will prove this theorem will be proved using the Theorem 8, so the two 
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conditions (88)and (89) must be proved. From Lemma 9, the condition (88) is verified, it 
remains to verify the condition (89), note that proving this condition is equivalent to: let O 
>0 there exists a constant C, > O such that the solutions of (15)-(21) for IA > à satisfy the 
inequality 


AU II < CollFIaU Ia. (96) 


tis not difficult to see that this inequality (96) arises from the inequalities of the 
Lemmas 4 and10, so the proof of this theorem is finished. 


3.2 Regularity of the second system 


Next, two fundamental Lemmas for the determination of the Gevrey class and the 
proof of the analyticity of the semigroup S(t) = e? are presented. 
Lemma 12 Let S(t) = e”! be a C,-semigroup on contractions on Hilbert space H=H,, 
the solutions of the system (18) (22) satisfy the inequality 
lim sup |(1AT — Bo) "Hle(n) < 00. (97) 
|A|— 00 
Proof: To show the (33) inequality, it suffices to show that, given ô > 0 there exists 
a constant C, > O such that the solutions of the system (18)-(22) for IAI > ô satisfy the 
inequality 


JU lx 
HP Il 


< Os [UI = IICA — Bo)! FI, < CollFIn|Uln- (98) 


IF =(f1,f2,fº,f?,fs,f9) E Hthen the solution U=(d, &, y, IP, 0, 0) E D(B,) the 
resolvent equation (iA/ - B)U = F can be written in the form 


no-6 = fl in D(A?), (99) 
no- Elo tWú)+r to, +PAaO = fin D(A?), (100) 
P1 P1 Pa 
Nv-v = fin D(A?), (101) 
b ; E 
nus Av+ (dr +y)- Lora = fin D(A?), (1092) 
P2 P2 P2 P2 
iNN-0 = f in D(A?) (103) 
a) 3 Ê et . 
no+-40+ (6, +W)+ 40 = fim D(A). (104) 
P3 Ps P3 


Using the fact that the operator is dissipative B,, result in 
mi |AF 62 + poll AFWI|2 + y] AO]? = Re((i— B9)U,U) = Re(F,U) < |Fllu]Ullu. (105) 


On the other hand, performing the duality product of (100) for p,9, and as the 
operators A”for all vE R are self-adjoint, lead to 


(dr +), 60) = pB + pi(D, 1) im || AZ OP + mA PAO) + pilf2, O) + (O, do), 
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now performing the duality product of (99) for p,ty, and as the operators A” for all v 
E R are self-adjoint, result in 


2 = polfvIpP + par, 1º) 
+ u(O, 4) — po] AF]? + po(AF 3, AT) + polft ay). 


(de +), 4) + DIAZ 


Adding the last 2 equations, result in 


ellóz + VIP + db] AZU]? = pillBI2 + palv|2 — Mn || AF ul)? + po] AFI?) 
+ pi(D, SD) + polr, 18) + (AF FI ATO) + polAZ PS, ATO) 
+ plf,4) + pal!) + O, dae +). (106) 


Taking real part, applying Cauchy-Schwarz and Young inequalities, for > 0 exists 
C. >0 such that 


e? + Cs 


ellôs + W]2 + dIAZ W]? < pad]? + palv|2 + ellós + 4]? + Co FlfaU la. 


From estimative (102) and the fact O < 3, 0 ss and 0 < 5 the continuous embedding 


D(Aº2) > D(Aº'), 6, >6,, resultin 
nllde + VIP + bJAZV|P < ColFIalIU Ia: (107) 


On the other hand, performing the duality product of (101) for p,9, and remembering 
that the operators A for all v E R are self-adjoint and from (100), lead to 


o] 429]? = pa]? + pa(9, Fº) + 1(D, 0.) — n(7,0) — iv] AÉO|P + (42 15, 20) + pa($º,0). 
Taking the real part, using the inequalities Cauchy-Schwarz and Young, for e >0 
exists C, >0 which does not depend on À such that 
oI42o|? < cqllel? + Ile|ll$|] + 142 FS] A26]] + fo!) 
+ CAIS? + INPI) + eflloj? + 11429]). 
Then, from estimative (102) and the fact O < 3,0 < 5 and O < 5 < > the continuous 
embedding D(Aº2) — D(4º!),0,>0,, resultin 
Brô]|A2O]P < Co]lF la Ulha- (108) 


Finally, ftom estimates (102), (104) and (105), the proof of this lemma is finished. 
Lemma 13 Letô >0. There exists a constant C, >0 such that the solutions of (18)- 
(22) forlAl > 0 satisfy the inequalities 


[Alfa IPI? + pa]? + CIO|?) + CollPIulUla. (109) 
CIaIO]? + [19]2) + Cs|lFIa UI. (110) 


(à) Mlxllde + 4]2 + 6] 42]2] 


Es 
(ii) OMJIAZO]|2  < 


Proof: Item (1): Performing the duality product of (100) for o Ag, and remembering 
that the operators A” for all ve R are self-adjunct, lead to 
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eM(dr + UU), dr) = paIB|2 + pi(AD, Ff!) — im || AF O]? 
—inBuATO, ID +ipilf2,D) +ipilf2, SN) + ulO, Adr). 
Now performing the duality product of (99) for p,AWy, and remembering that the 
operators A” for all v E R are self-adjunct, resulting in 


A((dr +10), 1) + DIAZ 02] = pa AE]? + po (AP, 13) — ipa || AFI — ipa (AT, 1º) 
+ipa(F8 0) + ipalf!, 1º) + ulO, Ay). 


Adding the last 2 equations, result in 


Alsilóz + 4] + d] AFWIP] = Apr [BI + po]P|]2] — ifp || AF OI)? 
+ ua] AF WI) + pr(AD, Ff!) + po(AP, 8) — ifua( ATO, Fº) 
+ po( AD, PH ipilf?,D) + ipalft, VD) + ipulf2, F!) 


E A 
+ ipa! 1º) + uv IAO, = (de +). (111) 


val 


On the other hand, from (100) and (102), result in 


puAD, 1) + paço, 1º) = ifn((do +10), f7) — (O, 17) 
+ n(AT O, 1) — pr(f2, 1) + (AZ, AS FS) + nl(do +), Fº) 
— n(O, 1º) + ma( AC, 1º) — palfº, Fº). (112) 


Using the identity (112) in the (111) equation and simplifying, leads to 
Aleillós + W]2 + dIAZ 02] = Moi BI + pa]PI?) — if || AF O] + po] AF W]2) + in(óo, f2) (113) 
+in(y, 5) — iulO, f)) + ib( AZ, AZ PS) + inlóa, Fº) + inf, 1º) 
-in(O, 1) + ipi($2,D) + ipolf" o) + u(v/DIO, Ea +). 


vH 


Asfor e >0, exists C, >0, such that 


A, ' 
nv Io, és +4))| < CeAIJOIP + elMIlóz + VIP, 
taking the real part of the equation (113), applying the Cauchy-Schwarz and Young 
inequalities and estimative (98) of Lemma 12, the proof of item (i) this lemma is completed. 
On the other hand, performing the duality product of (104) for p,ÃO, and remembering 
that the operators A” for all v E R are self-adjunct, result in 


ôA 430]? = paA|JO]2 + in(D, O) — in] AÍO]2 + ips(/º,0) — in(V,0). (114) 
As 
i(D, O) = (DÃO, — fº) = ul |O, Im —iuld, fº). (115) 


Then, using (115) in (114) and taking real part, lead to 
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: 2 2 À 5 E TIRO 
ô|AI||430]]2 = pao? + u(v|AD, 0,5) — inld, 2) + ipa(fº,0) — in(P,O), (116) 


VD 


applying Cauchy-Schwarz and Young inequalities, for e > 0 exists C, >0, such that 
ôjalIl429]]2 < pyJalle|2 + ejxl||420]]2 + Co |||? + CEGA? FS] + | fSIINO] + [9/2 + |j9]]?). 


Then, from Lemma 12, the proof of item (ii) of this lemma is finished. 


3.2.1 Gevrey class of the second system 


Theorem 14 The semigroup S(t) = e?! associated to system (18)-(22) is of Gevrey 
class v for every s > E forr=mintiiz» Ti: tizle (0, 1), with (T,0, E) EE Rçç= (0, 1P. 
Proof: From the resolvent equation F = (/A/B,)U for ER, resultin U = (iAI-B,)"'F. 


Furthermore to show (57) this is theorem(6) it is enough to show: 


AF IUI <CotllFIulUln) for re(0,1). (117) 


T o E 


o & 
l4+7º l+oº 1+€ 


RR A 

Next, IATE |[6]], |A[TÊ= |V] and |»/7*z |[O |will be estimated. 

Let's start by estimating the term /À| Tr ||: tis assumed that IN >1. Set p= 
&, + 4, where À, E D(A) and &, E D(4º), with 


id + Ad = f2, Ad = -DAd6 + Edo O — FLAG + Ad. (118) 
p1 p1 p1 P1 


Firstly, applying the product duality the first equation in (118) by à,, then by AQ,, and 
recalling that the operator Ais self-adjoint, resulting in 


MB 2 + | 426,2 = (92,0,). (119) 


Taking first the imaginary part of (119) and in the sequence, the real part and applying 
Cauchy-Schwarz inequality, results in 


DIB < CIP [4394] < C|Fla- (120) 
Applying the 4- 5 operator on the second equation of (118), lead to 


AA td, = TAió + T Adv a Ate, Daio + Ato, 
p1 p1 p1 p1 


then, as [A-2w,|2 = ||? and 7-3 <5 taking into account the continuous 
embedding D(4º2) >D(Aº!), 6, >6,, resultin 


DPIIATS da]? < CHI A 9]? + AF GIP + |U]? + |6]P) + [14264]? (121) 
Using (98) and second estimative of equation (120), resultin 


4725]? < CINE UI + FIO (122) 
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On the other hand, from &,= 4 - 4, (105) andas 5 < > the inequality (120), resultin 
IAM]? < C(AFO+ AO) < CEF Ux + IFIZO: (123) 

Now, by Lions' interpolations inequality O E [-3. 3], result in 
Ba]? < C(||A72 Da]]2) 7 (147 do]2) 1. (124) 


Then, using (122) and (123) in (124), resultin 


1Da||2 < CIALTES (FIO e + FIO) (125) 


Therefore, as || 6]? = || || +|| ||? from inequality (120), (125) and IN? <IM +E, 
we have 


-2 


TB] < Cola + DIS < CAIU Ia+FIA) for 0<7<1. (126) 


On the other hand, let's now estimate the missing term AI += |W ||: Set W =W, 
+ UU, where W E D(A) and W, E D(4º), with 


b a 
NV, + Ab, = ft INV, Ay — Ed, - Ep+ Lg - Pam Av. (127) 
P2 Pa P2 P2 P2 


Firstly, applying the product duality on the first equation in (127) by WW,, the by AW, 
and recalling that the operator Ais self-adjoint, resulting in 


AP |? + AU? = (18,09). (128) 


Taking first the imaginary part of (128) and in the sequence, the real part and applying 
Cauchy-Schwarz inequality, lead to 


AP <CIFla and [AW||<C]Flyn. (129) 


Itarises from the second equation in (127) that 


t v U 2 
DA ty, =-— Aty- LA to, - LA tp Lado PA dys Aty,, 
P2 P2 P2 Pa Pa 


then, as [|4-> 4,|º= 6]; |4-> 9,]|*= ||0]|? and o- 5 <5 taking into account 


the continuous embedding D(Aº2 ) — D(Aº!),0,>6,, resultin 
PIAS Val? < CHIJoIP + IAZU]P + 4-6] + AF UIP) + (1420, |P (130) 
Using estimative (98) of Lemma 12 and estimative (129), yields 
42 Wa |? < CIPA U Ia + Flo) (131) 


On the other hand, from W, = W — W (105) and as O < 4 < >, the inequality (129), 


result in 


AF WaP < CqpASv|?+ ASP < CHPIaU la + FIG) (132) 
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Now using interpolation inequality O E [- És 1. Since 


1 o . o 
0=n( 5) +q Da for = 
using (131) and (132) we get that 
[PalP < Co(AP a P(A? Wo)! < [AS QEU] + FIO. (133) 


Also, as |W||2 = cw |2+ ||ty ||2) from inequality (129), estimative (133) and as 
INI? < INT$S, resultin 


) 


ANTE E] < Cola + IMJVIÊ < CONT (FI UIa + FI) for 0<0<1 (134) 


£ 


Finally, let's now estimate the missing term IAI 1% | 
O,€E D(A) and O, E D(Aº9), with 


Ol. SettO =0,+0,, where 


JO, + 40, = fº AO, = A9 —- L4, —- Ly - TAB + 46, (135) 
P3 P3 P3 P3 


Firstly, applying the product duality the first equation in (135) by O,, and as the 
operators A" for all v E R are self-adjoint, result in 
MO]? + [| 4304]? = (18,01). (136) 
Taking first the imaginary part of (136) and in the sequence, the real part and applying 
Cauchy-Schwarz inequality, lead to 
MISi<ClFla and 4%] <C]Fla (137) 
It arises from the second equation in (135) that 


DA to, =-2Ato- Lato, - Haty- Tattoo Ato, 


P3 P3 P3 P3 


then, as [4-3 4,.|2= ||], =; <O and E- < 5 taking into account the continuous 
embedding D(4º2) >D(4º'), 6, > 60, and using estimative (98) and second estimative of 
equation (137), lead to 


472 6a]? < CINE] llU e + LF ao). IR) 


On the other hand, ffom O,=O0 — O,, (105) and as el the inequality (137), result in 


lsfes|P < crl4éol? + AO) < CEF U|y + FB) (139) 
Now we are going to use an interpolation inequality O E [- ã , a. Since 
O) É 1 6 for = 
= 3 +( Do [0 1= Tê 


using (138) and (139) we get that 
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leal? < co(|A-36a PAPO] < [AI (UFI UI + FIO (140) 


Also as, |o ||? < ctllo,||º + |0,1|3 from first inequality of (137), (140) and as IN 
< INTÉ, result in 


£ 2€ 


AEIO<ClFlx <> INOJê< CATE QEIaUIa+ FIGO for 0<E<1 (141) 
Finally, using the estimates (126), (134) and (141) in the inequality item (i) of Lemma 
13, resultin 


2 1 2 —27 —20 —2€ E 
Bylsllde + WI]? + dIAZ?) < CLA + JApTES + JA[TEE HIFI |U a 
+ COMI Ia lU In. (142) 


And, using the estimates (126) and (141) in the inequality item (ii) of Lemma 18, 
result in 


27 


o] AZO|? < Cf + ATÉ 


Fal Ia + Co JAP Ie UI. (143) 


T a 


As for (T, 0, 8) E [0, 1), we have —1 2 max(-55., sto. —rie)= 2 max (15185 ré) 
forlA>1. 
Finally, summing the the estimates (126), (134), (141), (142) and (143), we obtain 


|jU!2, < Co|Aj2mintrEsvtfs The) [I F F 


nUlu + Fla). 


Equivalently: 
[De To k 


Ullu < CollFlly for (7,0,€) e [0,1]º. 


As r= mini 15; rc) E (0, 1) for (T, 0, E) E (0, 1]º the proof of this theorem is 
finished. 
3.22 Analyticity of S(t) = e?! for (T, o, E) [>, 1]º such that T = E. Fig. 02 


Before proving the main result of this section, the following lemma will be proved. 
Lemma 15 Letô >0. There exists a constant C, >0 such that the solutions of (18)- 
(22) forlAl > 0 satisfy the inequality 


. : 2 IL O! r : l 
(à) IMISIP+]O]<OIFIalUlu for s<T=E<1 (144) 


2 ' , 1 - 
(ii) [MVP <CFlullUlu for S<o<l (145) 


p 
Proof: Item (i): Realizing the duality product of (100) with p AT AQ and using the 
property that the operator A“ is self-adjoint for all v E R and using the equation (99), result in 
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ABI = “(M-AG+ Ur), ATO) + UPA TTO) ci RIA TOP (Ou, MATA) 


= iD(AD+ AS!),ATO) =iS(U, + f3,470) — RIA Ed]? 


1 
+P,inAlTTO — RAT br + ind — ip AT fo) + (O, AMT Ba). 


Then 


ISA GIP + (ABS! AÉTO) di (0, AB) = UE ATO) 


Ed 
ai GAcÊ RIP (146) 


HP) UP, ALTO) AIM TP 
x 4-39]2 + É -(8, JA D,). 


On the other hand, now performing the duality product of (104) for É “A “AO, using 


(101), and as the operators A“ for allv E R are self-adjoint, result in 


e du o PS £ “ ô I-E wo o) 1 «rm fo) . £ 4 
ael?= Ate? +i-l4 ep + (Ai, AÉ-tO) io (fº, Al-E9) — É (19, A-0,) 


=D (9º, AEW) — i(f8,0) +if8 71 Aê fº|2 — — A Ed, 6). (147) 
F 


Imposing the condition T = & and since Re (EO, AA” 9,E-AA” b, OEJ=-0, adding the 
equations (146) and (147) and then taking the real part, result in 


À [Este + IP E Es ATO) +70) — (SE, ATD) (6) 
/ 1 


EP, ALTó JAIATP, va) HAS fº, Ate) 


ER Al-S9) — É (SS, AEB) — E (SS, ATE) —i($º,0) 
From lis(W,ATW,)|<C(v|B+|A E O2)and É <T=E<1, applying Cauchy- 
Schwarz and Young inequalities, norms ||F ||, and |U ||, the proof of item (j) of this lemma 


is finished. 
RI : : 2a : 
Item (ii): Similarly, performing the duality product of (102) for poA AU, using (101), 


and recalling the self-adjointness of 4”, v E R, lead to 


+ (At f8, AS) — E (9,4""W,) 


; q & o ab 
Agp? = RAT VI +i 


nm Dia ae patup si UA -ep) 


PRA A . b ; 1- KR 4 é 
(VDlo, ==A"") =i—(f4,Al-0p) if, 464) 

rr ZÉ o) Ho 

il ACOU) + GUI SAO) “UV + IP. 


Noting that: For 3 < ae. > -0<0and = <0,onthe other 
hand as [A W,||2= ||A 12º W||2, taking real part and considering that 5 < o < 1 and using 


Cauchy-Schwarz and Young inequalities, norms ||F|), and |U||,, for=>0, exists C>0 
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independent of À such that 


MIO < ColPIalUlx + CAD]? + |V?) + Ce 


AJjej? + ejaphy. 


Using estimative (98) and estimate (i) of this lemma, the proof of item (ii) of this 
lemma is finished. 

The main result of this subsection is the following theorem 

Theorem 16 The semigroup S(t) = e*2 associated to the system (18) (22) is analytic 
when the three parameters tT, O and E vary in the interval [ : MW witht = E: 


É 
Ss 


FIG. 02: Region R,, of Analyticity de S(t) = e? 


Proof: This theorem will be proved again using now the Theorem 8, so the two 
conditions (88) and (89) must be proved. 

Here the condition's test (88) will be omitted because it is completely similar to the 
test already conducted for the first the first system. 

Next, the condition (89) is proved, note that proving this condition is equivalent to 
show, letô >0. There exists a constant C, >0 such that the solutions of (18)-(20) and (21)- 
(22) for IAl> ô satisfy the inequality 

UI < CollFIalU Ia. (148) 

Itis not difficult to see that this inequality (96) follows from the inequalities of the 
Lemmas 13 and 15, so the proof of this theorem is finished. 

Remark 17 (Asymptotic Behavior) lt is emphasized that thanks to the Lemmas 
3 and 9 the exponential decay of the first system is obtained when the 3 parameters 
take values in the closed interval [0, 1], also for the second system using the Lemmas 
12 and proving in a completely similar way to the proof of Lemma 9 which iR cp( B,) the 
corresponding semigroup S(t) = e?2, will be exponentially stable when the 3 parameters 
take values in the closed interval[0, 1]. 
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ABSTRACT: In this work, we study the 
stability and regularity of the system formed 
by the third-order vibration equation in 
Moore-Gilson-Thompson time coupled with 
the classical heat equation with Fourier's 
law. We consider fractional couplings. He 
then fractional coupling is given by: nAºo, 
anAºBu, and nAºBu, where the operator Aº is 
self-adjoint and strictly positive in a complex 
Hilbert space H and the parameter p can 
vary between O and 1. When & = 1 we have 
the MGT-Fourier physical model, previously 
investigated, see; 2013[1] and 2022/9], 
in these works, the authors respectively 
showed that the semigroup S(t) = ef 
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associated with the MGT-Fourier model 
are exponentially stable and analytical. The 
model abstract of this research is given by: 
(3)-(5), we show directly that the semigroup 
S(t) is exponentially stable for p E [0, 1], we 
also show that for q = 1, S(t) is analytic and 
study of the Gevrey classes of S(t). We show 
that for p e ( 5 1) there are two families of 
Gevrey classes: st > 2 when q € (1/2, 2/3] 
and s2 >= When pEl2/3, 1), inthe lastpart 
of our investigation using spectral analysis 
we tackled the study of the nonanalyticity 
and lack of Gevrey classes of S(t) when q 
e [0, 1/2]. For the study of the existence, 
stability, and regularity, semigroup theory 
is used together with the techniques of the 
frequency domain, multipliers, and spectral 
analysis of a system, using the property of 
the fractional operator Aº for p E TO, 1]. 
KEYWORDS: Asymptotic behavior, 
Stability,  Regularity,  Gevrey Class, 
Analyticity, MGT-Fourier System. 


11 INTRODUCTION 


Various researchers year after 
year have been devoting their attention 
to the study of asymptotic behavior; the 


model Moore-Gilson-Thompson equation 
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appearing in the context of acoustic wave propagation in viscous thermally relaxing fluids, 
in [14], they studied the model abstract: 


Tum + QU + C2AU + bAuU, 


where A be a self-adjoint positive operator on H with a dense domain D(A) C H, 
show exponential stability requires y = a — E > 0. In the complementary region of the 
parameters, the system is unstable (y < 0) or marginally stable (y = 0). Quite interestingly, it 
can be used as a model for the vibrations in a standard linear viscoelastic solid, for it can be 
obtained by differentiating in time the equation of viscoelasticity with an exponential kernel 
(see [8] for more details). 

In [19], studies the decay of the energy of the Moore-Gibson-Thompson (MGT) 
equation with a viscoelastic term, this work is a generalization of the previous one (Part 1) 
see [20], in the sense that it allows the kernel memory to be more general and shows that 
power decays the same way kernel memory does, exponentially or not. The model is given 


by 


t 
Tum + Qua + 2 Au + bAu, — / g(t — s)Au(s)ds = 0, (1) 


0 


where A is a positive seft-adjoint operator defined in a real Hilbert space H. The 
convolution term [. g(t — s) Au(s)ds reflects the memory effect of viscoelastic materials; the 
“memory kernel” g(t): [0,0) — [0,º) directly relates to whether or how the energy decays. 
Without this memory term, it is known the MGT equation has exponential energy decay in 
the non-critical regime, where 7 =“ E >0 see [14]. 

One of the first investigations of the coupled MGT-Fourier model, was published in 
2013, see Alves et al. [1] studied the exponential decay of the MGT-Fourier physical model 
given by 


( aus + Um — a2 Au — a28Au, + nÃO = 0. 
0, — AO — anus — nÃu, = 0. 

where x € O, t E (0,90). The function u = u(x, t) represents the vibration of flexible 
structures, and 6 = 9(x, t) is the difference of temperature between the actual state and a 
reference temperature. The constants a and B are positive, n is the coupling constant. O is 
a bounded open connected set in Rn(n > 1) having a smooth boundary T = 90. The initial 
conditions are given by (5) and the boundary conditions:u = 0, ut=0,0=0,on dO =T. 

In more recent research from 2022 [9] they studied the model abstract MGT- Fourier 
system: 


um + aum + BAºu + yAºu = n40, 
0, + nÃO = —n Aug — na Au, 
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in the subcritical regime B >, where the operator A is strictly positive selfadjoint. For 
any fixed p € [1, 2] the authors showed that the solution associated to the semigroup S(t) 
is analytically and exponentially stable as well. Also in 2022 M. Conti et al [5], studied the 
model: 


um + aum + BAZu, + yA2u = —nÃO (2) 
0, — nÃO = nÃu, + andAu,. 


His research is focused on analyzing the energy transfer between the two equations, 
mainly when the first one is in the supercritical regime and exhibits an antidissipative 
character. The leading actor becomes the coupling constant n, governing the competition 
between Fourier damping and MGT antidamping. In fact, they showed that a large enough 
n is always capable of stabilizing the system exponentially fast. One of the characteristics 
of this model is the presence of the bilaplacian in the first equation. With respect to the 
analogous model with the Laplacian, this forced to adjust the mathematical techniques. On 
the one hand, the energy estimation method does not seem to be applied directly, on the 
other hand, there is a gain in regularity that allows us to rely on analytical techniques for 
the properties of the semigroup associated with the model. In light of the above discussion, 
the natural question to be addressed is how the dissipation produced by the heat equation 
influences the asymptotic dynamics of the system. In the subcritical case, both equations 
in (2) are dissipative (actually exponentially stable) and good stabilization properties are 
expected. This idea was confirmed in [9], where itis shown thatifu = y — aB < O then the 
semigroup associated to (2) is exponentially stable and also analytic. This means that the 
coupling allows not only a dissipation transfer between the equations, but also a regularity 
transfer. This research deals with the behavior of the system in the critical and supercritical 
regimes, that is, when yu > O. In these cases, the system (2) consists of a conservative 
system (if u = 0) or antidissipative (if yu > 0). More recent research in this direction can be 
found at [4, 3, 11, 23]. 

In the last decade, the number of investigations focused on the regularity of the 
semigroups S(t) associated with coupled systems has increased considerably. The interestis 
centered on Differentiability, Gevrey class, and/or Analyticity, among the models considered 
systems with fractional couplings, different types of dissipation (internal, localized), and 
fractional dissipations. One of the motivations for this new research is that the analytic 
semigroups associated with linear systems imply that this type of semigroup admits Gevrey 
classes, which in turn implies the semigroup's differentiability. And depending on the 
regularity, the semigroup will have the best asymptotic behavior (Exponential Decay). More 
recent research in this direction can be found at [2, 6, 7, 13, 16, 17, 18, 24, 25, 26, 28, 29]. 

The paper is organized as follows. In section 2, we study the well-posedness of the 
system (3)-(5) through semigroup theory. We leave our main contributions for the last two 
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sections; in the third we dedicate to stability, and in the fourth, to regularity. in the section (3) 
we showed that the semigroup S(t) = is exponentially stable as q E TO, 1]. In the last section 
(4) of the regularity, we studied the analyticity when q = 1, determine two family of Gevrey 


classes: st >2 when q E (5 ; 5) and s2 > e when pel é 1). And when the parameter & 
assumes values in the interval [0, 5] it shows that S(t) does not admit Gevrey classes and 
is not analytic. We prioritize direct proofs using semigroup-theory characterization together 
with frequency domain methods, multiplier, and spectral analysis of the fractional operator 
Agp for pelo, 1]. 


21 WELL-POSEDNESS OF THE SYSTEM 


In this section, we will use the semigroup theory to assure the existence and 
uniqueness of strong solutions for the system (3)-(5) where the operator A is defined to 
follow, let A: D(A) CH — H be a strictly positive selfadjoint operator with compact inverse on 
a complex Hilbert space H, with a continuous embedded domain D(A),-> H. Consider the 
abstract linear system: 


QU + Um + a2 Au + a2B Au, — nÃºO = 0, (3) 
0, + AO + anAºu, + nÃºu, = 0. (4) 


The initial conditions are given by 

u(x, 0) = uO(x), ut(x, 0) = u1(x), utt(x, 0) = uZ(x), 9(x, 0) = 00(x), in O. (5) 

Remark 1 Itis known that this operator A is strictly positive, selfadjoint, has a compact 
inverse, and has compact resolvent. And the operator Ar is self-adjoint positive for allt E R, 
bounded for T < 0, and the embedding 


DA) > DAM), 


is continuous for 71 > T2. Here, the norm in D(A?) is given by HullD (ao:= HATull, u E 
D(A?), where (:, -( and Il - Il denotes the inner product and norm in the complex Hilbert space 
H = D(Aº). Some of the most used spaces at work are D(A?) and (A?). 
Taking v = aut + u, the system (3)-(4) can be rewritten as 
va + a2 Av — a2(a — B) Au, — nÃº0 = 0, (6) 
0, + AO + nÃºv, = 0. (7) 


This system (6)-(7), for Pp = 1 is the model of a system of the coupled viscoelastic 
equation coupled with the heat equation given by Fourier's law. Taking the duality product 
between equation(6) and vt, and (7) with 0, taking advantage of the self-adjointness of the 
powers of the operator A and using the identity v = aut + u, Let B - a > 0, for every solution 
of the system (3)-(5) the total energy E: R* > R+is given in the t by 
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1 E “ , 1 1 ç € 
E(t) = 3 law + 0]? + aZa(8 — a)JjA?v]? + 2a 420 + AZul2 + |Jo]2| (8) 


and satisfies 


Set) = —0?(8 — a)Atol? — [1 4tojP (9) 


Taking ut = v and vt = w, the initial boundary value problem (3)-(5) can be reduced to 
the following abstract initial value problem for a first-order evolution equation 


d cons 
aU (b) =BU(t), U(0) = Us, (10) 


where U(t) = (u, v, w, 8, UO = (uO, u1, u2, 00 and the operator B: D(B) CH — 


T 
BU := (v. w, E [a? Au + a28Av — nAº0 + w] — [nAºv + anAw + A6)) ; (11) 


a 


forU=(u,v w 0 forU=(u,v, w, 6). This operator will be defined in a suitable 
subspace of the phase space 


H := [D(A%P x [D(AP, 
itis a Hilbert space with the inner product 
(U,,Uo)ja := aZa(B — o)(AZw, Ava) + a2(A2uy + aAiv, A? uz + aA2v,) 
+(v1 + aus, vz + awa) + (01, 05). 
for Ui = (ui, vi, wi, Oi) E H, i = 1, 2 and induced norm 
“ º ; 1 « ç 1 1 ç p p 
JUIã = 020(8 — o)JA?o]P + 2 AZu + a 420] + o + aw|2 + ol. (12) 
In these conditions, we define the domain of B as 
D(B) := (u eH: (v,w) € D(AMP, aZu + a2Bv — n4º *0 E D(A), 
nAº Iv + anAº tw + 0 E D(A) + (13) 

To show that the operator B is the generator of a CO-semigroup, we invoke a result 
from Liu-Zheng' [22]. 

Theorem 2 (see Theorem 1.2.4 in [22]) Let B be a linear operator with domain D(B) 
dense in a Hilbert space H. If B is dissipative and O E p(B), the resolvent set of B, then Bis 
the generator of a CO-semigroup of contractions on H. 

Proof: Let us see that the operator B given in (11) satisfies the conditions of this 


theorem(2). Clearly, we see that D(B) is dense in H. Taking the inner product of BU with U, 
we have 


Re(BU,U)g = —aX8 — a)|A?v|2— |A20]?<0, VU CD(B), (14) 
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and since B > a, the operator B is dissipative. 

To complete the conditions of the above theorem, it remains to show that O E p(B). 
Therefore we must show that (Ol — B)! exists and is bounded in H. We will first prove that 
(OI — Bexists, then it must be proved that B is bijective. Here we are going to affirm that B is 
surjective, then for all F = (f!, f2, fº, f)! & Hthe stationary problem BU = F has a solution for U 
=(u,v, w, 6)T E D(B). From definition of the operator Bin (11), this system can be written as: 

v=f,w=f?and 


a Au = -nAº!h-m AI - apAtIf- afS-aBAf!—f2, (15) 
AS = —f!-nAºf! — anAºf2. (16) 


From these equations, this problem can be placed in a variational formulation: v, 
w E ID(G)2 and we write the last two equations (15)-(16) in a variational form, using the 
sesquilinear form b : 


b(u, 0; 2) = (gr, 2; br, ba) (17) 
Where 


ga =-—nAº!fy—-n AF ans Iip- af-aBAf!- fgo=-—f!— 
nAº f! — anAº f? 


and the sesquilinear form b is given by 
b(u, 0; 19) = a2(A? TA AZ) + (A20, AZ»). (18) 


as D(A),—- D(A >), — D(Aº), — D(A >) — D(A'), this sesquilinear form is coercive 
in the space [D(A e. Asgl,g2 E [D(A->)I from Lax-Milgram's Lemma the variational form 
has unique solution (u, 0) E [D(AS) x D(AS) )] and it satisfies (15)-(16), such that (17) is 
verify for (1, W2) E ID(AS)P. From (18) and for all (1, y2) E ID(AD)P, we have 


(Alu, AM) + (A20, Aba) =(-f!-nAºf! — anAº f2, 1») 
o (-n 4º! fa E no AZ fl = ant A? p? ar afé = aBAf! - f2, vi), (19) 


choosing y2 = 0 E D(AS) from (19) and Vy1 E D(AS), we have 
a(AZu, Ay) = (-nA* If —- ASIA opa Ip af aBAf— Pb), 
then 


a2 Au Ee —nAº fa o nº A26-1 1 = an 421? Ee afê 


“BA in [D(AM =D(A?). (20) 


On the other hand, choosing Wi = O E D(A;) from (20) and Vy2 E D(AS), we have 


(40,12) = (—f! — nAºf! — anAº f2, 4») tnen 
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A0=-f!-nAºf!-onAºf? in [D(AM] = D(A?). (21) 


(20) and (21) are solutions in the weak sense. From F = (f,f, f,f)'CHandO<q 
2p—1 


< 1 implies that -> < 2 <:, we have D(A;),—> D(A),— D(A-=) ,—> D(A O — D(AS)), 
from (20), we have 


l Dm 4 7 2d —s 2 2-2 42 pm : 
u=—|-— na? fa —- n AA — ams? p - Ap 
a 


—aBft- ATP] in D(A?). (22) 


2p-1 2p-1 


Now ffom0<p<1timplesthaiZz2-Ppz1>— <: we have D(AS),— D(A 5) 
> D(A),—> D(A?),—= D(A?'),—> D(A''), from (21), we have 


0=-A1f! mA If! -onAP in D(A?). (23) 


Then alu + a2Bv — nAº-10 E D(A) and nÃ? Iv + and lw + 0 E D(A). 

The injectivity of B follows from the uniqueness given by the Lemma of Lax-Milgram's. 
It remains to show that is a bounded operator. From v = f!, w = f!, equations (22) and (23) 
and (12), we get 


2 à 1 —3 26-8 2 A26-B fp 

JU = aa(8-a)||Aà Fº |P+a?]a As fm [-n4º à fa Ai fan A? pf? 
—aA cê — a!gAsp' - Ap + las? + pI 

HA Pons fc anAt A PP (24) 


Using norm IIFIIH and applying inequalities Cauchy-Schwarz, Young and applying 
continuous embedding D(ArT2),— D(ArT1), T2>T1, after making various estimates, we obtain 


NU la = IB Fllã < CIlFlã- 


Therefore B' is bounded. So we come to the end of the proof of this theorem. 
As a consequence of the previous Theorem 2, we obtain 


Theorem 3 Given UO € H there exists a unique weak solution U to the problem (10) 
satisfying 


U e C([0, +00), H). 


Futhermore, if UO E D(B"), k E N, then the solution U of (10) satisfies 


A 
U e (Jc i(lo, +00), D(B). 


j=0 
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Theorem 4 (Lions' Interpolation) Leta< B <y. The there exists a constant L = L(a, 
B, y) such that 


1 4%ull < LIASuj=o | ATul|= (25) 
for every u E D(AY). 
Proof: See Theorem 5.34 [10]. 


31 STABILITY 


3.1 Exponential decay of the semigroup S(t) = e!” 


In this section, we will study the asymptotic behavior of the semigroup of the 
system (3)-(5). We will use the following spectral characterization of exponential stability of 
semigroups due to Gearhart[12](Theorem 1.3.2 book of Liu-Zheng [22]). 

Theorem 5 (see [22]) Let S(t) = e!” be a CO-semigroup of contractions on a Hilbert 
space H. Then S(t) is exponentially stable if and only if 


p(B) O (N/A ERj=iR (26) 
and 
lim sup ||(iX7 — B)!lc(m) < 00 (27) 
|A|-+00 
holds. 


Remark 6 Note that to show the condition (27) it is enough to show that: Let O > 
O. There exists a constant Cô > O such that the solutions of the system (3)-(5) for IAI > 6, 
satisfy the inequality 


JU la 
Fla 


<OG +» |]Ulk<GslFIuUla < oo. (28) 


In view this theorem(5), we will try to obtain some estimates for the solution U = (u, v, 
w, O of the system (AI -BJU =F where A E R and F = (f',f2, fº, f)T E H. This system, written 
in components, reads 


du-v = fin D(A?) (29) 
iv-w = fin D(A?) (30) 
iakw +a2Au+alBAv— nAO+w = afé in D(A?) (31) 
iMO +nAºv+anA?w+ AO = fin D(A?) (32) 
From (14), we have the first estimate 
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a(8 — a)JJA2ulP + 420]? = |a*(8 — a)]Aiv])? + |] 429]| = | — Re(BU,U)] 
|Re(iNUIP — (BU, U))| = IRef((iM — BJU,U)) 


(FU) < Ia Io. 


IA 


Therefore 


; 1 1 

"(8 — o) AZvl/ê + || 4201] < [FU Ie. (33) 

Next, we show some lemmas that will lead us to the proof of the main theorem of 
this section. 


Lemma 7 Let ô>0. There exists Cô > 0 such that the solutions of the system (3)-(5) 
forlAl>, satisfy 


(i) IMJAZuP < ColFlalUla for 0<6<1, (34) 
(ii) JAPIAZU]? < ColFluUla+IFIã) for O<6<1, (35) 
(iii) |, Au) < CaFlaUla+IIFIã) for O<ó<1. (36) 


Proof:(i): Taking the duality product between equation(29) and Au, taking advantage 
of the self-adjointness of the powers of the operator A, we obtain 


A AZul2 = (AZv, AZu) + (AZ FI, AZu). (37) 


Applying Cauchy-Schwarz and Young inequalities and norms IIFIIH and IUIIH, for 
>0, exists Ce > 0 such that 


MI 4Zu|2 < C.|/A2v]]? + e] Au]? + Cs] FlhalJU 


H- (38) 
Therefore for IA > 1, we have ellAzul)? < elAIA? ul? and using (33) finish 
proof of item (i) this lemma. 
Proof:(ii) On the other hand, multiplying by À the equation (37) and applying Cauchy- 
Schwarz and Young inequalities, for e > 0, exists Ce > O such that 
APIJAZul]? < Co Ao]? + clAPIA Zu]? + CoPi, (39) 


from (33) finish proof of item (ii) this lemma. 
Proof:(iii) Using identity ww = —AZu — iAf! — fº we have 


(w, Au) = (—Au — iXf! — f2, Au) = —AM || A?ul/2-i(A2 f!, AAZu)—(A3 f2, AZu). 


Applying Cauchy-Schwarz and item (ii) this lema, finish proof this item. 


Theorem 8 The semigroup S(t) = e is exponentially stable when the parameter q 
assumes values in the interval [0, 1]. 
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Proof: Let's first check the condition (28), we will start by proving in the following 
lemma the estimate of the term a?/|42u + a 42v||2 and Ilv + awil2of IIUI: 

Lemma 9 Let ô>0. There exists Cô > 0 such that the solutions of the system (3)-(5) 
for lAl >, satisfy 


|A%u+aA?v|2 < CãlFlsUlz for 0<ó6<1 (40) 


Proof: As ||42u + a Av]? < CIJAZul? + ||A2v|2] and for IM > 1 from item (i) of 
Lemma 7 and estimate (33), we finish the proof of this lemma. 

Lemma 10 Let ô > 0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAl >, satisfy 


Iw|Ê<ColFlalUla for 0<ó<1 (41) 


Proof: Taking the duality product between equation (30) and w and using (31), taking 
advantage of the self-adjointness of the powers of the operator A, we obtain 


2 1 . 2 a? 1 1 
rw!" = ——(v,iahw) — (fº,w) = —(AZv, AZu) 
a a 
28 1 2 1 4-1 1 < 2 
+ Ildo am (Adu, Aº-30) +>(,uw)- (0, f)— (Pow). (42) 


A A 


For e >0, exists Ce > O such that 1 al(v, w)l < Cellvll + cllwll2. Applying estimates 
(33) and (49) of Lemma 14 for lAI > 1, we finish to proof this is lemma. 

Finally the following lemma estimates the term Ilv + awll2 of IIUII2 

Lemma 11 Let ô >0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAl >, satisfy 


lo+ow|b<C;lFlsUlu for 0<6<1L (43) 


Proof: As Ilv + awll2< C(Ilvll2+ a2llwll2)2and as O < applying continuous embedding, 
estimative (33) and Lemma 10, we finish to proof this is lemma. 
Finally, using the Lemmas 9-11 and estimates (33), we get 


| 
Therefore the condition (27) for p E [0, 1] of Theorem 5 is verified. Next, we show 
the condition (26) of Theorem 5. 
Lemma 12 Let p(B) be the resolvent set of operator B. Then 


U 


É < CollF lx |U Ilza for 0O<qG<1. (44) 


iR Cc o(B). (45) 

Proof: Since B is a closed operator and D(B) has compact embedding into the energy 
space H, the spectrum o(B) contains only eigenvalues. Let us prove that iR C p(B) by using 
an argument by contradiction, so we suppose that iR/C p(B). As O E p(B) and p(B) is open, 
we consider the highest positive number AO such that the ] — iÃO, iÃO[C p(B) then iÃO or -iÃO 
is an element of the spectrum o(B). We Suppose iÃO E o(B) (if -iÃO E o(B) the proceeding is 
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similar). Then, for O < ô < ÃO there exist a sequence of real numbers (An), with O <An< AO, 
An —> O, and a vector sequence Un = (un, vn, wn, On) E D(B) with unitary norms, such that 
H(iAnl — BJUnIlH = HFnllH — 0, 


as n — e, From estimative (44), we have 


'n H aa B o) 2Um 2 a? Aun + a AZv, 
l H ] + 
vn own ||? 116 n|| < Csllf nm |U nllH Csll! nl H > 0. (46) 


Therefore, we have IIUnllH — O but this is absurd, since HUnllH = 1 foralln e N. 
Thus, iR C p(B). This completes the proof of this lemma. 

Therefore the semigroup S(t) = e is exponentially stable for p E [O, 1], thus we finish 
the proof of this Theorem 8. 


41 REGULARITY 


4.1 Study of Analyticity of the semigroup S(t) = e! 
In this subsection, we will show that the semigroup S(t), is analytic for the parameter 
& = 1. The following theorem characterizes the analyticity of the semigroups S(t). 
Theorem 13 (see [22]) Let S(t) = e!” be CO-semigroup of contractions on a Hilbert 
space H. Suppose that 
p(B) D (iN/AERJ=iR. (47) 
Then S(t) is analytic if and only if 
lim sup |JA(iA — B)lcq < 00, (48) 
|A|-+00 
holds. 
Before proving the main theorem of this subsection, we will prove some lemmas that 
will be used for analyticity for p = 1 and the determination of Gevrey classes para & E (G 1). 
Lemma 14 Let ô > 0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAl > 5, satisfy 
(i) IAIAZuIP 
E ç 24 
(ii) JAPIA 


CollFlaUla for 0<ó<l, (49) 
ColFlalUla + Fla) for 0<ó<1. (50) 


a 
To < 
Proof: (i) Taking the duality product between equation (30) and Aº v, we obtain 


24—1 24 


Aiv|P = (AT w, A2v) + (AS PAT 0). 


iÀ 


Applying Cauchy-Schwarz and Young inequalities and estimative (51) of Lemma 15, 
we finish the proof of item (i) of this lemma. 


20-1 


Proof: (ii) Applying the operator 4 2 in the equation (30) and then using Young 
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49 291 241 
7 vlÊ=CsNAT wP+ AT PP). 


Using item (i) of the Lemma 15, we finish the proof this item. 


Lemma 15 Let ô > 0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAl >, satisfy 


(DIA up < GlFlalUlz for 0<ó<1, (51) 
(ii) IXPIJAS o]? < Gs a +) for 0<ó<7 (52) 
(iii) MIA w]P < CglFlaUla for 0<6<1, (53) 
(10) 1006, 41-60) < ColPlaUla for 5<6<1 (54) 


Proof: (i): Performing the duality product between equation (32) and Aº 'w, taking 
advantage of the self-adjointness of the powers of the operator A, we obtain 


Tu) + (1,4) 


(55) 


“ul? = =(4* 'B,iadkw) — n(APT 0, APT w) — (A3O, 


O AO Ada) -2O(ASE O, Adv) + 
x 


—-Lg-to, AT w) + (Aº-g, f3) — MAE v, AE w) 
a 
w) + (f4, A). 


AppIying Young inequalities and from -:i<Z <: and q — 1 <0, using continuous 


26 1 


—(A&O,A 


embedding, for e > 0, exists Ce > O such that 


+ 4292) 
“E wlfê + Conf + CIfIIho]. (56) 


T ul <CHAZu]? + | Atul)? 


For lAl> 1, we have IAGul2s INIAGulE, from Lemma 7, estimative (33), using norms 
HFIIH and IUIIH finish proof of item (i) this lemma. 

Proof (ii): On the other hand, performing the duality product between equation (31) 
and Aw, taking advantage of the self-adjointness of the powers of the operator A, we obtain 
iolAZI| A w|? = —a 24 Tu, AS E) — aí 28(RAFE “u, A E y) — AJA? aw||? 

+n(DO, E w) + FS, AAAÉ-D ea 
= a (AF 7 “) +ia 28 AH Tu 2 + a(f3, A A2$-Dap) 
+ia2 (A? 2, A22w) + (AO, A8*2w) — AJ AS ml, 


then, taking imaginary part, for e > 0 exists Ce > 0, such that 


IAPIA?w]P < Gs 


T ul + AE w]2 + || 42 PI 42-00] (57) 
aii 4% -âw|PP + Co FP + el 26). 


Therefore, as 0 < 4 < i= ps Micic 5 f item (ii) Lemma 7, item (i) this lemma 


and estimative (33), finish e this a 
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Proof (iii): Multiplying by Aº'w the equation (31) and by nA?º2w the equation (32) 
and adding the results, we have 
iaA AT w||2 + en ido — f2) + a28(Aºv,ikv — f2) — n(Aº0, A Im) 
+ IA ul + in(ÃO, 222) + nº (Aºv,iXAZ2y — (202 2) 
+ am PIA” 7 “Pulp? + n(4º0, Au) = a(f8, Au) + n(f!, A). 


Taking imaginary part, we have 


aMIA“ wl2 = Rea + PASS v||2 + Imfa25( AFR uy, A f2) 


+ ia? (VlAZu, 
a 


Ra VIA 422) + a(f8, Aºw) + m(f1, Am). (58) 
vHI 
30-2 201 0 
25 


Aiv) + (APT u, AZ f2) + m( ASR, Aà f2) 
+ in( 


of (58), asfor0Ospsiwehave >< <>< 


embedding and aplying Cauchy-Schwarz and Young inequaliies for e>0 exists Ce >0, 


= > and 3p — — <> > Using continuous 


such that 
gel ç 6 E. if 
AJA vw] < CotIFIhalU Ia + |A/]|420]]2 + ]AI] 42 u]/P) 
+ C.I9||2 + ejA||| A2*-2w]|? for 0O<ó<1. 


Finally, as 2p — 2 < — using continuous embedding and item (i) Lemma 7 and item 


(i) of Lemma 14, finish proof this item. 


Proof (iv): Just observe that for > E p, we have > >— p< -- < - Then 


(19, Al-v)| = (420, A2*(-iw — if2))| < |(420, 42%w)| + |(420, AZ [2]. 
Lemma 16 Let ô > 0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAI > 5, satisfy 
(i) MOI? <CollPIallUlx for 0<ó<1, (59) 
(ii) al Bida < CotllFlalUla + Fla) for 0<ó 


E 


(iii) IAPIJAS 


< CotllPlaUla + Fla) for 0<ó<1 (61) 


Proof: (i): Performing the duality product between equation (32) and À 0, taking 
advantage of the self-adjointness of the powers of the operator A, we obtain 


aljojf? 


Il 


—iM||A-20]P — n(=iw — if2, A 19) — n(adw, ATO) + (f! MATO) 
= iMIjAZOI + in(w, A?19) + in(f2, Aº19) — ina?(AZu, A? 20) 
—ina2B(AZv, Aº20) + ipa 9 — intw, Aº*0) + ina(fê, AO) 
+in(f4, AS to) + ian( ft, AS Iwo) + (94,0) — JA Ê FIIP. (62) 


Taking the real part of (62) and simplifying, we get 
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Ajoj? = Refin(f2, Aº19) — ina2(AZu, Aº-20) — ina2B(A?v, A?-30) 
+ina(fº, 10) + im(f!, Av) + ian( ft, Aºw) + i(f',0)). (63) 

As |Ref-ina?(AZu, Aº-20))| < C4|A?ul)? + ||4º-20]2). And forlal> 1 ando < 4 <1 
we have, P-1<0and q - -<5 applying continuous embedding D(A!?),—> D(A!'), T2>T1, 
estimative (33) and applying Cauchy-Schwarz, we finish the proof of item (i) of this lemma. 

Proof of (ii): Multiplying by À the equation (63), we have 


AB ||9]]2 = Refin( A?! 2,10) — ima2(AA2u, Aº-26) — na? B(ikv, Aº0) + i(f*,A0) 
+ina(f*, NAO) + in(ft, AA Io) + ian( fi, AA). (64) 


For e >0, exists Ce >0, such that 


Retin( 4º! 2, A9))| < ejAPIo|]P + CojpAt PP, (65) 
|Re(-ina2(XA?u, Aº-20)| < Co(IAP|AZ ul)? + || Aº-20]]2), (66) 


on the other hand, from equation (30), we have iAv = w +f? and using Lemma 15, we 
have 


26-1 


na B|Re((AT w, 20) + (A? f2, Aº-20))| 
CF IalU la + [FÃS (67) 


Ref-na?B(iXv, Aº0))| 


I| 


IA 


besides, 
IRefian(fº,A4º"!9))] < C.llfº]P + cIAPI A? to], (68) 


using equation (30), we have 


|Ref-n(14, Au + ASP) 
CAUSAS] + SUAS PI 
CotILFI UI + [IPI (69) 


Ref=n(1*, Ativ) 


IA IA 


and 
Ref-an(f!, AºiAw) | < C.||F||2 + elAPI| A? w||? for 0O<ó<1, (70) 


orfor0<q=< E of item (ii) Lemma 15 


Reft-an(f*, AºiAw))] < ColP la + (FP IalU Ia + PES 
Co IP IeslIU has + [1 Fã) (71) 


IA 


And besides, for0 <p <1 
|Rei(fº,20)] < C.)lF lã + elAJ2]|O]|P. (72) 


Therefore, from the estimates (65)-(69) together with the estimates (71) and (72) we 
conclude the proof of the item (ii) of this lemma. 
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Proof. (iii): Performing the duality product between equation (32) and AA?º30, taking 
advantage of the self-adjointness of the powers of the operator A and using (30), we obtain 


26—-3 1 


iJAPI|A 2 g!|? me —n(=iw = if2, 13-39) o emtadr ww, A 426-36) 
— AJ AS-to]2 + (14,1 42-39). (73) 


Asfor0<p=<1 wehave 30 -3<0< > 2q i< =— and 2p — 3 < E , taking imaginary 


part in (73) and using continuous immersions, for e > 0 exists Ce > 0 such that 
24—3 


— 29-3 
2 


BI2 < Csf]wl|? + || 43012 + | IO] + elAPIA"T 0] 
+ CIA Tu] + Cf. (74) 


INPILA 


Finally using estimative (33) and Lemma 7, we finish the proof of the item (iv) of this 
Lemma. 

Lemma 17 Let ô > 0. There exists Cô > O such that the solutions of the system (3)- 
(5) for lAI > 5, satisfy 


ã 1 
(0) Alo|Ê < GotFIalUla + Fla) for 5<ó<1, (75) 
(ii) |AIAZU|? < Cota O|? + [Alo]? + [PIU e) (76) 


1 
+C.[AIJA2Zul)? + ela 4?v]]? for 5<6<1 


Proof: (i): Multiplying by nAº'w equation (31) and by nw the equation (32) and 
adding the results, we have 
ioA|pw]|2 + a2(Au,ikv — f2) + a28(Av,w) — n(4º0,w) 
+ IJw]]2 + in(ÃO, A Iw) + nº (Aºv,iAA to — A? 2) 
+ arAST w|? + n(4º0,w) = a(f,w) + m(f4, A*w). (77) 
On the other hand. Carrying out the following duality products: first between equation 


(32) and CÊ qo-a u, taking advantage of the self-adjunction of the powers of operator A, we 
have 


ao) 28 28 28 
2B(w, Av) +22 (O, Al-év) + CL pato? + CL (Ag, Al-óy) = CE (pA Al-óy) 
am a am am 
(78) 
2/ 2 , 
CO (Ai-óy Ag) + il (Mm, AL-S9) + (Au, Al-$9) — Lj Atojp 
am n am a 
1 e 
+—(w, Al-?9) = —(f3, Al-?9) (79) 
am n 
a? a? a2, a 1 : 
— (Al-?9, Au) +i— (AO, Al ?u) + —(A2v, AZu) + a2(ikv — f2, Au) 
am am a 
= att! Aly). (80) 
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Adding the equations (77), (78), (79) and (80), we get 
ioAjw|? = aNif(A?u, AZv) — (Adu, AE AR) = IoÊ 
—a2B((Av, w) + (w, Av)) — in( (DuAtE 70, VA DI ER “w) 


VD 


+in T ol + mA 20, AS f2) — am 214º E “ul? 
2 28.4 
+m(f, At) — E (0, AI) — SEA? v||? (81) 


y 

+a(Pu) — TÉI(4O, Att) + (Aivto, AO) + TÊ(pS Aid) 
1 2 

cial, Al-$9) — cof Al-69) + (AL-S, Au)) + =IAtoI? 


qt vVIAlA2Zu) 


apta Pe w, AÉ-289) + ii Al-6g = 


-C (Aiy A? 2/4142 AÍ, MO sd qiod 
E (Ato, Abu) + 0!(Atfº, Abu) + Cfi,A u). 


From the identities: YZ E C, it is verified Imfz + z;= O and Imfi(z — z)+ = O, taking the 
imaginary partin (81), we get 


oAIw|? = Imfa? (Aiu, A2 2 f2) — in( 


a, DIA 


uy + AA o 


ZA 


+Imfny (42-20, A? e je cê qo, AAIto) + a(fê,w) 
+n($S, Au) + ET At-ty) — HA ai (82) 
a2 
—i-(Aw, Al) + =(f3, Al-ó9 il 420, VAA?u) 
a ) + S(, ) Am | 


4 
—T (Aju, Au) + a2(AZ f2, Adu) + Cf A-óW). 
a am 


For -< p<1,we have > — (Oo) qi from (54) of Lemma 15, we have 
28 1 lg 2 1 ; 
Im$ — O AA2?v) < CllFllelUla + ||Flpy for 3 <4<1(83) 


Besides, using (32) and considering <<: , from estimative (52) Lemma 15, 
estimative (60) Lemma 16 and estimative (49) ai 14, we have 


Inf (44%, 40) < CE(MAw,-iO — nAºv — anA?w + fº)|) 
< C(APIAw||? + JAPIJoIj? + Iwo]? + 1941/73 
+Cc|AIho]P + ejAI]faoi]? 
1 3 
< CotPIalUla + PIE) for 5<ó<7 (84) 
on the other hand, using (31), we have 
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Im (Ay, 49) <C|(-a? Au — 22840 +n4º0 — w + afº, AA-*6) | 


n 
. 1 A 
<C| — a2(|A|A2u, 
| vHI 


A : 
- a -4 3 E) 
+nA!|/6]| (viu, 574 9) + ad(fº, AA 0)|, 


3 


is 3 1 3 1 , 
Considering ;<y<Lwehave -—q<0,-p<0,--2p<and-p<y— applying Cauchy- 


A-90) — 2 B(XAT vw, AÉ-269) 


Schwarz and Young inequalities, continuous immersion and estimates (33), (34) of Lemma 
7, (50) of Lemma 14, (59) and (61) of Lemma 16, for e > O exists Ce > 0 such that 


Im [ (04, 40) |] < CstlFellUla + IPI) 


<ó<1. (85) 


mic 


+elAwl? for 


And, from -< p<1,we have 1 -q<o, then 


1 p = é 1 
[a lm(º, A=0)] < CAN fIP + 420") 
5 1 
< Co1lPlallUla + FI) for 5<ó<1 (8) 


Applying Cauchy-Schwarz and Young inequalities in (82) and using estimates (84)- 
(86), for e > 0 exists Ce > 00, such that 


1 Z 26-1 49 —3 1 7 - 
Av? < Cof Au A? PI + MIA Tv] + 4220] A? PP + FA] 
q =» g 26-1 ç 3. e ro 6 
+] + IF A! 0]]P + AT ||? + 1] 42-2º0]PP + [MJ A2-*0]P 
À is 1 1 1 1 
+ Azul]? + || A20]]? + || Azul? + 42 PIAU] + NA! ul) 


é 1,2 4-1 s g g 
+CEIMIA 7 0] + eJAKIA Tu] + ]o][2) + Cota UI + [1 PIA). 


(0) 1 0-1 13 11 1 
E I<5,q i< <0<5,5-20<5,5—- 0 <s0,/1 ss 


2.1 
For ;SP<1 wehave:2p-<-, E 


using continuous immersion, finish proof this item. 
Proof. (ii): Performing the duality product between a2BAv and w e using equation 
(30) and advantage of the self-adjointness of the powers of the operator A, we obtain 


2 B(Av,w) = a2B8(A?v,iXAZv — A? f2) = —ia2 BA AZv]2-a2B(A?v, A? f2). (87) 
On the other hand, multiplying by w equation (31) and by nAº'w the equation (32) 
and adding the results, we have 
io Ju]? + a? (Au, ido — f2) + a28(Av, w) — n(Aº0,w) 
+ Ihw||2? + in(ÃO, Aº law) + nº (APv,iXA Iv — AI f2) 
+ aÃ T w+n(Aº%0,uw) = (fw) + n(f!,Aºw). (88) 


Taking imaginary part in (88) and using (87), we have 
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aXIjwl|? = a28A/| 420]? + mA “ul + Im(a28(A2v, A? f2) 


A 
+ia2(n/ A[AZu, 
V É 


(om 


Av) + a2( Au, Aê f? d+ mº (4263 2u, A? f2) 


—0, V|AjA* tw) + a(fS,w) + n(f!, A w)). (89) 
Then 
a2BA|A?v E E v||2 = aAljw||? — Im[a25(A2v, A? f2) 


| A 
+ ia(v|AjA?u, — 
VI 


) 
— in(-—=0, 
K Di 
26-1 


Therefore, of (90), as for : <p<1lwelhave Eco, 2p- <>, p-1<0andgp-1 


Av) + (Au, AZ f?) + mn (42-20, A? f?) 


AJA) + (PE) + m(f!, Aw)). (90) 


< a, applying continuous immersion, Cauchy-Schwarz and Young inequalities, for e > O 
exists Ce > 0, finish proof this item. 

Theorem 18 The semigroup S(t) = e“ is analytic for q = 1. 

Proof: From Lemma(12), (47) is verified. Therefore, it remains to prove (48), for that 
it is enough to show that, let O > 0. There exists a constant Cô > 0 such that the solutions of 
the system (3)-(5) for IAI > 5, satisfy the inequality 


IU la 
Fla 


A| <OG ++ [JUl <Cs|FlalU] < oo. (91) 


Finally, considering & = 1 and using item (i) the Lemmas 7, 14 and 16, Lemma 17 
and item (iii) estimative (75) of Lemma 17, finish the proof of this theorem. 


4.2 Gevrey Class 

Before exposing our results, it is useful to recall the next definition and result pre- 
sented in [15] (adapted from [27], Theorem 4, p. 153]). 

Definition 19 Let tO > O be a real number. A strongly continuous semigroup S(t), 
defined on a Banach space H, is of Gevrey class s> 1 fort>t0, if S(t) is infinitely differentiable 
for t>t0, and for every compact set K C (t0,%0) and each u > 0, there exists a constant C = 
C(u, K) >0 such that 


ISP (Bllca <Cu"(n!)', foral te Kn=0,1,2.. (92) 


Theorem 20 ([27]) Let S(t) be a strongly continuous and bounded semigroup on a 
Hilbert space H. Suppose that the infinitesimal generator B of the semigroup S(t) satisfies 
the following estimate, for some O< W <1: 


Eai sup |A[“||(:AZ — B)!leqm < 00. (93) 
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Then S(t) is of Gevrey class s fort>0, for every s > = 

Our main result in this subsection is as follows: 

Theorem 21 Let S(t) = e? strongly continuos-semigroups of contractions on the 
Hilbert space H, the semigroups S(t) is of Gevrey class sii = 1, 2, for every si > CR 


(0) E 
for p E (5 5) and s2 > mas 5071 for P E É, 1), such that we have the resolvent estimative: 


(i) lim sup XP II(5A —-B)!Ilem <00, for : << a (94) 
(ii) lim sup [AE |I(ÕA — B) lc < 00, for : Ep, (95) 
se 
Proof: (i): Note that the estimate 
DEÉNGA = Bj" Pla = [NÉJUI < CollPlls for 5<6<I (99) 


implies the inequality (94). Therefore ffom now on we will show (96), for this purpose 
let us estimate the term IAI IIUII2H: Using estimative (75) Lemma 15 and (59) of Lemma 16 
in (76), we have 


DIABO < Col IelUla+ FI) for S<ó<1 (9 
Using estimates (34), (75), (59) and (97), for E > 0 exists Ce > 00, such that 
DU < CAPA for S<ó<i 
Hence, in particular for p E (5 5). we arrive at 
IO <ColFla for S<6<5 (98) 


Proof: (ii): We assume À E R with IA > 1, we shall borrow some ideas from [21]. Let 
us decompose w as w = w1 + w2, such that 


iodw, + Au, = afº in D(Aº) 
and 
iakw = —a2 Au — a284v +n4º0 —- w + Aw, in D(Aº). (99) 


Firstly, applying the product duality the first equation in (99)1 by w1, then by Awi, 


we have 


Alo] + A 42 w1]] + | Au4]] < CIP ha. (100) 


1 us 1 = ai pm " 
Asfor;<p<1 wehave: 0 < =— s ;and from AZEiwa = AZ w- AZ w1, using 


estimative (100) and estimative (51) Lemma 15, we have 
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CotIA Two] + A am]? 
CotILF Ia |hex + INF IES, 


E ua) 


IA IA 


asforisq<1,we have —1 < *1, then 


IA ua? < 
< o (NE PIU + FI) for 5 << (101) 
Now applying the operator A - : in the second equation of (99) and using Lemmas (7) 
and (15), estimates (33) and (100) and as —1 < — we have 
PIA dual < O(NAtujP + AtoJP + NAC hONÊ + A tu]PA + Abe 
< NTE UPIalUla + FI) for 5<6<1 
Therefore 
A-Iwo||2 < Co for : <ó<1. (102) 
By interpolations inequality Theorem 4, since — i< 0=< a, using estimates (101) 


and (102), we derive 


wa] <C(IA2wa|2) E (| AT wo]2)35 
< eq a te [2 PI] Ia + | FI] 
1 
< GN ÉS PaUla + PIE] for 5<6<I (103) 


On the other hand, as Ilwll2 < C(llw1II2 + Ilw2II2) and as for is P=<1we have -2 < 


29-1 using estimates (100) and (103), we get 
1 
E NE for =<6<1 
ol? < os a] for 5<ó<1, 
equivalently 
2- 1 
pio? < CNS [NS [Ulla + Fla) for 5<ó<1 (104) 


Now applying the estimate (104) in estimative (76) of Lemma 17 for is p<s1,weget 


DIIAZo|ê < Cs (INIJJo]? + JA JU la + [171 + 11 7]la[U ha) 
1 
+ CoMIA ul? + clAtv|? for 5<6<1 (105) 


for other hand, using (34) and (59), Lemmas 7, 16 respectively, we get 


2- 


DIIIAZo]ê < JAP? 
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On the other hand, from estimative (35) of Lemma 7 and 


2 < PAP Au]? < Cota + FI j 
< Cota 


0<ó<1, 


itis equivalent to 


<6<1, (107) 


AO | tem 


DIIAZu]? < 


then from (106), (107) and 
a|AI|AZu + aA2v]? < a2A||AZul2 + aZo|A||AZo]f2, 
we get 


af AZutaA?v I|2 < DR 


Besides, from (59) of Lemma 16 


DRE 


< co 
< Co qalÉ 


[Ula+IlFlã) for 0<9<1. 


Equivalent to 


1 
JUla+|Flã) for 5<6<1 (109) 


Ijoj? < (x 3 


From 
Alto + avolf? < Alo? + ajMIolfP, 


Finally, adding the estimates (106), (107), (108) and considering estimative (104), 


we get 


JUIZ < s(NÉ <ó<1 (110) 


Therefore, finish the proof of inequality (96). 

We finish the proof of this theorem. 
4.3 Lack Gevrey Class and Analiticity 

Theorem 22 The semigroup S(t) = e? does not support Gevrey classes and is not 
analytic when q E o, 3 

Proof: Since the operator A is strictly positive, selfadjoint and it has compact 
resolvent, its spectrum is constituted by positive eigenvalues (on) such that on — o as 
n-> e, Forn € N we denote with en an unitary D(Aº)-norm eigenvector associated to the 


eigenvalue on, that is, 


Ciências exatas e da terra: teorias e princípios 2 Capítulo 2 


50 


Aen = Ontn, Ilenlloçso =lenl=1, nen. (111) 
sa the eigenvalues and eigenvectors of the operator A as in (111). Let Fn 


=(0,0, 4,0) EH. The solution U = (u, u, w,, 8,) of the system (iA,I-B)U = F, satisfies 


1 a) 


va = Al, Un = — Aun, Aºen = oYen fory E R and the following equations 
“iai un + a2 Au, + ia? BA Au, — Aun — n4º0n = —en, 
inAn Aun — anA2 Aun + iAnOn + A0n = 0. 
Equivalently, 
(i[aX2I — 28h A] +AI — 02 Aju, +nAº0, = en 


LinAn A? — anAZ AºJu, + fil + AJ = 0. 
Let us see whether this system admits solutions of the form 
Un = UnCn; 0, = VnCn; 


for some complex numbers y | and v,. Then, the numbers y , v, should satisfy the 


algebraic system 


fifa — a28An0n] + A — onsun+ not = 1, (112) 
tinAnos - anão 04 )Un E tin St On == 0. (113) 


At this point, we introduce the numbers 


je [É 2 
ia [55] An E 
Then 
= 2 2 B—-a na? |. 
[on] Ea |An| , A — q? On = | B pão nos = ES Ar 
1+ó 
; 2 ja 142 O N| 22 
inAnos, — anA£o 04 = =1 F Es An +26 — Es AS ro 


and in +ton=iAh+ E | AZ. (115) 
Using identities of (114) and (115) in system (112)-(113), we obtain 


B-a,2 no? 26 
Em hm + ( 2689 An Vn 


1+ó 
. na? 1426 A 2424 E A 2 = 
[i ao pan 2509 ed hum + [im + PET hm = (. (117) 


R (116) 


Then 
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2 24 
A= N ceepo 
no" aro CM ao ; DV 
255 a2é 8% Am P ida + FETAL 
a(8— a) | 4 tin? o as PQ Po? agr 
Et + [apo pro SO = gente] nm 
Therefore: 
For pelo, +] andB-a>0, we have 
i 1 
JA] = [Anjmoxt48) = A]! for 0<6< 5 (119) 
Besides: 
into : 
| An] = | ab | = [in +=] =|) for 0<6<>. (120) 
O D+ A di ” 
a2 8 n 
Then 
Au o 1 
tn| = Ge x |An|? for 0<94< 3º (121) 
From the definition of IIUIIH and estimative (121), for K > 0 we have 
1 
[Una > K|A2v,] = |onliMun) = Kin)? for 0<64< 5 (122) 
and 
sia: oil 1 
|An||Unlla > KA | for 0<G< 3º (123) 


Remark 23 |) If there is a Gevrey class for S(t) when 4 E To, =), there must exist a 
E (0, 1) such that the identity 


lim sup ||" (1X — B)Hlgem) < 00. (124) 
An |=>00 
is verified. But, multiplying both sides of the inequality (122) by IAnIY for 0 <W <1, 
we have 
W d W l E 
In |U, ll > K|M| for 0 < Q < 2º (125) 
then, IAn!VIIUniIH—= o approaches infinity as IÂnl — o», Therefore for q in [0, 5] the 
(124) condition fails. Consequently the semigroup S(t) does not admit a Gevrey class for & 
elo, =). 
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Il) From inequality (123), HUnilH — co approaches infinity as lAÂnl = o, There- fore 
for p in [0, a) S(t) does not satisfy the condition (48) of Theorem 13, so S(t) is not analytic 
when 4 E TO, =]. 
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RESUMO: O conhecimento dos níveis de irradiância solar que atingem a superfície terrestre 
é um pré-requisito para várias aplicações. Uma das características da radiação solar é o 
alto grau de variabilidade temporal, que está relacionado aos padrões climáticos que atuam 
em cada local. Sendo assim, o presente trabalho propôs realizar um estudo de caso para 
a cidade de Natal (5,83ºS; 32,20ºW) sobre o comportamento da radiação de onda curta, 
utilizando a Irradiância Global Horizontal (GHI) e a Irradiância Difusa Horizontal (DHI), por 
meio de avaliações estatísticas para dados mensais de satélite disponibilizados pelo site 
SoDa no período de 2015 a 2020. Os resultados obtidos através dos parâmetros utilizados 
foram as componentes do boxplot, ou seja, máximos, mínimos, médias e medianas. 
Também foi empregada uma regressão linear, a fim de identificar tendências positivas ou 
negativas. Os valores obtidos tanto para o GHI quanto para o DHI mostraram que ambas as 
séries estudadas apresentaram uma sazonalidade bem definida. O GHI apresentou pouca 
dispersão de dados com seu valor máximo em torno de 220560 Wh/m2 e seu mínimo em 
135481 Wh/m2. Quando comparados os valores de verão e inverno notou-se que há uma 
diferença de até 21,7% nos valores mínimos. Já para DHI os dados apresentaram maior 
dispersão, principalmente entre os meses de janeiro e março, o valor máximo ficou em 13918 
e o mínimo em 46944 Wh/m2. Quando comparadas as duas estações, verificou-se que a 
maior diferença fica nos valores médios com 19,1%. No que diz respeito à regressão linear, 
verifica-se a presença de uma tendência negativa para o GHI e a completa ausência de 
tendência no que se refere ao DHI. 

PALAVRAS-CHAVE: GHI, DHI, parâmetros. 


Ciências exatas e da terra: teorias e princípios 2 Capítulo 3 


56 


STATISTICAL ANALYSIS OF THE BEHAVIOR OF SHORT WAVE RADIATION 
FOR THE CITY OF NATAL, RN 


ABSTRACT: Knowledge of the levels of solar irradiance that reach the earth's surface 
is a prerequisite for various applications. One of the characteristics of solar radiation is 
the high degree of temporal variability, which is related to the climatic patterns that act in 
each location. Thus, the present study proposed to conduct a case study for the city of 
Natal (5.83ºS; 32.20ºW) on the behavior of short wave radiation, using Horizontal Global 
Irradiance (GHI) and Horizontal Diffuse Irradiance (DHI), through statistical evaluations for 
monthly satellite data made available by the SoDa website in the period 2015 to 2020. 
The results obtained through the parameters used were the components of the boxplot, 
i.e., maximums, minimums, averages and medians. A linear regression was also used in 
order to identify trends. The values obtained for both GHI and DHI showed that both series 
studied presented a well-defined seasonality. The GHI showed little data dispersion with 
its maximum value around 220560 Wh/m2 and its minimum in 135481 Wh/m2. When the 
summer and winter values were compared, it was observed that there is a difference of up to 
21.7% in the minimum values. For DHI, the data showed greater dispersion, mainly between 
January and March, the maximum value was 13918 and the minimum at 46944 Wh/ma2. 
When comparing the two seasons, it was found that the largest difference is in the mean 
values with 19.1%. With regard to linear regression, there is the presence of a negative 
trend for GHI and the complete absence of trend with regard to DHI. 

KEYWORDS: GHI, DHI, parameters. 


11 INTRODUÇÃO 


Um dos desafios da sociedade contemporânea, segundo Holf (2011), é suprir 
as demandas de água e energia, decorrentes de diversos fatores, entre eles está o 
desenvolvimento econômico, o crescimento da população, processos de degradação 
ambiental e as mudanças climáticas em curso. De acordo com o demonstrado por Chen 
e Chen (2016), observa-se que as regiões urbanas utilizam 70% da energia global e 
necessitam de uma grande quantidade de água para o consumo humano e geração 
de energia. Segundo De Oliveira e Trindade (2018), nos moldes atuais, a crescente 
demanda e a exploração intensa de combustíveis fósseis tornam a produção de energia 
insustentável. 

Diante desse cenário, no qual há uma grande demanda energética somada 
à insustentabilidade da produção de energias não renováveis, a utilização de fontes 
renováveis, tais como eólica, fotovoltaica e termelétrica, tem sido cada vez mais recorrente. 
Sendo assim, em consequência a esse fato, estudos que visam identificar e apresentar 
locais com grandes potenciais na geração desses tipos de energia, têm tornado-se cada 
vez mais comuns, de modo que a quantidade tem crescido de forma excepcional. Dessa 
forma, entender o comportamento das variáveis que influenciam e geram essas fontes 


energéticas é de suma importância. Como exemplo a isso, os estudos sobre os índices de 
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irradiância solar são importantes na geração de energia solar fotovoltaica, a qual é gerada 
em condições de altos níveis dessa variável. 

Assim, o conhecimento sobre os níveis de irradiância solar que atingem a superfície 
terrestre é um pré-requisito para várias aplicações da energia solar na indústria, na agricultura 
e em outras atividades relevantes sob o ponto de vista socioeconômico. Para se ter esse 
conhecimento é necessário mapear o comportamento dos dados disponíveis na superfície 
terrestre. Varela (2007) discorre que esse procedimento de análise em todo o território 
nacional é inviável, pois exige incentivos tecnológicos que nem sempre são atendidos, 
além da mão de obra especializada para instalação e manutenção dos instrumentos de 
medição de irradiância solar. Segundo Pereira (2006), o Brasil, por situar-se na região 
tropical próximo à linha do Equador, recebe grande incidência de radiação solar durante 
todo o ano ao longo de seu território, além de uma menor variabilidade interanual. 

Uma das características da radiação solar é um alto grau de variabilidade temporal, 
que está relacionada aos padrões climáticos que atuam em cada local. Sendo assim, 
o objetivo deste trabalho é avaliar o comportamento da radiação de onda curta (ROC), 
especificamente a Irradiação Global Horizontal (GHI) e a Irradiação Difusa Horizontal 
(DHI), as quais cnegam à superfície da cidade de Natal entre janeiro de 2015 e dezembro 
de 2020, fazendo uma comparação dos resultados obtidos para o inverno e verão. 


21 MATERIAL E MÉTODOS 


2.1 Área de Estudo 


A área de estudo refere-se à capital do Rio Grande do Norte, Natal-RN (5,83º8, 
32,20º W)(Fig. 1). Pertence à mesorregião do Leste Potiguar (IBGE, 2020), possuindo uma 
área territorial de 167,401 km2 com uma altitude média de 46 m em relação ao nível médio 
do mar. Além de uma costa litorânea, a cidade possui pontos de vegetação de Mata Atlântica. 
Recentemente, Alvares (2014) elaborou, utilizando os critérios de Kôppen, a classificação 
climática do Brasil, de modo que Natal é dito como um local com o clima do tipo As: tropical, 
com o verão seco. Segundo o INMET, 2021 apresentou as médias de precipitação anual e 
de temperatura compensada iguais a 1721,4 mm e 26,4 ºC, respectivamente. 
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———— Limite dos Municípios 


Figura 1: Mapa referente ao município de Natal, RN. Fonte: AMORIM, 2020. 


2.2 Dados Utilizados 


Foram coletados dados das variáveis lrradiância Global Horizontal (GHI, Global 
Horizontal Irradiance) e Irradiância Difusa Horizontal (DHI, Diffuse Horizontal Irradiance), 
propondo uma análise para o intervalo entre 2015 e 2020, seguindo o calendário juliano. 
As coordenadas de latitude e longitude citadas no tópico anterior serviram como ponto de 
coleta inserido no sistema on-line do Solar Energy Services for Professionals, o qual é um 
serviço disponibilizado pelo site Solar Radiation Data (SoDa). Foram coletadas as variáveis 
na unidade Wh/m2. O SoDa fornece séries temporais de irradiações globais, diretas e 
difusas na superfície horizontal e irradiação direta no plano normal, para as condições 
climáticas reais, bem como para condições de céu claro. A plataforma disponibiliza como 
cobertura geográfica a Europa, África, Oceano Atlântico e Oriente Médio. 


2.3 Análises Estatísticas 


Na análise do banco de dados da pesquisa foram utilizados métodos estatísticos 
destinados a avaliar o grau de variação dos dados próximos à média. Contudo, utilizam-se 
medidas de dispersão, tais como: desvio médio, variância e desvio padrão. Além disso, os 
conjuntos de dados de média e mediana também são incluídos para representar o valor 
central do conjunto de dados total. 
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2.3.1 Média 


A média é dada pela soma de todos os elementos que estão sendo considerados, 
dividido pelo número de elementos (n). É dada pela Eq. (1): 


1 n “ 
X = Aid Eq. 1 


2.3.2 Variância e Desvio Padrão 


Variância e desvio padrão são medidas de dispersão de uma amostra, e segundo 
Gurgel (1958), servem para dar uma ideia o quão distante a amostra está da média. O 
desvio padrão amostral é dado pela Eq.(2) e a variância amostral é o resultado obtido pelo 
desvio padrão ao quadrado. 


s=|-|3 (K-*? Eq. 2 


n-1 
2.3.3 Mediana 


A mediana é uma medida de posição que separa o conjunto de dados ordenados 
exatamente na metade, deixando 50% dos valores abaixo e 50% acima. 


2.3.4 Quartis 


Em um trabalho recente, Morettin (2017) propôs que os quartis são medidas de 
posição que separam o conjunto de dados ordenados em 4 partes iguais. O primeiro 
quartil (Q1) indica o ponto de separação onde 25% da amostra está abaixo desse valor. O 
segundo quartil (Q2) é a mediana, com 50% acima e abaixo. Já o terceiro quartil (Q3), é o 
ponto onde 75% da amostra está abaixo de seu valor. 


2.3.5  Boxplot 


De acordo com Morettin (2017), o gráfico Boxplot ajuda a visualizar como se portam 
os dados em relação à simetria e dispersão. Na sua construção é considerado o eixo 
vertical (y, ordenadas) do gráfico o valor dos dados. Primeiro, é traçado um retângulo onde 
a parte superior é o Q3 e a inferior o Q1. Um traço indica a mediana (Q2). Em seguida, é 
criado uma linha do Q3 até o limite superior (valor máximo amostrado) e do Q1 até o limite 
inferior (valor mínimo). Os pontos discrepantes, chamados outliers, são valores que ficam 
abaixo do valor Q1-1.5*(Q3-Q1) ou acima de Q3+1.5*(Q3-Q1). 
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Figura 2 - Informações contidas no Boxplot. Fonte: Oper*, 2021. 


2.3.6 Regressão linear simples 


A análise de regressão é considerada simples quando envolve duas variáveis: a 
variável dependente/resposta e a variável independente/explicativa. Para sua aplicação 
é necessário que os dados sejam oriundos de uma amostra aleatória. Como parte de 
avaliação da técnica se utiliza o gráfico de dispersão, o qual é utilizado para mostrar a 
relação entre as duas variáveis, podendo ela ser linear ou não (DRAPER, SMITH, 1998; 
KUTNER et al., 2004 ). Segundo Azevedo (2012), quando se tem apenas uma variável 
independente, o modelo de regressão linear é da forma: 

Yi=a+ BhXi Eq. 3 
Sendo: 
Y, - variável resposta (dependente); 
a - constante que representa a interceptação da reta com o eixo vertical; 
X, - variável independente; 
P - parâmetro 


31 RESULTADOS 


Nesta seção serão apresentados os resultados referentes às análises das séries 
de Irradiância Global Horizontal e Irradiância Difusa Horizontal (Wh/m2) para a cidade de 
Natal. Com auxílio da Fig. 3, são apresentadas as características do GHI para o município 
de Natal. Dessa forma, foi possível identificar a existência de uma tendência linear para a 
série, aplicando uma suavização de sua variabilidade, que ilustra o comportamento médio 
ao longo do tempo. 
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Através de uma análise visual da Fig. 3a observa-se a existência de uma sazonalidade 


muito bem marcada, ou seja, os dados mostram que há períodos com maiores incidências 


e outros com menores. Durante a observação notou-se que no mês de outubro foram 


registradas as maiores incidências, enquanto que no mês de junho as menores, sendo 


correspondentes a 220560 e 135481 Wh/m2, respectivamente. Para melhor entender 


o comportamento do GHI, foi gerada uma linha azul (Fig. 3a, Atlas) para acompanhar a 


transição entre os meses. Quando a análise é voltada para o período chuvoso que ocorre 


entre os meses de março e maio, nota-se um período de transição dos maiores para 


os menores valores. Entretanto, para a transição inversa, ou seja, dos mínimos para os 


máximos, ocorre uma mudança abrupta no período seco, o qual decorre entre os meses de 


setembro e novembro. 
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Figura 3 - Comportamento da Irradiação Global Horizontal (GHI): a) Gráfico com destaque ao boxplot 
referente ao comportamento do GHI dos 12 meses para um período de 5 anos. b) Gráfico referente 
com a série de 5 anos com destaque a regressão linear. 


ESTATÍSTICAS RESULTADOS (Wh/m2) 
Mínimo 135481 
1º Quartil 172611 
Mediana 192207 
Média 187785 
3º Quartil 204513 
Máximo 220560 


Tabela 1 - Resultados do boxplot referente ao GHI. 


GHI DIFERENÇA (%) 
Mínima 21,7 
Média 15,8 
Máxima 9,8 


Tabela 2 - Valores das diferenças entre o período do verão e inverno para o GHI 
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Na Fig. 3b foi aplicado um cálculo de regressão linear simples, a fim de verificar 
se os dados apresentavam algum tipo de tendência. Dessa forma, foi possível identificar 
a ocorrência de uma mudança mínima negativa em relação aos valores, e o ponto de 
alteração na tendência foi identificado no mês de julho no ano de 2017. Quando realizada 
a comparação entre o inverno e o verão austral, compreendidos entre os meses de junho 
e agosto e entre dezembro e fevereiro, respectivamente, é observado uma diferença com 
uma taxa de até 21,7% entre eles, como é possível observar para os mínimos valores 
(Linha 1, Tab. 2). A média para a análise do período completo ficou em torno de 187785 
Wh/m2 (Tab. 1). 

Os dados de Irradiação Difusa Horizontal (Fig. 4a) apresentaram o seguinte 
comportamento, a sazonalidade assim como o da Fig. 3 manteve um comportamento bem 
definido em relação ao seu vale (menor valor da série) que ocorreu no mês de junho com o 
valor de 46494 Wh/m2, mas quando observamos o pico (maior valor da série), observou-se 
a ocorrência em maio com 73421 Wh/m2, mas com um adendo, os maiores valores foram 
observados no mês de janeiro. Quanto ao box plot é possível observar que eles apresentam 
uma maior dispersão de valores, principalmente entre janeiro e março. No período seco e 
chuvoso, onde acontecem as transições de máxima para mínima e de mínima para máxima 
nas observações, os boxplots apresentam as mínimas dispersões da série. 

Sobre os valores encontrados nos métodos estatísticos, verificamos que a média 
ficou em torno de 59000 Wh/m2 e que quando foi aplicado à análise de regressão linear não 
foi possível identificar nenhuma mudança na sua tendência seja ela positiva ou negativa 
(Fig. 4b), apesar de não apresentar tendência, a série mostra valores de diferença acima 
dos 14% e com o seu máximo chegando a 17,4% no valores mínimos quando comparamos 


os resultados apresentados entre o inverno e verão. 


75000 
75000 


— -»- Atas 


70000 
70000 


65000 
1 
65000 


60000 


Irradiação Global Difusa - ( Wrym?) 
60000 
1 


55000 
1 
55000 


50000 
50000 


J FM AMI J AS O ND 2015 2016 2017 2018 2019 2020 2021 


Figura 4 - Comportamento da Irradiação Difusa Horizontal (DHI): a) Gráfico com destaque ao boxplot 
referente ao comportamento do DHI dos 12 meses para um período de 5 anos. b) Gráfico referente 
com a série de 5 anos com destaque a regressão linear. 
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ESTATÍSTICAS RESULTADOS (Wh/m2) 
Mínimo 46944 
1º Quartil 54113 
Mediana 59449 
Média 59284 
3º Quartil 62921 
Máximo 73918 


Tabela 3 - Resultados do boxplot referente ao DHI 


GHI DIFERENÇA (%) 
Mínima 17,4 
Média 19,1 
Máxima 14,4 


Tabela 4 - Valores das diferenças entre o período do verão e inverno para o DHI 


41 CONCLUSÃO 


Neste trabalho foi apresentada a caracterização e o comportamento para os dados 
de Irradiância Global Horizontal e Irradiância Difusa Horizontal para a cidade de Natal. 
Os resultados apresentados demonstram valores significativos quanto à comparação do 
período do inverno e verão, com diferenças de até 21,7 % para o GHl e 19,1% para o DHI. 
Outro ponto destacável é a presença de uma tendência negativa para o GHI e a completa 
ausência de tendência no que se refere ao DHI. Porém, para determinar o que ocasionou 
isso, é necessário pesquisas adicionais para determinar a causa desses comportamentos. 
Esses estudos são complexos devido ao grande número de combinações de fatores 
envolvidos. Propõe-se que, em trabalhos futuros, sejam aplicados métodos mais robustos 
para testes de tendências, além de pesquisas de cunho meteorológico para complementar 
os resultados que venham ser obtidos. 
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RESUMO: Seguindo o Acordo de Paris 
na COP21, 197 países têm como objetivo 
reduzir significativamente as emissões de 
gases de efeito estufa em todo o mundo. 
Como as cidades são grandes emissoras 
de gases de efeito estufa, 25cidades 
americanas estão lançando o American 
Cities Climate Challenge. O objetivo é 
reduzir as emissões de gases de efeito 
estufa em 32% até 2025. A luta contra a 
mudança climática é conduzida em várias 
frentes, tanto para mitigar os efeitos do 
aquecimento global quanto para adaptar 
nossas comunidades a um estilo de vida 
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em mudança. As medidas escolhidas 
dependem, naturalmente, dos problemas 
específicos das diferentes regiões do mundo 
e o impacto previsto certamente levará ao 
desenvolvimento sustentável. São medidas 
estratégicas para aumentar a resiliência e a 
inteligência das cidades, ou seja, monitorar 
e controlar em tempo real a infra-estrutura e 
os serviços gerados pelas cidades, levando 
a uma redução no consumo de energia e na 
poluição. A instalação de uma extensa rede 
de sensores e dispositivos para monitorar 
a qualidade do ar. Estabelecimento de 
ciclovias para incentivar o uso da bicicleta. 
Dar prioridade aos veículos elétricos sobre 
os veículos a gasolina. Tornar o transporte 
público acessível enquanto reduz as 
viagens de pessoas. Este artigo é parte de 
um paradigma construtivo de adaptação 
e mitigação da mudança climática, com 
o objetivo de apresentar estratégias para 
combater a mudança climática no âmbito 
do curso Inovação e Novas Dinâmicas 
Sociais e Econômicas nas Cidades: 
Cidades Inteligentes e Cenários de 
Inovação Urbana, do Programa de Pós- 
Graduação em Desenvolvimento Regional 
e Meio Ambiente, no âmbito do programa 
bolsas Brasil — PAEC — OEA — GCUB, da 
Universidade do Estado de Minas Gerais. 
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PALAVRAS-CHAVE: Mudança climática; Mitigação; Adaptação; Cidades inteligentes; 
Desenvolvimento sustentável. 


MITIGATION AND ADAPTATION: TWO STRATEGIES TO COMBAT CLIMATE 
CHANGE 

ABSTRACT: Following the Paris Agreement at COP21, 197 countries are aiming to 
significantly reduce greenhouse gas emissions worldwide. Because cities are large emitters of 
greenhouse gases, 25 American cities are launching the American Cities Climate Challenge. 
The goal is to reduce greenhouse gas emissions by 32% by 2025. The fight against climate 
change is conducted on several fronts, both to mitigate the effects of global warming and to 
adapt our communities to a changing lifestyle. The measures chosen depend, of course, on 
the specific problems in different regions of the world and the expected impact will certainly 
lead to sustainable development. These are strategic measures to increase the resilience 
and intelligence of cities, i.e. to monitor and control in real time the infrastructure and 
services generated by cities, leading to a reduction in energy consumption and pollution. The 
installation of an extensive network of sensors and devices to monitor air quality. Establishing 
bicycle lanes to encourage bicycle use. Giving priority to electric vehicles over gasoline- 
powered vehicles. Making public transportation affordable while reducing human travel. This 
article is part of a constructive paradigm of climate change adaptation and mitigation, with the 
objective of presenting strategies to combat climate change in the course Innovation and New 
Social and Economic Dynamics in Cities: Smart Cities and Urban Innovation Scenarios, of the 
Postgraduate Program in Regional Development and Environment, under the Brazil - PAEC - 
OAS - GCUB grant program, of the Minas Gerais State University. 

KEYWORDS: Climate change; Mitigation; Adaptation; Smart cities; Sustainable development. 


11 INTRODUÇÃO 


Durante milhares de anos, o clima da Terra tem variado com o tempo e o lugar. 
As mudanças observadas são geralmente espalhadas por longos períodos de tempo, o 
que torna difícil para os seres humanos percebê-las em um determinado momento. Nas 
últimas décadas, no entanto, a mudança climática parece ter se acelerado. De fato, a 
mudança climática é uma realidade inegável. Devemos agir urgentemente para evitar uma 
acumulação irreversível de gases de efeito estufa (GEE) na atmosfera e o aquecimento 
global que poderia ter consequências colossais para a economia e a sociedade mundial. 
Para este fim, 197 países assinaram e adotaram o Acordo Paris no final da Conference Of 
the Parties 21 (COP 21) em Paris. O acordo visa a reduzir significativamente as emissões 
globais de gases de efeito estufa e limitar o aumento da temperatura global a 2 ºC neste 
século, enquanto busca formas de reduzir este aumento para 1,5 “C (Nações Unidas, 2015). 
O acordo estabelece compromissos para que todos os países reduzam suas emissões e 
cooperem para se adaptar aos efeitos da mudança climática, ao mesmo tempo em que 
apela para que fortaleçam seus compromissos ao longo do tempo. Ela proporciona uma 
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oportunidade para os países desenvolvidos ajudarem os países em desenvolvimento em 
seus esforços de mitigação e adaptação climática, ao mesmo tempo em que estabelece 
uma estrutura para o monitoramento transparente e a elaboração de relatórios sobre as 
metas climáticas que foram estabelecidas. Já 25 cidades americanas estão no caminho 
certo para reduzir coletivamente suas emissões de gases de efeito estufa em 32% até 
2025, superando as metas dos acordos internacionais através do programa American 
Cities Climate Challenge. 

É neste contexto que o desenvolvimento das atividades humanas modernas, 
incluindo transporte, indústria, desmatamento e agricultura, é responsável pela emissão 
massiva de três principais de gases de efeito estufa (GEE), dióxido de carbono (CO,), 
metano (CH,) e óxido nitroso (N,0). O objetivo geral deste artigo é apresentar estratégias 
para combater a mudança climática e as condições para implementar essas estratégias 
para reduzir as emissões de gases de efeito estufa. Por este fato, resume o trabalho 
recente sobre os vínculos entre clima e clima. Primeiro, analisa os desafios da mudança 
climática e nos concentra nos mecanismos implementados pelas 25 cidades americanas 
que estão no caminho certo para superar as metas climáticas de Paris até 2025, de acordo 
com Udasin (2021). Finalmente, mostra como a adaptação e a mitigação efetiva são dois 
estratégias similares que podem combater a mudança climática. 


21 REFERENCIAL TEÓRICO 


2.1 Gases de Efeito Estufa (GEE) 


Desde que nós humanos aparecemos no planeta pela primeira vez, temos 
mudado nosso meio ambiente e explorado os recursos naturais à nossa disposição para 
melhorar nossas condições de vida. Atualmente, muitas atividades humanas, destinadas 
principalmente a garantir um padrão de vida confortável, estão emitindo muitos gases 
de efeito estufa (GEE) na atmosfera. Estes gases de efeito estufa estão naturalmente 
presentes na atmosfera. Eles formam uma camada ao redor da Terra, permitindo-lhe reter 
seu calor: este é o efeito estufa. O sol aquece a Terra, que então irradia parte de seu 
calor de volta para o espaço. Os GEE na atmosfera retêm parte desse calor, impedindo 
que ele escape de volta para o espaço. Isto mantém a temperatura média em nosso 
planeta em 15 ºC. Sem ela, a temperatura seria de cerca de -18ºC, o que não suportaria 
a vida como a conhecemos (Chapron e al., 2010). Desde o início da revolução industrial, 
por volta de 1750, o efeito estufa foi amplificado pela liberação de grandes quantidades 
de GEE na atmosfera. O uso massivo de combustíveis fósseis como petróleo, carvão 
e gás natural, o desmatamento, certos processos industriais, a deposição de resíduos 
em aterros, o uso de CFC em sistemas de ar condicionado, a agricultura intensiva, os 
resíduos animais, a inundação de vales, os incêndios florestais, etc., desempenharam um 
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papel importante no aumento das emissões de gases de efeito estufa desde o início da 
era industrial. Como consequência do aquecimento global, a temperatura na superfície 
do globo está aumentando muito rapidamente, o que pode levar a mudanças climáticas 
significativas com repercussões sobre os seres humanos e seu meio ambiente. A fim de 
combater o fenômeno global da mudança climática, todos os níveis de governo precisam 
ser mobilizados para reduzir e mitigar ao máximo as emissões de GEE e proteger nosso 
planeta tanto quanto possível. 


2.2 Mitigação e adaptação 

Como as mudanças climáticas e sua atribuição sobre as atividades humanas 
são consideradas uma realidade por quase todos os cientistas que estudaram estas 
questões, a mobilização internacional sobre este assunto se acelerou. Mitigação e 
adaptação são duas estratégias para lidar com a mudança climática. A mitigação propõe 
ações para engajar os emissores de GEE no interesse coletivo do mundo, reduzindo 
qualquer risco de mudança climática. O conceito de adaptação é definido pelo Terceiro 
Relatório de Avaliação do GIEC como o ajuste em sistemas naturais ou humanos em 
resposta a estímulos climáticos atuais ou futuros ou seus efeitos, a fim de mitigar efeitos 
adversos ou explorar oportunidades benéficas. A adaptação à mudança climática é uma 
abordagem proativa para reduzir a vulnerabilidade aos impactos diretos e indiretos da 
mudança climática (Glossário GIEC): uma ampla gama de medidas de adaptação pode 
ser implementada para lidar com as consequências da mudança climática. Entretanto, 
para ser eficaz, a adaptação deve ser feita simultaneamente em vários níveis. É 
fundamentalmente local, pois os impactos diretos da mudança climática são e serão 
sentidos localmente. De acordo com uma fórmula comum, a mitigação visa evitar o 
incontrolável e a adaptação para administrar o inevitável. Como Holdren (2010), ex- 
presidente da Associação Americana para o Progresso da Ciência, graficamente coloca, 
temos três opções: mitigar, adaptar ou sofrer. Faremos um pouco de todos os três. A 
questão é qual combinação escolheremos: quanto mais mitigarmos, menos teremos que 
investir na adaptação, e menos sofrimento haverá (Holdren, 2010). 

Na implementação de políticas de combate à mudança climática, é necessária 
uma certa vigilância, pois a mitigação e a adaptação podem, em alguns casos, revelar- 
se contraditórias. Um exemplo frequentemente citado de uma medida de adaptação que 
contradiz os objetivos de mitigação é o uso de ar condicionado durante o tempo quente, 
com suas emissões de GEE associadas. Da mesma forma, uma política de densificação 
urbana para reduzir as emissões relacionadas ao transporte que não leva em conta o risco 
de aumentar o efeito ilha de calor urbana é um exemplo de conflito entre adaptação e 
mitigação. É por isso que é essencial implementar as duas estratégias de forma integrada, a 
fim de limitar os conflitos e buscar sinergias. Existem muitas dessas sinergias: por exemplo, 
a ecologização urbana limita o efeito ilha de calor urbana, melhora a drenagem urbana e 
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promove a biodiversidade (adaptação), ao mesmo tempo em que reduz a necessidade de 


resfriamento e armazenamento de carbono (mitigação). 


2.3 Entendendo a adaptação e mitigação como uma estratégia de 
gerenciamento de risco climático 

Basicamente, os dois principais cursos de ação, mitigação e adaptação, são 
estratégias para gerenciar os riscos associados à mudança climática. Ambos visam a 
reduzir os gases de efeito estufa, ou seja, a probabilidade de impactos adversos induzidos 
pelo clima nas sociedades humanas e nos ecossistemas. Para reduzir os riscos climáticos, 
as políticas de mitigação tentam limitar os distúrbios humanos ao sistema climático (por 
exemplo, reduzindo as emissões de GEE), ou mesmo alterar deliberadamente o sistema 
climático para contrariar a mudança climática antropogênica ou compensar alguns de seus 
efeitos (por exemplo, através da possível implementação de técnicas de geo-engenharia 
emergentes e altamente controversas). As estratégias de adaptação, por outro lado, 
visam a agir sobre os sistemas humanos, em particular, tentando reduzir sua exposição e 
vulnerabilidade. Alguns autores preferem falar sobre o aumento da resiliência, ou seja, a 
capacidade de um sistema de absorver distúrbios enquanto mantém sua estrutura básica e 
seus modos de operação (Van Gameren; Weikmans; Zaccai, 2014). 

No entanto, outras soluções podem se tornar necessárias quando as mudanças 
climáticas se tornarem significativas ou quando uma consideração de longo prazo se tornar 
necessária. Neste caso, a revisão das vulnerabilidades existentes torna-se insuficiente 
para construir uma estratégia de adaptação ou mitigação e medidas específicas devem ser 
postas em prática para lidar com novos impactos potenciais. Este é o caso, por exemplo, de 
projetos de desenvolvimento em larga escala que criam irreversibilidade a longo prazo. A 
curto prazo, é útil considerar primeiro a origem e a evolução da vulnerabilidade. O aumento 
do risco que estamos observando atualmente e que explica o aumento das perdas ligadas 
aos desastres naturais tem causas específicas, ligadas aos atuais desenvolvimentos 
socioeconômicos e demográficos (migração para áreas de risco, escassez de terra, custo 
econômico e político de medidas de redução de risco, consumo excessivo de água, etc.). 
Ao invés de simplesmente procurar reduzir os riscos, é necessário considerar os fatores 


que explicam o aumento dos riscos e atacar suas causas originais. 


2.4 Planos estratégicos de adaptação e mitigação da mudança climática 


Uma das formas pelas quais as autoridades públicas podem integrar a mudança 
climática é desenvolvendo estratégias e/ou planos de adaptação destinados a reduzir 
a vulnerabilidade (de um país, região ou cidade) à mudança climática, inclusive à 
variabilidade climática e aos extremos. Com o apoio de diferentes medidas e políticas, 
estes instrumentos de coordenação podem construir uma visão global para todo o território 
e os setores envolvidos ou focar em certas vulnerabilidades mais específicas. Este é um 
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momento crucial para planejar e tornar as cidades mais inteligentes com o objetivo de utilizar 
e inserir novas tecnologias de informação e comunicação nos diferentes setores, a fim de 
otimizar o uso das infra-estruturas existentes para enfrentar o novo desafio global. Seja no 
transporte, construção civil, governança ou meio ambiente e energia, as novas tecnologias 
podem ajudar a enfrentar os atuais desafios urbanos que estão aumentando os gases de 
efeito estufa. Com base nas estratégias de 25 cidades americanas que deverão reduzir 74 
milhões de toneladas métricas de emissões de carbono até 2030, enquanto cumprem ou 
excedem as metas estabelecidas no Acordo de Paris de 2015 em meados da década 2010 
a 2020, de acordo com o relatório publicado pela Bloomberg Philantropies. As 25 cidades 
aprovaram um total de 54 grandes políticas verdes que podem reduzir nossa pegada de 
carbono coletiva: Construção Civil, Transporte e Energia (Udasain, 2021). Outras ações 
podem ser acrescentadas, como o desenvolvimento do teletrabalho e a implementação de 
cidades esponjosas. 

Assim, segundo Egli, Hasler & Probst (2016), as abordagens e estratégias para 
responder às mudanças climáticas e para iniciar e implementar medidas podem ser 
agrupadas em três grandes categorias. 

1. A estratégia de eficiência, frequentemente chamada de estratégia de eco- 
eficiência, visa uma maior produtividade dos recursos, ou seja, o uso de um 
mínimo de recursos por unidade de produção. Ela também procura minimizar 
as emissões de gases de efeito estufa, em particular de CO2. Esta abordagem 
é favorecida no mundo econômico, pois visa e permite a introdução de novos 
produtos e processos. Um exemplo clássico desta estratégia é o desenvolvimento 
da mobilidade inteligente (carros elétricos, bicicletas elétricas, etc.). Em uma 
abordagem eco-eficiente, o uso de energia renovável e de matérias-primas 
produzidas de forma sustentável é garantido. As inovações que limitam as emissões 
de CO2 e controlam a qualidade do ar são encorajadas, tanto na produção de um 
bem como em sua utilização. 


2. A estratégia de suficiência se baseia na suposição de que é possível limitar o 
consumo de recursos e a degradação ambiental sem reduzir a qualidade de vida. A 
qualidade de vida poderia até mesmo ser melhorada consumindo menos. O slogan 
“qualidade de vida em vez de crescimento econômico” poderia resumir a visão dos 
defensores desta abordagem. 


3. A estratégia de coerência visa limitar o uso dos recursos ao que a natureza 
pode renovar. O conceito de reciclagem está relacionado a esta estratégia, e 
os biocombustíveis são outro exemplo. A idéia é que no final de um processo 
praticamente só há matéria-prima para novos produtos, e nenhum desperdício. 
Esta abordagem implica o envolvimento de produtores e consumidores na busca 
de soluções relevantes. 
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31 MÉTODO DE PESQUISA 


Este artigo faz parte de um paradigma construtivo de adaptação e mitigação da 
mudança climática. É o resultado de uma pesquisa bibliográfica dentro da disciplina de 
Inovação e Nova Dinâmica Social e Econômica em Cidades: Cidades Inteligentes e Cenários 
de Inovação Urbana do programa de pós-graduação em Desenvolvimento Regional 
e Meio Ambiente da Universidade do Estado de Minas Gerais (UEMG). As informações 
que formam a base deste trabalho provêm dos mais recentes artigos científicos, teses e 
relatórios sobre o cenário climático, a fim de produzir um documento técnico-sintético que 


atenda às normas da comissão científica. 


41 RESULTADOS E DISCUSSÕES 


4.1 Tecnologias e práticas de mitigação na construção civil 


É bem conhecido que os combustíveis fósseis consumidos pelo setor de construção 
(para a operação de edifícios e a construção de novos) é a principal fonte de emissões de 
GEE. O Quarto Relatório de Avaliação Intergovernamental do IPCC identifica as principais 
fontes de emissões de GEE associadas aos edifícios da seguinte forma: aquecimento de 
salas, refrigeração de salas, aquecimento de água, iluminação artificial e uso de aparelhos 
(Levine et al., 2007). O relatório também identifica três categorias de medidas para reduzir 
as emissões de GEE dos edifícios: 


1. Reduzir o consumo de energia dos edifícios e sua energia incorporada 


2. Usar mais energia com baixo teor de carbono e promover uma maior participação 


de energia renovável 


3. Controle de emissões de GEE não CO, (Levine et al., 2007) 


Neste contexto, os enfoques de mitigação no setor da construção civil são sobre a 
implantação e concepção de estratégias, tecnologias e práticas que Reduzir a demanda 
e o consumo de energia associados aos edifícios - desde o projeto, construção, descarte, 
até a renovação e o fim de vida útil. Avançar em direção a uma energia de baixo carbono. 
Maximizar as oportunidades de sequestro de carbono nos edifícios. Promover mudanças 
comportamentais para alcançar um estilo de vida sustentável. Tudo isso resume que os 
edifícios precisam ser mais inteligentes, utilizando tecnologias de ponta. 


4.2 Mitigação na mobilidade inteligente 


A contribuição do setor de transportes para as emissões de gases de efeiti estufa 
tem sido amplamente demonstrada, e o importante papel que ele pode desempenhar em 
uma estratégia eficaz e eficiente para reduzir as emissões globais de GEE é cada vez mais 
compreendido. Como resultado, o transporte é um ativo chave na transição energética e 
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ambiental, com mais de 53% do consumo global de petróleo atendendo a 94% da demanda 
energética do setor em 2010. O setor de transportes é uma área chave para questões 
de segurança energética e uma importante fonte de poluentes atmosféricos, tais como 
ozônio, óxido de nitrogênio e material particulado, bem como (Universidade de Cambridge, 
2015) dióxido de carbono. A luta contra a mudança climática requer uma ação urgente e 
sustentável. Neste contexto, a mobilidade é um setor prioritário na redução dos gases de 
efeito estufa (GEE). Esta é uma das principais razões que o mundo tem para mudar de 
veículos a gasolina para veículos elétricos. Teremos inevitavelmente uma mobilidade mais 
limpa e ambientalmente mais amigável. Além disso, as 25cidades americanas que estão 
implantando mais de 11.000 veículos elétricos no “American Cities Climate Challenge * 
estão demonstrando a eficácia dos veículos elétricos na luta contra a mudança climática 
porque, ao seguirem esta rota, ultrapassarão de longe os objetivos estabelecidos pelo 
Acordo de Paris até 2025 (Udasain, 2021). 

Outro aspecto chave da estratégia de combate à mudança climática é o ciclismo, que 
é uma das maiores esperanças da humanidade. As 25 cidades americanas desenvolveram 
mais de 510 milhas de ciclovias, e cada viagem adicional feita por bicicleta reduz as 
emissões de CO2 em 14%. Substituir o carro por um passeio a pé ou de bicicleta, mesmo 
apenas um dia por semana, pode reduzir nossa pegada de carbono em até meia tonelada 
de CO? ao longo de um ano. Ao fazer isso nas cidades, podemos reduzir significativamente 
nossa pegada de carbono coletivamente. 


4.3 Desenvolvimento do teletrabalho 


Já vimos anteriormente que o setor de transporte é o principal emissor de gases 
de efeito estufa, grande parte dos quais se deve a viagens entre a casa e o trabalho. 
O desenvolvimento do teletrabalho é outra arma fundamental na luta contra a mudança 
climática. Este pressuposto de um impacto ambiental positivo do teletrabalho é parte de 
uma crença mais geral na capacidade da tecnologia através de Wi-Fi ubíquo e dispositivos 
móveis como smartphones, computadores ou tablets para resolver problemas ambientais 
por si só (Trambley e al., 2007). 


4.4 Ambiente inteligente 


A gestão da água, do lixo e da energia estão no centro das preocupações ambientais 
de uma cidade. Em uma cidade inteligente, várias ferramentas tecnológicas permitem 
a proteção e preservação de nossos recursos naturais e ambientes naturais, tais como 
sensores para medir os níveis de poluição do ar. Estas são novas tecnologias que fornecem 
uma gama de informações em tempo real. Em energia, as Smarts Grids, uma tecnologia 
informática para redes inteligentes de distribuição de eletricidade, podem otimizar a 
produção e distribuição de eletricidade enquanto se ajustam à demanda. Economizar 
energia através de novas tecnologias também é um ambiente inteligente. Equipar a infra- 
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estrutura da cidade e implementar novas tecnologias de informação e comunicação no 
campo do meio ambiente visa a proteger o meio ambiente, garantir o uso sustentável dos 
recursos e melhorar o meio ambiente natural (Simard, 2015). 


4.5 Desenvolvimento e implementação de cidades esponjosas 


Inevitavelmente, da forma como o planeta Terra está aquecendo, algumas regiões 
têm que desenvolver estratégias eficazes para se adaptarem à nova mudança agora. 
Por exemplo, conforme a temperatura aumenta com o aquecimento global, mais e mais 
umidade evapora para a atmosfera, resultando em chuvas mais pesadas e violentas em 
todo o mundo. A implementação de mecanismos de mitigação não poderá impedir o mundo 
de enfrentar esses momentos de chuvas fortes. Mas diante desta situação, estratégias de 
adaptação podem ser desenvolvidas, tais como o estabelecimento de cidades esponjosas. 
Este conceito foi concebido pelo eminente professor Yu Kongjian (Tessa, 2021), um dos 
principais planejadores urbanos da China, com o objetivo de gerenciar as inundações 
nas cidades chinesas. O professor propõe não construir nas áreas baixas e deixar esta 
parte como uma espécie de reservatório com lagos naturais. Assim, a partir daí, em vez 
de canalizar a água rapidamente em linhas retas, haverá rios sinuosos com vegetação ou 
planícies de inundação que retardarão a água. Outra vantagem da cidade esponjosa é a 
criação de espaços verdes, parques e habitats animais, enquanto purifica a água superficial 
com plantas que removem toxinas e nutrientes poluentes. 

No final, adaptação e mitigação são duas estratégias que andam de mãos dadas. 
A mitigação reduz a gama de riscos futuros. A adaptação nos permite lidar com mudanças 
reais, é óbvio que algumas mudanças ainda ocorrerão, portanto teremos que estar 
mais preparados para aproveitá-las ao máximo. Especialistas dizem que as medidas de 
adaptação e mitigação são apoiadas por fatores de capacitação: instituições e governança 
eficazes, inovação e investimento em tecnologias e infra-estrutura verdes, meios de vida 
sustentáveis, estilos de vida e comportamento. Em resumo, quando agimos para reduzir 
os GEE, também nos capacitamos para sermos mais eficazes na adaptação. Para que 
a adaptação e mitigação sejam eficazes, precisamos de padrões de planejamento de 
uso do solo (construção e aumento das ciclovias, instalação de sensores para monitorar 
a qualidade do ar, etc.) e padrões de construção inteligentes, ou seja, com tecnologia 
podemos otimizar os edifícios para reduzir nossa pegada de carbono e, finalmente, a 
transição para energia sustentável. 


4.6 Limitações e perspectivas das políticas nacionais e subnacionais de 
adaptação e mitigação 

As perspectivas de ação pública sobre adaptação e/ou mitigação da mudança 
climática dizem respeito a diferentes instrumentos políticos através da intervenção do 
Estado. Mas os resultados das práticas atuais são muito limitados quando comparados aos 
impactos projetados. De acordo com Biesbroeck et al. (2011), existem sete (7) categorias 
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de barreiras. Primeiro, escalas temporais de conflito entre os impactos de longo prazo da 
mudança climática e a visão de curto prazo da política e dos processos de tomada de decisão. 
Segundo, as incertezas das complexidades nas decisões políticas (incertezas subjetivas, 
estratégicas e institucionais). Uma terceira categoria de barreiras é a falta ou, pelo menos, 
o excesso de instituições. Quarto, a fragmentação ou falta de coordenação entre políticas 
e instituições em diferentes níveis de governança é problemática, especialmente em um 
campo tão multissetorial quanto a adaptação. Um baixo nível de consciência e comunicação 
sobre o tema da adaptação e/ou mitigação é outro obstáculo ao desenvolvimento de 
políticas, assim como a falta de motivação e vontade de agir, ligada, entre outras coisas, 
aos determinantes psicológicos e cognitivos. Finalmente, a falta de recursos humanos, 


financeiros e tecnológicos é obviamente uma grande limitação operacional. 


51 CONSIDERAÇÕES FINAIS 


O estudo do clima usando conjuntos de modelos climáticos é um método semelhante 
ao usado pelo GIEC na preparação de seus relatórios de avaliação sobre a mudança 
climática. Assim, a parametrização do modelo para estudar o impacto da mudança climática 
se revela uma verdade absoluta sobre este fenômeno climático e já está causando mudanças 
desproporcionais no clima. Já estamos vendo os efeitos, com eventos ocorrendo em uma 
velocidade e magnitude maiores do que o esperado. A mudança climática está afetando o 
meio ambiente, o planeta Terra e os seres humanos que o habitam de muitas maneiras. O 
mais preocupante é que cada grau de aumento na temperatura aumenta o risco de eventos 
climáticos mais destrutivos, causando danos significativos e possivelmente irreversíveis 
em todo o mundo. É necessária uma ação urgente para mitigar o aumento dos gases de 
efeito estufa e para se adaptar às mudanças que ainda terão que ocorrer. Há um consenso 
internacional de que a ação sobre a mudança climática requer uma abordagem dupla 
com o objetivo de reduzir as emissões de gases de efeito estufa (mitigação) e reduzir 
a vulnerabilidade dos sistemas naturais e humanos aos impactos da mudança climática 
(antecipação e adaptação). Adaptação e mitigação são duas (2) estratégias para lidar com 
o aquecimento global e como as cidades são fontes de emissão de gases de efeito estufa, 
elas estão na linha de frente para agir usando ferramentas tecnológicas para otimizar 
os edifícios tornando-os mais inteligentes, para assumir o controle da qualidade do ar 
com a instalação de sensores, Incentivar o ciclismo, aumentando o número de ciclovias, 
favorecendo os carros elétricos em vez dos carros a gasolina e limitando as viagens das 
pessoas ao trabalho através do desenvolvimento do teletrabalho. Na verdade, quanto mais 
inteligente for uma cidade, melhor ela é capaz de lidar com o aquecimento global. Como 
a mudança climática é um fenômeno global, as ações de adaptação e mitigação devem 
ocorrer em todo o mundo. Portanto, será necessário aumentar a resiliência e a inteligência 


de todas as cidades do mundo. 
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RESUMO: Esteroides são compostos 
orgânicos classificados como hormônios 
e são produzidos pelas diversas espécies 
do reino animal, vegetal e fungi protista. 
São moléculas formadas por cadeias 
carbônicas tetracíclicas, lipossolúveis, 
estáveis e apresentam diferentes tipos 
de grupos funcionais (álcoois, ésteres, 


cetonas, aldeídos, ácidos carboxílicos, 
núcleo 'benzenoide) os quais são 
responsáveis pelas suas diferentes 


funções e propriedades físicas e químicas. 
Dentre os esteroides, a progesterona é um 
hormônio sexual produzido pelos ovários 
e é fundamental para a regulação do ciclo 
menstrual e para a manutenção da gravidez. 
Com os estudos e o desenvolvimento dos 
contraceptivos, observou-se que pequenas 
variações, nas estruturas e nas funções 
orgânicas dessas moléculas, podem levar a 
grandes alterações da biodisponibilidade e 
de suas atividades biológicas. O consumo 
de contraceptivos, remete posteriormente 
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na liberação dessas moléculas em corpos hídricos, através do descarte ou da eliminação 
natural humana. Por serem moléculas lipofílicas e recalcitrantes, essas tendem a se 
acumular nos tecidos de animais aquáticos (exs. peixes, anfíbios, etc.) causando danos em 
seus sistemas reprodutivos, o que às classificam como poluentes disruptores endócrinos. 
Uma das maneiras de se promover alterações nessas moléculas no meio ambiente se dá 
através de reações naturais, catalisadas pelos sistemas enzimáticos de fungos e bactérias. 
Assim, este Capítulo apresenta os estudos de biotransformação da progesterona com 
o fungo de ambiente marinho Penicillium oxalicum CBMAI 1996. Este fungo promoveu a 
hidroxilação da progesterona, formando os derivados 15P-hidroxiprogesterona e a 7,15P-di- 
hidroxiprogesterona, o quais foram isolados, identificados e caracterizados por CG-EM, RMN 
e por difração de Raios-X. Comumente, enzimas do complexo do citocromo P-450, produzidas 
pelo fungo P. oxalicum CBMAI 1996, foram as que promoveram as reações de hidroxilação 
da progesterona. Além disso, há uma hipótese que tais moléculas esteroidais, presentes no 
meio ambiente terrestre e marinho, estão sujeitas às diversas biotransformações microbianas. 
Assim, podem ocorrer alterações na lipofilicidade, e consequentemente, de bioacumulação 
podendo levar alterações metabólicas em organismos vivos alvos que passam a ter contato 
com esses compostos xenobióticos e seus derivados. Desta forma, este trabalho apresenta 
alguns dos estudos recentes de biotransformação da progesterona com o fungo. P oxalicum 
CBMAI 1996. 

PALAVRAS-CHAVE: Biotransformação, Esteroides, Fungo Marinho, Penicilium oxalicum. 


RECENT STUDIES OF THE BIOTRANSFORMATION OF PROGESTERONE BY 
MICROORGANISMS 


ABSTRACT: Steroids are organic compounds classified as hormones and are produced by 
different species of the animal, vegetable and protist fungi kingdoms. They are molecules 
formed by tetracyclic carbon chains, liposoluble, stable and have different types of functional 
groups (alcohols, esters, ketones, aldehydes, carboxylic acids, benzene ring) which are 
responsible for their different functions and physical and chemical properties. Among steroids, 
progesterone is a sex hormone produced by the ovaries and is essential for regulating the 
menstrual cycle and maintaining of the pregnancy. With the studies and development of 
contraceptives, it was observed that small variations in the structures and organic functions 
of these molecules can lead to major changes in bioavailability and their biological activities. 
The consumption of contraceptives subsequently refers to the release of these molecules 
into water bodies, through disposal or natural human elimination. Because they are lipophilic 
and recalcitrant molecules, they tend to accumulate in the tissues of aquatic animals (e.g. 
fish, amphibians, etc.) causing damage to their reproductive systems, which classify them 
as endocrine-disrupting pollutants. One of the ways to promote changes in these molecules 
in the environment is through natural reactions, catalyzed by the enzymatic systems of fungi 
and bacteria. Thus, this Chapter presents the studies of progesterone biotransformation 
with the marine fungus Penicillium oxalicum CBMAI 1996. This fungus promoted the 
hydroxylation of progesterone, forming the derivatives 15B-hydroxyprogesterone and 7P,15B- 
dihydroxyprogesterone, which were isolated, identified and characterized by GC-MS, NMR 
and X-ray diffraction. Commonly, enzymes of the cytochrome P-450 complex, produced by 
the fungus P. oxalicum CBMAI 1996, were the ones that promoted progesterone hydroxylation 
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reactions. In addition, there is a hypothesis that such steroidal molecules, present in the 
terrestrial and marine environment, are subject to various microbial biotransformation. Thus, 
changes in lipophilicity and, consequently, bioaccumulation may occur, which may lead to 
metabolic changes in target living organisms that come into contact with these xenobiotic 
compounds and their derivatives. Thus, this work presents some of the recent studies of 
progesterone biotransformation by the fungus. P oxalicum CBMAI 1996. 

KEYWORDS: Biotransformation, Steroids, Marine Fungi, Penicillium oxalicum. 


11 INTRODUÇÃO 


1.1 Progesterona 


A progesterona é um hormônio esteroidal, lipofílico, derivado do colesterol, 
caracteriza-se por ser o único progestágeno natural, ou seja, é sintetizada pelo corpo 
humano, especialmente nos ovários, na placenta e nas glândulas suprarrenais!!, Sua 
estrutura química é formada a partir de um núcleo base composto por 17 átomos de 
carbono, que compõe um sistema de anéis tetracíclicos (A, B, C e D), denominado de 
gonano, típico dos esteroides (Figura 1). 


H O 


5 6 


4 
(A) Núcleo Gonano (B) Progesterona 


Figura 1. (A) Representação do esqueleto gonano dos esteroides. (B) Representação da estrutura 
química da progesterona. 


A biossíntese da progesterona envolve o citocromo P-450, nas células do ovário, 
sob a influencia dos hormônios hipotalâmico-hipofisiários, com a conversão do colesterol 
à pregnenolona. Esta por sua vez, deve sair da mitocondria para sofrer a ação das 
enzimas 38HSD ou 17a-hidroxipregnenolona e 17a-hidroxilase, presentes no retículo 
endoplasmático, convertendo-se finalmente à progesterona. A maior parte da progesterona 
formada é secretada para o sangue, porém uma outra parte é oxidada no carbono C17, 
tornando-se um substrato de partida (17a-hidroxiprogesterona) para a formação de outros 
hormônios esteroidais (Figura 2)!I. 

A progesterona é um dos hormônios mais importantes da gestação e comunicação 
materno-fetal, estimulando a secreção uterina de nutrientes que influenciam diretamente 
no início do desenvolvimento embrionário!!, Ainda, a progesterona está relacionada ao 
processo de ovulação, e com isso, diversos progestágenos vem sendo (bio)sintetizados 
e comercializados, em combinação com outros hormônios como o estradiol, na forma de 
contraceptivos!'], 
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Durante todo o período menstrual, a expressão dos receptores endometriais para 
estrogênios e progesterona varia, o que resulta na ciclicidade e formação das condições 
para a concepção. A progesterona exerce um papel importante na transição da fase 
proliferativa, ao final da menstruação, quando o endométrio está descamado, com sua 
menor espessura, para a fase secretora, período fértil, quando o endométrio apresenta 
sua maior espessura, rico em nutrientes favorecendo a fixação e o sustento do embrião!º, 


A produção de progesterona nos ovários é induzida por hormônios produzidos 
no hipotálamo e pela glândula hipófise, na segunda metade do ciclo menstrual, antes 
da ovulação. O hormônio é capaz de ligar-se aos receptores de estrogênio, presentes 
nas células epiteliais do estroma e miométrio, de maneira supressora. Desta forma, na 
fase ovulatória, quando se observa uma alta liberação de estrogênio, tem-se de forma 
simultânea um aumento exponencial da concentração de receptores da progesterona, 
localizados nas células epiteliais do endométrio. Isso faz com que a progesterona esteja 
menos biodisponível e não realize sua função supressora. Durante a fase lútea o nível de 
progesterona no sangue encontra-se em torno de 20 a 30 mg/dia, enquanto fora desta fase 


fica próximo a 1mg/dia!i. 
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Figura 2. Biossíntese de hormônios esteroidais a partir das células do ovário. 


Fonte: Adaptado da ref. ORIZABA-CHÁVEZ, B; et. al. Farmacocinética de la progesterona. Revista 
Hospital Juarez México, México, vol. 60, 1, 59-66, 2013. 
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A progesterona é metabolizada de forma rápida, em sua primeira passagem pelo 
fígado. O principal metabólito da progesterona é o pregnanediol (Figura 3), que em 
sequência se conjuga com o ácido glicurônico, formando o pregnanediol glicuronídeo. Esse 
metabólito é então eliminado pela urina. Como esse processo ocorre na primeira passagem 
hepática, a progesterona possui um curto tempo de meia-vida de eliminação de apenas 5 
minutos. Seus derivados, em sua maioria, possuem um tempo de meia vida de eliminação 
muito maior, como o acetato de medroxiprogesterona (Figura 4) que chega a 24 horas, e 


por isso são mais utilizados na composição dos contraceptivos!!. 


HO 
UDP “ 
fe 
ane EE 
transferase) Ho À 
Ácido glicurônico Pregnanediol Pregnanediol glicuronídeo 


Figura 3. Reação de conjugação do pregnanediol glicuronídeo nas células hepáticas!º). 


Além do baixo tempo de meia vida de eliminação da progesterona, outros fatores 
levaram ao desenvolvimento de derivados esteroidais para aplicação em métodos 
contraceptivos. As progestinas são moléculas sintetizadas a partir da progesterona, sendo 
estas derivadas da 17a-hidroxiprogesterona e norprogesterona. Estas possuem atividades 
semelhantes à progesterona e podem ser administradas de forma terapêutica por via 
oral“, Estão dispostas na Tabela 1 as estruturas de algumas progestinas utilizadas em 
tratamentos contraceptivos e seu ano de introdução do mercado!”. 

As progestinas também estão sendo empregadas em terapias hormonais para 
mulheres que sofrem com os efeitos da menopausa. Estudos localizaram receptores 
de progesterona em osteoblastos (células onde se formam os óssos). A partir de testes 
in vitro verificou-se o efeito benéfico promovido pelas terapias com progestágenos, que 
aumentaram a proliferação de osteoblastos, levando a uma diminuição da perda óssea 
(consequência da menopausa), 
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Nome Estrutura Ano de introdução no mercado 


Clormadinona 1959 


Medroxiprogesterona 1957 


Medrogestona 1964 


Tabela 1. Classificação dos compostos derivados da progesterona utilizados nos tratamentos 
contraceptivos e seu ano de introdução no mercado. 


Fonte: REGIDOR, P. A. The clinical relevance of progestogens in hormonal contraception: Present 
status and future developments. Oncotarget, Ismaning, Germany, vol. 9, 77, 34628-34638, 2018. 
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1.2 Contaminação ambiental por compostos esteroidais 


O consumo dos contraceptivos hormonais cresce progressivamente desde a década 
de 1960. Segundo a declaração da posição das terapias hormonais de 2022, da Sociedade 
Norte Americana da Menopausa (The North American Society of Menopause), além do uso 
para prevenção da concepção, os hormônios apresentam resultados altamente eficazes 
para o tratamento de sintomas vasomotores, da síndrome geniturinária e da perda óssea, 
sintomas da menopausat!!, 

Dados da literatura tem mostrado, assim como a maioria dos fármacos, que apenas 
uma pequena parte do princípio ativo dos anticoncepcionais e reguladores hormonais 
é absorvida pelo organismo, enquanto uma outra parcela é metabolizada ou eliminada 
diretamente na urina. Estudos apontaram a presença em águas subterrâneas de até 66 ng 
L! de 17B-estradiol e até 50 ug L'! de 17a-etinilestradiol (ambos esteroides presentes em 
contraceptivos) 19, 

As águas residuais domésticas e industriais passaram a apresentar concentrações 
consideráveis de compostos hormonais e seus derivados de metabolização. Tal situação 
torna-se crítica e preocupante em função do forte potencial disruptivo endócrino desses 
esteroides, e da ausência de tecnologias eficientes para a remoção e a identificação de 
moléculas orgânicas em baixas concentrações nas estações de tratamento de águas 
residuais!" 

Devido a incipiência dos métodos de tratamentos de água e esgoto frente a tais 
substâncias, essas acabam influenciando o comportamento dos organismos vivos em vários 
aspectos. Estudos mostraram que peixes fêmeas, na presença de determinados compostos 
esteroidais, passaram a apresentar aspectos andrógenos, ou seja, desenvolvimento 
induzido de órgãos reprodutivos masculinos, mesmo quando presentes em concentrações 
mínimas (valores inferiores a ng/L)!'2, 

Foram realizados testes com peixes alevinos da espécie Oreochromis niloticus 
(tilápia do Nilo), os quais foram expostos por um período de 1 a 4 semanas aos hormônios 
esteroidais 17B-estradiol e 17a-etinilestradiol. As concentrações de esteroides avaliadas 
variavam de 250 ng L'! a 1 mg L'!. Concluiu-se que as concentrações hormonais do meio 
influenciaram diretamente no peso, comprimento e na sobrevivência dos indivíduos. 
Observou-se o desenvolvimento de indivíduos intersexos e com malformações, 
principalmente na região da cabeça ou encefálica!'?, 

O Fator de Bioconcentração - FBC (Bioconcentration Factor) trata-se de um número 
que remete a afinidade que um determinado composto apresenta em se bioacumular em 
tecidos animais ou vegetais!'!, Mediu-se o FBC de progestágenos (grupo de hormônios 
sexuais que atuam nos receptores da progesterona) em indivíduos da espécie Danio 
rerio (peixes-zebra) que possuíam de 9 a 12 meses de vida. Os testes foram realizados 
com um período de exposição aos esteroides de 28 dias. Os resultados apresentaram 
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um intervalo de FBC = 7, relativo ao hormônio contraceptivo Dienogeste, até FBC = 
128, relativo ao hormônio Acetato de Medroxiprogesterona. Avaliou-se, que a exposição 
aos hormônios levou alterações bioquímicas, gerando respostas não lineares dos 
biomarcadores analisados. Além disso observou-se uma intensa alteração na produção 
de proteínas do grupo das vitelogeninas, sugerindo um risco para o processo reprodutivo 
destes indivíduos !'SI, 

Os esteroides liberados no meio ambiente tendem a comportar-se de acordo com 
suas propriedades físico-químicas as quais são hidrofóbicos, apresentam baixa volatilidade 
e são bastante estáveis e recalcitrantes 18, Por consequência, causam acumulação em 
sedimentos ambientais, bem como em tecidos adiposos de animais e outros organismos 
presentes na biota local os quais têm contato com essas substâncias. 

Com o decorrer do tempo pode-se observar um processo de bioacumulação, 
em que essas substâncias esteroidais se concentram progressivamente nos tecidos 
dos organismos a eles expostos. Simultaneamente a esse processo, desenrolam-se os 
ciclos naturais predatórios das cadeias alimentares, nos quais surgem as tendências dos 
acúmulos dos contaminantes nos tecidos dos organismos consumidores do topo da cadeia 
alimentar via processo de biomagnificação!!?, 

Em contrapartida, os microrganismos presentes nos ambientes de disseminação 
de tais substâncias interagem com estes compostos podendo biotransformá-los em novas 
moléculas, comumente mais hidrofílicas. O etinil estradiol, por exemplo, é parcialmente 
metabolizado no corpo humano gerando derivados conjugados, os quais são regenerados 
por bactérias E. coli retornando a sua estrutura originallal, 

A contaminação ambiental por esteroides é emergente, sendo relevantes os estudos 
da ocorrência e a concentração de tais compostos e de seus derivados de biotransformação 
(mesmo que em escalas de ng L?) 9. Além disso, deve-se buscar compreender como a 
biota local interage com os esteroides, para se avaliar as possíveis formas de biodegradação 
ou a obtenção de novos derivados e as suas possíveis consequências. Contudo, é um 
trabalho complexo, de difícil interpretação e de conclusão das hipóteses, devido às 
inúmeras variáveis que envolvem todo o processo de biotransformação dos esteroides no 
meio ambiente. 

Ressalta-se que os produtos de biotransformação dos esteroides, que estão sendo 
obtidos no meio ambiente, podem ser classificados como mais ou menos tóxicos, ou 
ainda, como novos compostos com potencial atividade para uso em humanos. Diversas 
moléculas comerciais foram obtidas a partir do metabolismo natural de plantas, animais 
e fungos, ou por processos de biotransformação??, Destaca-se, com isso, a importância 
e a relevância dos estudos de biotransformação de esteroides. 
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1.3 Biotransformação de esteroides 


Durante o desenvolvimento e o crescimento de microrganismos (fungos, bactérias, 
algas) no meio ambiente ou em condições de laboratórios são biossintetizadas uma 
diversidade de enzimas, as quais podem ser excretadas, isoladas, imobilizadas e utilizadas 
para promover a biotransformação e a biodegradação de compostos xenobióticos. Em 
estudos laboratoriais, para promover a biotransformação de compostos esteroidais, 
conhecer os metabólitos produzidos e as enzimas envolvidas, podem ser utilizados os 
sistemas biológicos constituídos de células totais ou o caldo enzimático obtido dos meios 
de cultura em que foram cultivados e crescidos os microrganismos 2", 

Os microrganismos empregados nas biotransformações encontram-se distribuídos 
no meio ambiente em todos os tipos de biomas. Consequentemente, as condições 
empregadas nos meios reacionais para a biotransformação de xenobióticos são brandas 
e próximas das CNTP, além do meio líquido de cultivo laboratorial utilizar a água e sais 
minerais. Com isso, estudos nessa área convergem com os princípios da Química Verde, 
uma vez que levam ao desenvolvimento de técnicas e metodologias que exigem um baixo 
consumo de energia e reduzem o uso de substâncias nocivas à saúde humana e ao meio 
ambiente e que sejam de fontes renováveis 22, 

A biotransformação remete a qualquer mudança promovida em moléculas na 
presença de um agente biológico vivo, a qual leva às alterações nas propriedades físicas e 
químicas do material de partida, cujo produto(s) pode(m) ser extraído(s) após o processo. A 
biodegradação, em contrapartida, remete às transformações que causam a fragmentação 
total dos compostos (comumente de interesse os compostos xenobióticos tóxicos) em 
estruturas químicas mais simples como gás carbônico, água e outras moléculas voláteis 2º, 

Estudos na literatura, por exemplo, mostraram que o fungo Aspergillus niger 
N402 promoveu a obtenção de uma mistura de produtos via reações de hidroxilação da 
progesterona, produzindo a 11a-hidroxiprogesterona, a 6P,11a-di-hidroxiprogesterona e a 
21-hidroxiprogesterona (desoxicorticosterona). A adição de íons metálicos (Cu?*, Cd? e 
Co?*) aos meios reacionais levou a inibição das enzimas que promoveram a hidroxilação, o 
que remete na possibilidade de se controlar a biotransformação!!, 

A biotransformação da progesterona pelo fungo Aspergillus brasiliensis por 7 
dias levou a obtenção da 11a-hidroxiprogesterona, da 14a-hidroxiprogesterona e da 
21-hidroxiprogesterona. Ressalta-se que a 21-hidroxiprogesterona é utilizada como 
fármaco, na forma do seu derivado acetato, para o tratamento de doenças de Addison, 
Miastenia Gravis, dentre outras! 

As monoxigenases promovem reações de oxidação (ex. hidroxilação) dos 
esteroides, gerando uma vasta gama de produtos, sendo estes muitas vezes regioisômeros. 
A consequência disso, é que muitas vezes, dificulta a purificação do(s) produto(s) alvo(s). 
Nestes casos, pode-se utilizar microrganismos como plataformas para a expressão 
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heteróloga das enzimas. Por exemplo, empregar células microbianas hospedeiras, como 
de bactérias (Escherichia coli) e de leveduras (Saccharomyces cerevisiae, Pichia pastoris). 
Assim, obtêm-se as enzimas fúngicas de interesse com um alto nível de expressão e com 
melhor rendimento de biotransformação frente ao(s) produto(s) desejado(s)!º]. 

A origem dos microrganismos que podem ser utilizados nas reações de 
biotransformação é altamente diversa, já que estes ocupam todos os ambientes, 
independente das condições de temperatura, pressão e composição nutricional. Fungos 
endofíticos, por exemplo, vivem dentro dos tecidos vegetais, sem causar danos à planta 
hospedeira. O fungo endofítico Talaromyces sp., isolado das folhas de Handroanthus 
impetiginosus (Ipê-roxo) reduziu a ligação dupla sp? dos carbonos C-4 e C-5 (Anel A), 
além de promover oxidação no carbono C-17 com hibridação spº (Anel C) da progesterona 
produzindo derivados diferentes!?, 

Na literatura foi avaliado o potencial de biotransformação da progesterona com 
microrganismos presentes no lodo ativado de tratamento de efluentes. Observou-se várias 
reações (metilação, hidrogenação, hidroxilação, etc.), com uma amostra microbiológica do 
efluente. Neste estudo, realizou-se o teste do efluente resultante após a biotransformação 
da progesterona com peixes-zebra. Embora não tenha se observado a morte de indivíduos, 
verificou-se interferências na atividade de transcrição gênica evidenciando que os resíduos 
da biotransformação continuaram possuindo atividade endócrina e biológica, possivelmente 
devido à presença dos esteroides biotransformados pela microbiota!), 

O potencial potencial de biotransformação da progesterona por 12 linhagens de 
fungos isolados de ambiente marinho foram investigados entre 3 a 7 dias. Os melhores 
resultados foram obtidos após 7 dias, cuja conversão da progesterona foi superior a 99%, 
em diferetnes produtos de biotransformação. Tais resultados foram promovidos pelas 
linhagens de Aspergillus sydowii CBMAI 935, que produziu a testosterona, o Penicillium 
oxalicum CBMAI 1185, que produziu a testololactona e o Cladosporium sp. CBMAI 1237, 
que produziu a 11-hidroxiprogesterona!??l, 


21 BIOTRANSFORMAÇÃO DA PROGESTERONA PELO FUNGO PENICILLIUM 
OXALICUM CBMAI 1996 

Em nossos estudos foi avaliado o fungo o Penicillium oxalicum CBMAI 1996, frente 
às reações de biotransformação da progesterona. Esse fungo foi isolado de uma esponja 
(Chelonaplysilla erecta) coleta do litoral norte do estado de São Paulo, O objetivo de 
estudar os fungos de ambiente marinho em biotransformações, pois o mar é um dos 
destinos finais dos esteroides e muitas outras moléculas. 

Para realização da reação, primeiramente cultivou-se o fungo P oxalicum CBMAI 
1996 em placas de Petri contendo o meio de cultura sólido de extrato de malte esterilizado, 
por 7 dias à temperatura de 32ºC. O meio de cultura foi preparado a partir de uma solução 
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de água do mar sintética composta de sais [CaCI,.2H,0 (1,36 g/L), MgCI,.6H,0 (9,68 g/L) 
NaCl (309/L), NaH,PO, (1,4x10“ g/L), Na,SO, (3,47 g/L), NaHCO, (0,17 g/L), KBr (0,10 
g/L), SrCI,.6H,0 (0,04g/L) e do ácido HBO, (0,03 g/L)]. Além destes componentes, foi 
adicionado o extrato de malte (20 g/L), agar (20 g/L) e o pH do meio foi ajustado (pH = 7). 

Após o crescimento do P. oxalicum CBMAI 1996 em meio de cultura sólido, foram 
inoculados os esporos/micélios transferindo-se 7 discos em meio de cultura líquido 
esterilizado e o pH foi ajustado para 7,4. Os frascos Erlenmeyer foram tampados com 
estofos de algodão/gaze e mantidos para o crescimento celular em agitador orbital (7 dias, 
130 rpm, 32ºC). 

Após os 7 dias de crescimento do P oxalicum CBMAI 1996 separou-se, através de 
uma filtração em funil de Buchner, a massa micelial fúngica do caldo enzimático (produzido 
durante seu período de crescimento). Para as reações de biotransformação da progesterona 
adicionou-se 100 mL do caldo fúngico cultivado (filtrado) em um frasco Erlenmeyer de 250 
mL e adicionou-se 1 mL de solução de progesterona dissolvida em dimetilsufóxido (50 mg/ 
mL). Enquanto que para as reações com a massa micelial fúngica umidecida adicionou-se 
3 g (micélios/esporos) em 100 mL de meio de cultura líquido de extrato de malte recém 
preparado e e adicionou-se 1 mL de solução de progesterona dissolvida em dimetilsufóxido 
(50 mg/ml). 

Realizou-se o acompanhamento das reações a cada 24 h, coletou-se alíquotas 
de 2 mL dos meios, adicionou-se 2 mL de acetato de etila em tubos Falcon e agitou-se 
por 2 minutos em Vórtex. Os tubos Falcon foram levados para centrífuga (6000 rpm, 6 
minutos). Coletou-se a fase orgânica, a qual foi filtrada em presença de sulfato de sódio 
anidro (Na,SO,). Transferiu-se 2 mL desta solução para frascos de 2 mL os quais foram 
injetados e analisados por Cromatografia a Gás acoplada a Espectrometria de Massas 
(CG-EM). Através da integração dos picos obtidos nos cromatogramas, pode-se observar a 
evolução da reação tanto no caldo enzimático (14 dias) quanto com os micélios/esporos ( 
21 dias). Os dados estão apresentados nas Tabelas 2 e 3 e nas Figuras 4 e 5. 


caldo enzimático 
P oxalicum 


CBMAI 1996 
pH 7,4, 32ºC, 
130 rpm, 14 dias 
Progesterona 15B-Hidroxiprogesterona 
Tempo de reação Progesterona cui : 
(dias) (%) 15-Hidroxiprogesterona (%) 
2 93 7 
3 87 13 
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E 86 14 
5 83 17 
6 81 19 
7 75 25 
8 64 36 
9 64 36 
10 68 32 
N 65 35 
12 Té 23 
13 56 44 
14 56 44 


Tabela 2. Monitoramento da reação de biotransformação da progesterona pelo caldo enzimático 
produzido pelo fungo Penicillium oxalicum CBMAI 1996 por 14 dias. 
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Figura 4. Cromatogramas obtidos a partir das análises por CG-EM (IE, 70 eV) da reação de 
biotransformação da progesterona pelo caldo enzimático do fungo P. oxalicum CBMAI 1996 por 14 
dias. (a) Padrão da 7B,15P-di-hidroxiprogesterona. (b) Padrão da 15-hidroxiprogesterona. (c) Padrão 
da progesterona. (d) Reação em 14 dias. (e) Reação em 8 dias. (f) Reação em 4 dias. 
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massa micelial 
P. oxalicum 


CBMAI 1996 = 
pH 7,4, 32ºC, 
130 rpm, 21 dias 
Si uai 15P-Hidroxiprogesterona 7B,15f-di-Hidroxiprogesterona 
Tempo de reação Progesterona 15-Hidroxiprogesterona 7,15-di-Hidroxiprogesterona 

(dias) (%) (%) (%) 
2 93 7 0 
3 60 38 2 
4 16 58 26 
5 7 35 58 
6 1 23 76 
7 1 19 81 
8 0 19 81 
9 0 12 88 
10 0 9 91 
1 0 9 91 
12 0 9 91 
13 0 8 92 
14 0 6 94 
15 0 6 94 
16 0 6 94 
17 0 6 94 
18 0 6 94 
19 0 6 94 
20 0 6 94 
21 0 7 93 


Tabela 3. Monitoramento da reação de biotransformação da progesterona em meio reacional com 
massa a micelial/esporos do fungo Penicillium oxalicum CBMAI 1996 por 21 dias. 
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Figura 5. Cromatogramas obtidos a partir das análises por CG-EM (IE, 70 eV) da reação de 
biotransformação da progesterona em meio reacional com massa micelial do fungo P. oxalicum CBMAI 
1996 por 21 dias. (a) Padrão da progesterona. (b) Padrão da 15f-hidroxiprogesterona. (c) Padrão da 
7B,15B-di-hidroxiprogesterona. (d) Reação em 1 dia. (e) Reação em 7 dias. (f) Reação em 14 dias. (9) 
Extrato final/total da reação (21 dias). 


A partir dos dados obtidos das reações de biotransformação da progesterona 
com o P oxalicum CBMAI 1996, observou-se que com caldo enzimático obteve-se 
majoritariamente o produto mono-hidroxilado (15-hidroxiprogesterona) com conversão de 
44% no 13º dia (Tabela 2, Figura 4). Ao avaliar os dados da reação em meio aos micélios/ 
esporos, observou-se a formação do composto 15-hidroxiprogesterona no primeiro dia de 
reação (24 h), com concomitante decaimento do pico (13,7 min) referente a progesterona. 
O surgimento do pico (14,0 min) referente à formação do composto di-hidroxilado (7,15-di- 
hidroxiprogesterona) ocorreu no 3º dia e evoluiu de forma crescente até o 14º dia de reação, 
a partir do qual se manteve estável até se completar os 21 dias (Tabela 3, Figura 5). 

Na reação com os micélios/esporos a concentração relativa do produto mono- 
hidroxilado atingiu um valor de 58%, a partir do 4º dia. Porém, após este período passou a 
decair, constatando que este composto poderia ser um intermediário da reação na formação 
produto di-hidroxilado. Observou-se ainda que no 8º dia de reação todo o substrato 
(progesterona) foi biotransformado pelo P oxalicum CBMAI 1996 (Tabela 3, Figura 5). 

Concluiu-se que através deste conjunto de reações que o aumento do tempo de 
reação em meio micelial/esporos de 14 dias para 21 dias não promoveu uma melhoria no 
rendimento do produto di-hidroxilado. Constatou-se também que em 4 dias de reação, com 
o uso da massa micelial fúngica, obteve-se um rendimento relativo de 58% para o composto 
mono-hidroxilado, ou seja, 15% mais eficiente do que a reação com caldo enzimático, que 
levou 14 dias para obter uma conversão de 44% (Tabelas 3 e 4, Figuras 4 e 5). 
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Ao final das reações de biotransformação da progesterona (caldo enzimático e 
micélios/esporos) com o fungo P oxalicum CBMAI 1996 realizaram-se as extrações com 
acetato de etila (Figura 6). Os extratos finais obtidos foram purificados por cromatografia 
em coluna cromatográfica, cuja fase estacionária foi composta de sílica gel de porosidade 
de 230-400 mesh. A fase móvel se iniciou com uma mistura de hexano e acetato de etila 
(7:1), que com o decorrer da eluição dos compostos foi sendo enriquecida na composição 
de acetato de etila até ser composta de 100%. 

O processo de purificação por cromatografia em coluna foi monitorado por 
cromatografia em camada delgada (CCD), utilizando cromatofolhas em alumínio (TLC Sílica 
Gel 60 F 


continham os compostos de biotransformação para serem reunidos. Uma vez determinados 


254» Merck) a fim de se verificar os possíveis grupos de frações semelhantes que 
estes grupos de frações por CCD, foram analisadas por CG-EM e as frações obtidas com 
os produtos 7f,15B-di-hidroxiprogesterona e 15P-hidroxiprogesterona foram pesadas e os 
rendimentos foram determinados. As amostras foram dissolvidas em clorofórmio deuterado 
(CDCI,) e analisadas por RMN de 'H (400 MHz) e RMN de *ºC (100 MHz) (Tabelas 4 e 5, 
Figuras 7 e 8). 


a massa filtrada foi adicionada à fase aquosa e realizou-se o procedimento 2 x 


30 ml de = 
acetato de etila à 
—em 
fase aquosa 
OQ isa 


separação das 
agitação por 30 centrifugação a filtração a vácuo fases 


minutos 10.000 rpm, por 
20 minutos 

meio reacional com b fase orgânica 

micélios/esporos 


Coleta de alíquota do 


extrato concentrado e Rotaevaporação 
Análise CG-EM preparo de amostra 


Figura 6. Etapas do processo de extração dos compostos hidroxilados nas reações de 
biotransformação da progesterona com o fungo Penicillium oxalicum CBMAI 1996. 


Fonte: Ligia Breda e Vasconcelos. 


2.1 Caracterização dos produtos de biotransformação da progesterona 


Os espectros de RMN de !'H e de SC dos produtos isolados foram interpretados e 
comparados com os espectros da progesterona. Observou-se, a presença coincidente de 
todos os picos característicos dos carbonos da progesterona, exceto àqueles referentes 
aos carbonos onde ocorreu(ram) a(s) hidroxilação (des), que apresentaram valores com 
descolamentos químicos específicos devidos às inserções dos grupos hidroxilas (Figuras 
7e8). 
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Os valores dos deslocamentos químicos dos carbonos C-7 e C-15 na progesterona 


foram de 32,0 ppm e 24,4 ppm, respectivamente. O valor de deslocamento químico do 


carbono C-7 para o produto di-hidroxilado foi de 74,3 ppm. Enquanto para o carbono 


C-15 foi de 70,2 ppm e 71,6 ppm, para os produtos mono-hidroxilado e di-hidroxilado, 


respectivamente (Tabelas 4 e 5). 


A caracrerização espectral completa dos compostos mono-hidroxilado e di- 


hidroxilado por Ressonância Magnética Nuclear, Infravermelho e Espectometria de Massas 
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Figura 7. Espectro de RMN de 'H (400 MHz, CDCL) e suas ampliações da 7f,15P-di- 


hidroxiprogesterona. (a) Espectro completo. (b) Ampliaç 
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Figura 8. Espectro de RMN de 'ºC (100 MHz, CDCI,) da 7P,15f-di-hidroxiprogesterona. 


Spa duda pa d, 15B-Hidroxiprogesterona Hi pede 
IUPAC/UBMB ida (ppm) (ppm) 
1a 1,71 (m) 1,74 (dd; 4,8, 13,2) 1,68 (d; 4,1) 
1B 2,08 (dd: 3,3, 5,1) 2,05 (m) 2,04 (m) 
2a 2,33 (m) 2,33 (m) 2,41 (m) 
28 2,45 (im) 2,41 (dd; 4,9) 2,44 (m) 
4 5,76 (s) 5,75 (s) 5,77 (m) 
6a 2,32 (ddd; 4,1, 5,5, 14,7) 2,38 (m) 2,55 (m) 
6p 2,40 (m) 2,49 (m) 2,55 (m) 
7a 1,08 (ddd:; 1,2, 4,4, 12,7) 1,14 (dd; 4,0, 13,5) 3,63 (ddd: 6,6, 9,7 9,7) 
7B 1,85 (m) 2,15 (m) - 
8p 1,55 (dd; 3,8, 8,3 10,8) 2,05 (m) 2,02 (m) 
9a 0,99 (m) 1,05 (m) 1,06 (m) 
11a 11,68 (m) 1,64 (m) 1,70 (m) 
118 11,44 (m) 1,52 (m) 1,59 (m) 
12a 1,48 (m) 1,47 (m) 1,40 (m) 
12B 2,07 (m) 2,05 (m) 2,04 (m) 
14a 1,14 (m) 1,05 (m) 1,14 (dd; 5,6, 10,8) 
15a 1,71 (m) 4,32 (ddd; 3,1, 6,0, 6,1) 4,47 (m) 
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15B 
16a 
16B 
17a 
18 
19 
21 


1,28 (m) 
1,68 (m) 
218 (m) 
2,54 (t; 8,6) 
0,67 (s) 
1,19 (s) 
2,13 (s) 


2,26 (m) 
2,26 (m) 
2,49 (m) 
0,94 (s) 
1,22 (m) 
2,15 (m) 


2,34 (m) 
2,27 (m) 
2,44 (m) 
0,97 (m) 
1,25 (m) 
2,14 (s) 


* As multiplicidades dos sinais foram descritas como: “s” para singleto; “?” para tripleto; “dd” para duplo- 
dubleto; “ddd” para duplo-duplo-dubleto; “m” para multipleto. Os valores dispostos após a indicação da 


multiplicidade referem-se às constantes de acoplamento em Hz. 


Tabela 4. Dados de RMN de 'H (400 MHz, CDCI) da progesterona, da 15f-hidroxiprogesterona e da 
7P,15P-di-hidroxiprogesterona. 


Mad do da Progesterona à 15f-Hidroxiprogesterona |, di 
IUPAC/UBMB spp) ppm) (ppm) 
1 35,6 35,7 35,7 
2 34,0 33,9 33,9 
3 199,4 199,4 199,1 
4 123,9 123,9 125,0 
5) 170,9 170,9 166,5 
6 32,8 32,7 42,6 
A 32,0 Si, 74,3 
8 35,6 31,7 38,6 
S) 58;7 54,0 50,9 
10 38,7 38,7 38,0 
1 21,0 20,9 20,8 
12 28,6 40,12 39,6 
13 43,9 43,6 43,9 
14 56,0 60,3 60,1 
15 24,4 70,2 TÁ lo) 
16 22,9 36,2 34,2 
17 68,5 63,7 63,3 
18 13,3 15,9 15,6 
19 17,4 117058) 117,58) 
20 209,3 208,1 208,3 
21 Bio Bi1,9 31,4 


Tabela 5. Dados de RMN de **C (100 MHz, CDCI,) da progesterona, da 15P-hidroxiprogesterona e da 
7PB,15P-di-hidroxiprogesterona. 
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Ao final desse processo, obteve-se cristais dos produtos isolados das respectivas 
frações, os quais foram submetidos para as análises de Raios-X em colaboração com o 
Prof. Dr. Javier Alcides Ellena (IFSC-USP), coordenador do Laboratório Multiusuário de 
Cristalografia Estrutural (LaMuCrEs) que integra o Grupo de Cristalografia do Instituto 
de Física de São Carlos (IFSC). Os dados inéditos de cristalografia de Raios-X para os 
compostos 15P-hidroxiprogesterona e 7f,15f-di-hidroxiprogesterona serão publicados em 
um periódico específico. 

As enzimas envolvidas nas reações de hidroxilação da progesterona nos carbonos 
C-7 e C-15, não são totalmente elucidadas na literatura. Em geral, trata-se de hidroxilação 
enzimática de carbonos não ativados eletronicamente nos esteroides. As enzimas 
envolvidas são monoxigenases do citrocomo P-450, como a esteroide 15B-monoxigenase 
EC 1.14.15.8, que produz a 15f-hidroxiprogesteronal!l, 


2.2 Avaliação da lipofilicidade e ecotoxicidade da 15p-hidroxiprogesterona e 
da 7P,15P-di-hidroxiprogesterona 

O efeito de um composto químico sobre um determinado organismo pode ser agudo, 
quando sua dose letal é liberada em um único evento e rapidamente absorvida, ou crônico, 
quando o composto é liberado em eventos periodicamente sucessivos, em doses subletais, 
durante um determinado período de tempo!" 

A toxicidade de compostos é avaliada a partir das taxas de mortalidade, de 
imobilização ou de inibição de crescimento, dos indivíduos a eles expostos 2, Enquanto 
a ecotoxicidade avalia os efeitos que os compostos químicos apresentam sobre um 
ecossistema, considerando suas interações entre os indivíduos, as comunidades e as 
populações, inclusive ao homem. Compreende-se, dessa forma, como a molécula estudada 
(matéria e energia) interage com toda a cadeia trófica que coexiste no ambiente em que 
se encontra. 

A partir das ferramentas computacionais do site “Chemicalize” 83, estimou- 
se o logaritmo do coeficiente de partição octanol/água (log K,,) da progesterona e dos 
produtos de biotransformação obtidos nesse estudo (15f-hidroxiprogesterona e 7B,15B-di- 
hidroxiprogesterona). Esse coeficiente refere-se à razão entre a afinidade que um composto 
apresenta entre um meio orgânico e o meio aquoso. Quanto maior este valor, maior é a 
afinidade de uma molécula pelo octanol (composto apolar) e, consequentemente, mais 
lipofílica e bioacumulativa é o composto. Pela análise dos dados obtidos, pode-se concluir 
que a biotransformação promoveu uma diminuição da lipofilicidade. Os produtos mono- 
hidroxilado e di-hidroxilado foram aproximadamente 20 e 420 vezes menos lipofílicos, 
respectivamente, do que a progesterona (Tabela 6). 
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Composto log Ky 


Progesterona 4,15 
15B-Hidroxiprogesterona 2,84 
7B,15PB-Hidroxiprogesterona 1,53 


Tabela 6. Dados de log de K,, da progesterona e de seus produtos de biotransformação com o fungo 
Penicillium oxalicum CBMAI 199652, 


Chemaxon. Chemicalize. “Disponível em: https://chemicalize.com/app/calculation”. “Acesso em: 
07/06/2028.” 


Com o auxílio do Programa ECOSAR 2.0 (Ecological Structure Activity 
Relationships) da EPA (United States - Environmental Protection Agency; Em português 
-Agência de Proteção Ambiental dos Estados Unidos), foram realizadas simulações da 
ecotoxicidade dos compostos mono-hidroxilado e di-hidroxilado (Figura 9)º4, O programa 
gera dados baseando-se no efeito causado em uma cadeia trófica, a partir de algas verdes 
(produtores) dafinídeos (espécie de crustáceos pequenos, consumidores primários) peixes 
(consumidores secundários). 

Primeiramente, obtiveram-se os valores dos efeitos agudos: EC, (em mg L”), a 
concentração do composto capaz de gerar efeitos tóxicos a 50% dos indivíduos testados 
(algas verdes) e LC,, (em mg L"'), a concentração letal para 50% dos indivíduos testados 
(dafinídeos e peixes). Em seguida, obtiveram-se os ChV, ou seja, as concentrações dos 
compostos capazes de causar efeitos crônicos estatisticamente significativos 1, 

Pôde-se inferir com a simulação, que a biotransformação da progesterona 
promoveu uma diminuição do seu efeito ecotóxico (Figura 9). A 15B-hidroxiprogesterona 
e a 7B,15B-hidroxiprogesterona foram, aproximadamente, 10 e 100 vezes menos 
ecotóxicos (ecotoxicidade aguda), respectivamente, do que a progesterona. Em relação à 
ecotoxicidade crônica, observou-se, aproximadamente, o mesmo cenário. Como esperado, 
as concentrações estimadas que causam efeitos crônicos foram bem menores do que as 
necessárias para causar os efeitos agudos. 
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Figura 9. Simulação de ecotoxicidade da progesterona e seus produtos de biotransformação pelo 
programa ECOSAR 2.0. (A) Ecotoxicidade aguda. (B) Ecotoxicidade crônica. 


31 CONCLUSÃO E PERSPECTIVAS 


Pequenas mudanças em moléculas com uma atividade biológica intensa, como 
os esteroides, podem promover diversos resultados. As progestinas, derivadas da 
progesterona, são muito mais utilizadas em fármacos que a própria progesterona, devido 
ao aumento de sua biodisponibilidade e a forma como interagem com os receptores nos 
organismos. Dessa forma, os estudos de biotransformação podem levar a obtenção de 
moléculas com novas propriedades biológicas. 
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O emprego de fungos em reações de biotransformação, como o P oxalicum 
CBMAI 1996 foi promissor em nossos estudos frente ao esteroide progesterona. Tratam- 
se de reações altamente específicas e seletivas, pois são catalisadas por enzimas e que 
apresentam um amplo potencial biotecnológico, uma vez que funcionalizaram carbonos 
não reativos eletronicamente (hidridação spº). 

Vale ressaltar que as reações de biotransformação, na maioria das vezes, ocorrem à 
temperatura ambiente e em meio aquoso. Com isso o estudo da reação de biotransformação 
se mostra como potencial para a Química Verde e para a obtenção de novos compostos de 
difícil obtenção por vias químicas sintéticas. 

A partir destes estudos sugere-se que a progesterona pode continuar sendo 
biotransformada no meio ambiente, formando derivados hidroxilados. Assim, o fungo P. 
oxalicum CBMAI 1996 está sendo capaz de metabolizar a progesterona em derivados com 
menor potencial de bioacumulação e com menores efeitos tóxicos ao ecossistema. 

A atividade biológica dos derivados de biotransformação obtidos neste estudo será 
avaliada em trabalhos futuros. E, as análises de raios-X, dos produtos de biotransformação 
foram obtidos, confirmando inequivocamente suas esturutras, bem como a estereoquímica 
absoluta dos centros estereogênicos formados, cujos dados serão publicados na literatura 


científica. 
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RESUMO: Apresenta-se aqui uma 
introdução aos conceitos básicos das águas 
subterrâneas, descrevendo a ocorrência, o 
volume e a qualidade dessas águas, bem 
como o uso e a exploração das mesmas. 
São abordados os diferentes tipos de 
aquíferos e suas funções, além de se 
discutir as áreas de reabastecimento e 
descarga dos aquíferos. Apresenta-se ainda 
uma visão geral das principais ocorrências 
de aquíferos no Brasil. São apresentados 
conceitos teóricos fundamentais, como 
a definição de aquífero, recarga e 
descarga, bem como o papel da litologia 
na formação dos aquíferos. Este conteúdo 
é fundamental para o entendimento das 
águas subterrâneas como recurso hídrico, 
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fornecendo informações importantes para 
gestores públicos e privados, bem como 
para a população em geral, na busca por 
uma gestão adequada e sustentável dos 
recursos hídricos subterrâneos. 

PALAVRAS-CHAVE: Águas Subterrâneas, 
Aquíferos, Hidrogeologia, Brasil 


INTRODUCTION TO 
GROUNDWATER: CONCEPTS AND 
OCCURRENCES IN BRAZIL 
ABSTRACT: Here is an introduction 
to the basic concepts of groundwater, 
describing the occurrence, volume, and 
quality of these waters, as well as their 
use and exploitation. The different types of 
aquifers and their functions are addressed, 
as well as the recharge and discharge 
areas of the aquifers. An overview of 
the main occurrences of aquifers in 
Brazil is also presented. Fundamental 
theoretical concepts are presented, such 
as the definition of aquifer, recharge and 
discharge, as well as the role of lithology in 
aquifer formation. This content is essential 
for understanding groundwater as a water 
resource, providing important information 
for public and private managers, as well 
as the general population, in the search for 
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adequate and sustainable management of groundwater resources. 
KEYWORDS: Groundwater, Aquifers, Hydrology, Brazil 


11 ÁGUAS SUBTERRÂNEAS 


A água subterrânea é toda a água que ocorre abaixo da superfície da Terra, 
preenchendo os poros ou vazios intergranulares das rochas sedimentares, ou as fraturas, 
falhas e fissuras das rochas compactas, e desempenha um papel essencial na manutenção 
da umidade do solo e do fluxo dos rios, lagos e brejos, estando envolvida em uma fase do 
ciclo hidrológico ao constituir uma parcela da água precipitada (Borghetti et al., 2004). 

Após a precipitação, parte das águas que atinge o solo infiltra-se e percola no 
interior do subsolo, em períodos de tempo extremamente variáveis, que dependem de 
muitos fatores, tais como: 

* Porosidade do subsolo: a presença de argila no solo diminui sua permeabilidade, 
não permitindo uma grande infiltração; 


* Cobertura vegetal: um solo coberto por vegetação é mais permeável do que um 
solo desmatado; 


* Inclinação do terreno: em declividades acentuadas a água corre mais rapidamente, 
diminuindo a possibilidade de infiltração; 


* Tipo de chuva: chuvas intensas saturam rapidamente o solo, ao passo que chuvas 
finas e demoradas têm mais tempo para se infiltrarem. 

Durante a infiltração, uma parcela da água fica retida nas regiões mais próximas 
da superfície do solo, constituindo a zona não saturada, enquanto outra parcela atinge 
as zonas mais profundas do subsolo, constituindo a zona saturada, sob a ação da força 
de adesão ou de capilaridade e da gravidade, respectivamente (Tundisi, 2005). Como é 
possível observar no modelo conceitual da Figura 1. 


Cobertura do Solo 


eso 


Lençol Freático (Limite) 

Zona Insaturada 
————» 

Zona Saturada 


>) 


Ea 
Água Subterrânea 


Figura 1- Modelo conceitual das zonas não saturada e saturada no subsolo. 


Fonte: O autor. 
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Durante esse processo a água é mantida perto das partículas superficiais da rocha, 
pela atração molecular, sejam estas rochas sedimentares ou não. (Figura 2). O divisor limite 
do nível da linha de água, passa a ser chamado de lençol freático, e todas as aberturas 
abaixo deste limite estarão teoricamente cheias de águas subterrâneas. 


Figura 2 - Caracterização Esquemática do comportamento da Água em Zonas Saturadas. 


Fonte: O Autor. 


Zona Insaturada: Também chamada de zona não saturada, zona de aeração ou 
vadosa, é a parte do solo que está parcialmente preenchida por água. Nesta zona, pequenas 
quantidades de água distribuem-se uniformemente, sendo que suas moléculas aderem às 
superfícies dos grãos do solo. Nesta zona ocorrem os fenômenos de transpiração pelas 
raízes das plantas, filtração e autodepuração da água (Freeze et al., 1979). Dentro desta 
zona encontram-se: 

* Zona de umidade do solo: é a parte mais superficial, onde a perda de água de 
adesão para a atmosfera é intensa. Em alguns casos, é muito grande a quantidade 
de sais que se precipitam na superfície do solo após a evaporação dessa água, 
dando origem a solos salinizados ou a crostas ferruginosas (lateríticas). Esta zona 


serve de suporte fundamental da biomassa vegetal natural ou cultivada da Terra e 
da interface atmosfera/litosfera (Freeze et al., 1979). 


* Zona intermediária: região compreendida entre a zona de umidade do solo e a 
franja capilar, com umidade menor do que nesta última e maior do que a da zona 
superficial do solo. Em áreas onde o nível freático está próximo da superfície, a zona 
intermediária pode não existir, pois a franja capilar atinge a superfície do solo. São 
brejos e alagadiços, onde há uma intensa evaporação da água subterrânea (Freeze 
et al., 1979). 


* Franja de capilaridade: é a região mais próxima ao nível d'água do lençol freático, 
onde a umidade é maior devido à presença da zona saturada logo abaixo (Freeze 
et al., 1979). 


Zona saturada: é a região abaixo da zona não saturada onde os poros ou fraturas 
da rocha estão totalmente preenchidos por água. As águas atingem essa zona por 
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gravidade, através dos poros ou fraturas, até alcançar uma profundidade limite, onde as 
rochas estão tão saturadas que a água não pode penetrar mais. Para que haja infiltração 
até a zona saturada, é necessário primeiro satisfazer as necessidades da força de adesão 
na zona não saturada. Nessa zona, a água corresponde ao excedente de água da zona 
não saturada que se move em velocidades muito lentas (em/dia), formando o manancial 
subterrâneo propriamente dito. Uma parcela dessa água irá desaguar na superfície dos 
terrenos, formando as fontes e os olhos d'água. A outra parcela desse fluxo subterrâneo 
forma o caudal basal que deságua nos rios, perenizando-os durante os períodos de 
estiagem, com uma contribuição multianual média da ordem de 13.000 kmº/ano (Peixoto e 
Oort, 1990), ou deságua diretamente nos lagos e oceanos. A superfície que separa a zona 
saturada da zona de aeração é chamada de nível freático, ou seja, esse nível corresponde 
ao topo da zona saturada (IGM, 2001). 

Dependendo das características climatológicas da região ou do volume de 
precipitação e escoamento da água, esse nível pode permanecer permanentemente a 
grandes profundidades ou se aproximar da superfície horizontal do terreno, originando as 
zonas encharcadas ou pantanosas, ou convertendo-se em mananciais (nascentes) quando 
se aproxima da superfície através de um corte no terreno (Rebouças, 1996). 


21 OCORRÊNCIA E VOLUME DAS ÁGUAS SUBTERRÂNEAS 


Assim como a distribuição das águas superficiais é muito variável, a das águas 
subterrâneas também o é, uma vez que elas se inter-relacionam no ciclo hidrológico e 
dependem das condições climatológicas. Entretanto, as águas subterrâneas (10.360.230 
km?) são aproximadamente 100 vezes mais abundantes do que as águas superficiais dos 
rios e lagos (92.168 km?). (WSS 2019). 

Embora estejam armazenadas nos poros e fissuras milimétricas das rochas, as 
águas subterrâneas ocorrem em grandes extensões, gerando volumes significativos de 
água subterrânea da ordem de aproximadamente 23.400.000 km?º, distribuídos em uma 
área aproximada de 134,8 milhões de km? (Shiklomanov, 1998; Hanasaki et al., 2020), 
constituindo-se em importantes reservas de água doce. 

Alguns especialistas indicam que a quantidade de água subterrânea pode chegar a 
até 60 milhões de kmº, mas sua ocorrência em grandes profundidades pode impossibilitar 
seu uso. Por essa razão, a quantidade passível de ser captada estaria a menos de 4.000 
metros de profundidade, compreendendo cerca de 8 e 10 milhões de kmº (CEPIS, 2000), 
que, segundo Rebouças et al. (2002), estaria assim distribuída: 65.000 km3 constituindo a 
umidade do solo; 4,2 milhões de kmº desde a zona não-saturada até 750 m de profundidade 
e 5,3 milhões de kmº de 750 m até 4.000 m de profundidade, constituindo o manancial 
subterrâneo. 
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Além disso, a quantidade de água capaz de ser armazenada pelas rochas e pelos 
materiais não consolidados em geral depende da porosidade dessas rochas, que pode 
chegar a até 45% (IGM, 2001), da comunicação desses poros entre si ou da quantidade e 
tamanho das aberturas de fraturas existentes. 

No Brasil, de acordo com o Ministério do Meio Ambiente (2021), as reservas de água 
subterrânea são estimadas em 110.000 kmº (110 trilhões de m?) e a contribuição multianual 
média à descarga dos rios é da ordem de 2.850 kmº/ano. Nem todas as formações 
geológicas possuem características hidrodinâmicas que possibilitem a extração econômica 
de água subterrânea para atendimento de médias e grandes vazões pontuais. As vazões já 
obtidas por poços variam, no Brasil, desde menos de 1 m%/h até mais de 2.000 mº/h (ANA, 
2021). 


31 QUALIDADE DAS ÁGUAS SUBTERRÂNEAS 


Durante o percurso em que a água percola entre os poros do subsolo e das rochas, 
ocorre a depuração da mesma por meio de uma série de processos físico-químicos (troca 
iônica, decaimento radioativo, remoção de sólidos em suspensão, neutralização de pH em 
meio poroso, entre outros) e bacteriológicos (eliminação de micro-organismos devido à 
ausência de nutrientes e oxigênio que os viabilizem), que atuam sobre a água, modificando 
as suas características adquiridas anteriormente e tornando-a particularmente mais 
adequada para o consumo humano (Silva, 2003; Albuquerque et al., 2020). 

Sendo assim, a composição química da água subterrânea é o resultado combinado 
da composição da água que entra no solo e da evolução química influenciada diretamente 
pelas litologias atravessadas. O teor de substâncias dissolvidas nas águas subterrâneas 
aumenta à medida que elas se movem (SMA,2003; Sousa et al., 2018). 

As águas subterrâneas apresentam algumas propriedades que tornam o seu uso mais 
vantajoso em relação ao das águas dos rios: são filtradas e purificadas naturalmente por 
meio da percolação, o que determina excelente qualidade e dispensa tratamentos prévios; 
não ocupam espaço em superfície; sofrem menor influência das variações climáticas; são 
passíveis de extração perto do local de uso; possuem temperatura constante; têm maior 
quantidade de reservas; necessitam de custos menores como fonte de água; suas reservas 
e captações não ocupam área superficial; apresentam grande proteção contra agentes 
poluidores; o uso do recurso aumenta a reserva e melhora a qualidade; possibilitam a 
implantação de projetos de abastecimento à medida da necessidade (Werege, 1997; 
ABAS, 2016). 


41 USO E EXPLORAÇÃO DAS ÁGUAS SUBTERRÂNEAS 


A busca por água subterrânea continua sendo uma prática importante na atualidade, 
devido à crescente demanda de água e à necessidade de garantir a segurança hídrica. 


Ciências exatas e da terra: teorias e princípios 2 Capítulo 6 


106 


Os habitantes continuam buscando fontes alternativas, como poços e furos profundos 
(Trimmer, 2000; Foster et al., 2018). A exploração de água subterrânea é condicionada a 
fatores quantitativos, qualitativos e econômicos (Howard, 2017; Leal, 1999), como: 


* Quantidade: influenciada pela taxa de recarga dos aquíferos e pelas caracte- 
rísticas hidráulicas dos terrenos; 


* Qualidade: influenciada pela composição das rochas, pela presença de conta- 
minantes e pelas condições climáticas e de renovação das águas; 


* | Econômico: depende do custo da perfuração de poços, da profundidade do 
aquífero, da qualidade da água e das condições de bombeamento. 

O uso de águas subterrâneas é uma prática antiga e sua evolução tem acompanhado 
a própria evolução da humanidade. As técnicas de construção de poços e métodos de 
bombeamento têm sido aprimoradas ao longo do tempo, possibilitando a extração de água 
em volumes e profundidades cada vez maiores, atendendo às demandas de cidades, 
indústrias, projetos de irrigação, entre outros (Foster et al., 2018). 

A relação entre a demanda e o uso varia entre os países e, dentro deles, de região 
para região. No entanto, em geral, o abastecimento público continua sendo a maior 
demanda individual por água subterrânea (PROASNE, 2003; Barlow et al., 2019). 

Praticamente todos os países do mundo, desenvolvidos ou não, utilizam água 
subterrânea para suprir suas necessidades. Países como Alemanha, Áustria, Bélgica, 
Dinamarca, França, Holanda, Hungria, Itália, Marrocos, Rússia e Suíça atendem de 70% 
a 90% da demanda para o abastecimento público (Margat e Van der Gun, 2013). Outros 
países utilizam a água subterrânea para atender a todas as suas necessidades (como 
Dinamarca, Arábia Saudita e Malta) ou apenas como suplemento para o abastecimento 
público e atividades como irrigação, produção de energia, turismo e indústria. Na Austrália, 
60% do país depende totalmente do manancial subterrâneo e em mais de 20% o seu uso 
é predominante. Na região metropolitana da Cidade do México, cerca de 80% da demanda 
dos quase 20 milhões de habitantes é atendida (Margat e Van der Gun, 2013). 

De acordo com a UNESCO (2021), cerca de 2 bilhões de pessoas dependem da 
água subterrânea para sobreviver em todo o mundo. Regiões áridas e semiáridas, bem 
como algumas ilhas, têm a água subterrânea como sua única fonte hídrica disponível para 
uso humano. Mesmo regiões desérticas, como a Líbia, atendem à demanda de água em 
cidades e na irrigação por meio de poços tubulares perfurados no Saara (UNESCO, 2021). 

De acordo com a FAO (2021), estima-se que foram perfurados mais de 500 milhões 
de poços no mundo, a maioria dos quais nos países em desenvolvimento, para fins de 
irrigação e abastecimento doméstico. Nos Estados Unidos, por exemplo, são perfurados 
cerca de 150 mil poços por ano, com uma extração de cerca de 100 bilhões de m3/ano, 
atendendo a mais de 40% do abastecimento público e das indústrias (USGS, 2021). Em 
todo o mundo, cerca de 340 bilhões de toneladas de água são retiradas anualmente na 
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Átrica do Norte, China, Índia, Estados Unidos e Arábia Saudita, e não se renovam. Essa 
água poderia produzir comida suficiente para alimentar mais de 1 bilhão de pessoas por 
ano (Shah, 2016). 

A busca intensa por água, sobretudo potencializada pela expansão das áreas 
agrícolas tem provocado o uso intenso das águas subterrâneas, além das fontes superficiais. 
Existem diversos exemplos no mundo de esgotamento de aquíferos por sobrexploração 
para uso em irrigação (Gleick, 2020). Estima-se que existam no mundo 350 milhões de 
hectares irrigados com água subterrânea, sendo 21 milhões desses nos Estados Unidos e 
70 milhões na Índia (FAO, 2021). 

No que se refere ao Brasil, vários núcleos urbanos dependem da água subterrânea 
como recurso principal ou complementar, incluindo indústrias, propriedades rurais, escolas, 
hospitais e outros estabelecimentos que utilizam água de poços profundos com frequência. 
O abastecimento público também é bastante dependente dessa fonte de água, e importantes 
cidades brasileiras, como Ribeirão Preto (SP), Mossoró e Natal (RN), Maceió (AL), Região 
Metropolitana de Recife (PE) e Barreiras (BA), dependem integral ou parcialmente da água 
subterrânea para seu abastecimento (ANA, 2022). 

Portanto, a gestão sustentável dos recursos hídricos subterrâneos é essencial 
para garantir a segurança hídrica no país. As águas subterrâneas desempenham um 
papel fundamental no abastecimento de muitas cidades brasileiras, sendo que em alguns 
estados, como o Maranhão, São Paulo e Piauí, a dependência dessas águas chega a 
atingir percentuais superiores a 70% e 80%. Além de fornecer água potável, as águas 
subterrâneas também apresentam outras características únicas, como a presença de 
águas termais, que estimulam o turismo em cidades como Caldas Novas em Goiás, Araxá 
e Poços de Caldas em Minas Gerais. Outro uso importante das águas subterrâneas é 
a utilização da água mineral pelas populações dos centros urbanos devido à sua alta 
qualidade (ANA, 2022). 

Mesmo em áreas com elevado teor salino, como nas regiões de ocorrência dos 
sistemas aquíferos fissurados do semiárido nordestino, as águas subterrâneas muitas 
vezes são a única fonte de suprimento permanente, apesar dos desafios associados ao 
seu uso (Oliveira et al., 2019). 

Segundo o Censo do IBGE de 2010, aproximadamente 53% da população brasileira 
é abastecida, para fins domésticos, com água subterrânea, sendo que 5% se auto-abastece 
das águas de poços rasos, 10% de nascentes ou fontes e 38% de poços profundos. 
Portanto, o número de poços tubulares em operação no Brasil está estimado em cerca de 
450.000, com um número anual de perfurações de aproximadamente 15.000, o que pode 
ser considerado insuficiente diante das necessidades de água potável das populações e 
se comparado com outros países (ANA, 2022). Os estados com maior número de poços 
perfurados são São Paulo (60.000), Bahia, Rio Grande do Sul, Ceará e Piauí (Oliveira et. 
al., 2019). 
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A exploração de água subterrânea é realizada de várias maneiras, variando de 
métodos tradicionais a modernos. Tais métodos estão se desenvolvendo todos os dias 
através de novos meios e dispositivos. A utilização de tecnologias como a geofísica, que 
emprega técnicas como eletrorresistividade, sísmica de reflexão, sísmica de refração, entre 
outras, tem ganhado cada vez mais espaço na exploração de água subterrânea (ANA, 
2022). 

A exploração de água subterrânea é realizada de várias maneiras, variando de 
métodos tradicionais a modernos. Tais métodos estão em constante desenvolvimento 
através de novos meios e dispositivos. A utilização de tecnologias como a geofísica, que 
emprega técnicas como eletrorresistividade, sísmica de reflexão, sísmica de refração, entre 
outras, tem ganhado cada vez mais espaço na exploração de água subterrânea (ANA, 
2022; Meijirik, 2007) 

Várias técnicas podem ser utilizadas para explorar os recursos hídricos 
subterrâneos. A perfuração de teste e a análise estratigráfica são os métodos mais 
confiáveis e padrão para determinar a localização de um poço e a espessura da unidade 
aquífera (Sophocleous, 2018; Madan et al., 2010). No entanto, esses métodos de 
investigação da água subterrânea não são eficazes em termos de tempo e custo, e muitas 
vezes exigem habilidades específicas (Sophocleous, 2018; Roscoe, 1990; Fetter, 1994). 

As técnicas de prospecção geofísica também têm sido utilizadas por vários 
pesquisadores para explorar os recursos hídricos subterrâneos em diferentes tipos de 
terreno geológico (Singh et al., 2017; Zhao et al., 2021; Ako et al., 1989; Amadi et al., 1990; 
Olorunfemi et al., 1995; Olayinka et al., 2001 e Adiat et al., 2009). No entanto, devido à 
falta de precisão de uma análise no local, os resultados e interpretações de levantamentos 
geofísicos sempre requerem validação com dados de furos (Ghazanfari et al., 2018; Adiat, 
2009). 

O uso de técnicas modernas como Sensoriamento Remoto e Sistema de 
Informação Geográfica tem fornecido meios mais eficazes, rápidos e econômicos para 
avaliar e gerenciar os recursos hídricos subterrâneos (Sophocleous, 2018; Singh et al., 
2021; Zhao et al., 2021; Jha et al., 2007; Meijerink, 2007). Essas técnicas permitem a 
identificação de áreas promissoras para exploração de água subterrânea, com base na 
avaliação de parâmetros hidrológicos, geológicos e topográficos relevantes para a sua 
disponibilidade. Estudos têm sido realizados para investigar a aplicação dessas técnicas 
em diferentes regiões e condições geológicas (Ghazanfari et al., 2018; Jha et al., 2019). 


51 AQUÍFEROS 


Um aquífero é uma formação geológica do subsolo, constituída por rochas 
permeáveis que armazenam água em seus poros ou fraturas. Outro conceito se refere ao 


aquífero como sendo somente o material geológico capaz de servir como depositório e 
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transmissor da água armazenada. A água subterrânea é armazenada em espaços abertos 
e fraturas dentro de formações geológicas abaixo da superfície da terra, conhecidos como 
aquíferos (Lehr et al., 2005). 

Um aquífero é um leito saturado ou formação que não só armazena água, mas a 
produz em quantidade suficiente para ser importante como fonte de abastecimento. Assim, 
uma litologia só será considerada um aquífero se, além de ter seus poros saturados de 
água, permitir a fácil transmissão da água armazenada. 

A principal propriedade de um aquífero é a sua capacidade de liberar a água retida 
em seus poros e sua capacidade de transmitir o fluxo facilmente (Lehr et al., 2005). Essas 
propriedades essencialmente dependem da composição do aquífero e incluem a porosidade, 
rendimento, condutividade hidráulica, permeabilidade e coeficiente de armazenamento. 

Um aquífero pode ter extensão de poucos quilômetros quadrados a milhares de 
quilômetros quadrados, ou pode, também, apresentar espessuras de poucos metros a 
centenas de metros (Rebouças et al., 2002). Etimologicamente, aquífero significa: aqui = 
água; fero = transfere; ou do grego, suporte de água (Heinen et al., 2009). 

Os aquíferos podem ser formados por rochas consolidadas ou não consolidadas 
(Lehr et al., 2005). As rochas consolidadas são compostas de materiais como arenito, 
tufo, calcário e granito, já as rochas inconsolidadas são compostas por sedimentos não 
compactados, como areia, argila, cascalho e outras partículas que se acumulam ao longo 
do tempo e formam camadas sedimentares. 

Os aquíferos mais importantes do mundo, seja por extensão ou pela 
transnacionalidade, são: o Guarani — Argentina, Brasil, Paraguai, Uruguai (1,2 milhões de 
km?); o Arenito Núbia - Líbia, Egito, Chade, Sudão (2,5 milhões de km?); o Sistema Aquífero 
Guaraní - Argentina, Brasil, Paraguai, Uruguai (1,2 milhões de km?); o Digitalwaterway 
Vechte — Alemanha, Holanda (7,5 mil km?); o Slovak Karst-Aggtelek - República Eslovaca e 
Hungria; o Praded — República Checa e Polônia (3,3 mil km?) (UNESCO, 2021); a Grande 
Bacia Artesiana (1,7 milhões km?) e a Bacia Murray (375 mil km?), ambos na Austrália. Em um 
levantamento mais recente, a UNESCO identificou mais de 200 aquíferos transfronteiriços 
em todo o mundo, incluindo mais de 100 na Europa (UNESCO, 2021). 


61 TIPOS DE AQUÍFEROS 


A litologia do aquífero, ou seja, a sua constituição geológica (porosidade/ 
permeabilidade intergranular ou de fissuras), é que irá determinar a velocidade da água 
em seu meio, a qualidade da água e a sua capacidade como reservatório. A Iitologia é 
resultado da sua origem geológica, que pode ser fluvial, lacustre, eólica, glacial e aluvial 
(rochas sedimentares), vulcânica (rochas fraturadas) e metamórfica (rochas calcárias), 
determinando assim os diferentes tipos de aquíferos. Quanto à porosidade, existem três 
tipos de aquíferos, os fissurais, porosos e cársticos. (Sophocleous, 2018; Fetter, 2001; 
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Freeze e Cherry, 1979; Todd e Mays, 2005; Domenico e Schwartz, 1997; Bear e Cheng, 
2010). (Figura 3) 


o 
= 


= 


Figura 3- Tipos de Aquíferos quanto à porosidade. 


Fonte: O Autor. 


* | Aquífero poroso ou sedimentar — é aquele formado por rochas sedimentares 
consolidadas, sedimentos inconsolidados ou solos arenosos, onde a circulação 
da água se dá nos poros formados entre os grãos de areia, silte e argila de 
granulação variada. Constituem os mais importantes aquíferos, pelo grande vo- 
lume de água que armazenam e por sua ocorrência em grandes áreas. Esses 
aquíferos ocorrem nas bacias sedimentares e em todas as várzeas onde se 
acumularam sedimentos arenosos. Uma particularidade desse tipo de aquífero 
é sua porosidade quase sempre homogeneamente distribuída, permitindo que 
a água flua em qualquer direção, em função somente dos diferenciais de pres- 
são hidrostática ali existentes. Essa propriedade é conhecida como isotropia.( 
Fetter, 2001; Freeze e Cherry, 1979; Todd e Mays, 2005; Domenico e Schwartz, 
1997; Bear e Cheng, 2010) 


* — Aquífero fraturado ou fissural — é formado por rochas ígneas, metamórficas 
ou cristalinas, duras e maciças, onde a circulação da água se dá nas fraturas, 
fendas e falhas abertas devido ao movimento tectônico. Ex.: basalto, granitos, 
gabros, filões de quartzo, etc. (SMA, 2003). A capacidade dessas rochas de 
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acumularem água está relacionada à quantidade de fraturas, suas aberturas 
e intercomunicação, permitindo a infiltração e fluxo da água. Poços perfurados 
nessas rochas fornecem poucos metros cúbicos de água por hora, sendo que 
a possibilidade de se ter um poço produtivo dependerá somente de intercep- 
tar fraturas capazes de conduzir a água. Nesses aquíferos, a água só pode 
fluir onde houver fraturas, que, quase sempre, tendem a ter orientações pre- 
ferenciais. São ditos, portanto, aquíferos anisotrópicos. Um caso particular de 
aquífero fraturado é representado pelos derrames de rochas vulcânicas basá- 
lticas das grandes bacias sedimentares brasileiras. (Melo et. al. 2017; Frost e 
Zaghloul 2013) 


* — Aquífero cárstico (Karst) — é formado em rochas calcárias ou carbonáticas, 
onde a circulação da água se dá nas fraturas e outras descontinuidades (dia- 
clases) que resultaram da dissolução do carbonato pela água. Essas aberturas 
podem atingir grandes dimensões, criando verdadeiros rios subterrâneos. São 
aquíferos heterogêneos, descontínuos, com águas duras e com fluxo em ca- 
nais. As rochas são os calcários, dolomitos e mármores. (Doerfliger et. al. 2007, 
Ford e Williams 2007) 


Os aquíferos também podem ser classificados quanto à superfície superior, segundo 
a pressão da água (Figura 4), podendo ser classificado em aquífero livre, confinado e 
semiconfinados. 


Afloramento Elevação 


, do Lençol 
Solo | Recarga Freático 
Recarga FluxoemPoço Recarga ; 


ençol Freático 


Nível do 


Figura 4- Tipos de Águifero quanto a pressão. 


Fonte: O autor. 
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Aquífero livre: constituído por uma formação geológica permeável e superfi- 
cial, totalmente exposta em toda a sua extensão e limitada na base por uma 
camada impermeável. A superfície superior da zona saturada está em equilíbrio 
com a pressão atmosférica e a recarga ocorre diretamente pela chuva. Em um 
poço de aquífero livre, o nível da água varia de acordo com a quantidade de 
chuva e é conhecido como nível freático. É o aquífero mais comum e explorado, 
mas também é vulnerável à contaminação. (Amaral, 2012; Makridakis e Dianat, 
2020) 


Aquífero confinado: constituído por uma formação geológica permeável, con- 
finado entre duas camadas impermeáveis ou semipermeáveis. A pressão da 
água no topo da zona saturada é maior do que a pressão atmosférica, fazendo 
com que a água suba no poço acima do nível do aquífero. A recarga indireta 
ocorre preferencialmente nos locais onde a formação aflora à superfície. Em 
um poço de aquífero confinado, a água sobe acima do teto do aquífero devi- 
do à pressão exercida pelas camadas confinantes sobrejacentes. A altura que 
a água sobe é conhecida como nível potenciométrico, e o poço é artesiano. 
(Makridakis e Dianat, 2020; Fetter, 2001) 


Aquífero semi-confinado: limitado na base, no topo ou em ambos, por cama- 
das cuja permeabilidade é menor do que a do aquífero em si. O fluxo preferen- 
cial da água se dá ao longo da camada do aquífero. A recarga pode ocorrer de 
forma secundária através das camadas semi-confinantes, desde que haja uma 
diferença de pressão hidrostática entre a camada do aquífero e as camadas 
subjacentes ou sobrejacentes. Em certas circunstâncias, um aquífero livre pode 
ser abastecido por água oriunda de camadas semi-confinadas subjacentes, ou 
vice-versa. Zonas de fraturas ou falhas geológicas podem constituir pontos de 
fuga ou recarga da água da camada confinada. (Makridakis e Dianat, 2020; 
Fetter, 2001) 


71 ÁREAS DE REABASTECIMENTO E DESCARGA DOS AQUÍFEROS 


Um aquífero apresenta uma reserva permanente de água e uma reserva ativa ou 


reguladora que são continuamente abastecidas através da infiltração da chuva e de outras 


fontes subterrâneas. As reservas reguladoras ou ativas correspondem ao escoamento de 


base dos rios. 


A área por onde ocorre o abastecimento do aquífero é chamada zona de recarga, 


que pode ser direta ou indireta. O escoamento de parte da água do aquífero ocorre na zona 
de descarga (ANA, 2001). 


Zona de recarga direta: é aquela onde as águas das chuvas infiltram-se 
diretamente no aquífero, através de suas áreas de afloramento e fissuras de 
rochas sobrejacentes. Sendo assim, a recarga sempre é direta nos aquíferos 
livres, ocorrendo em toda a superfície acima do lençol freático. Nos aquíferos 
confinados, o reabastecimento ocorre preferencialmente nos locais onde a 
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formação portadora de água aflora à superfície (Dill et. al., 2012; Braga et. al. 
2011). 


* Zona de recarga indireta: são aquelas onde o reabastecimento do aquífero 
se dá a partir da drenagem (filtração vertical) superficial das águas e do fluxo 
subterrâneo indireto, ao longo do pacote confinante sobrejacente, nas áreas 
onde a carga potenciométrica favorece os fluxos descendentes (Dill et. al., 
2012; Braga et. al. 2011). 


* Zona de descarga: é aquela por onde as águas emergem do sistema, alimen- 
tando rios e jorrando com pressão por furos (Dill et. al., 2012; Braga et. al. 2011). 

As maiores taxas de recarga ocorrem em regiões planas, bem arborizadas, e nos 
aquíferos livres. Em regiões de relevo acidentado, sem cobertura vegetal, sujeitas a práticas 
de uso e ocupação que favorecem as enxurradas, a recarga ocorre mais lentamente e de 
maneira limitada (Rebouças et al., 2002). 

Sob condições naturais, apenas uma parcela dessas reservas reguladoras é passível 
de exploração, constituindo o potencial ou reserva explorável. Em geral, esta parcela 
é calculada entre 25% e 50% das reservas reguladoras (ANA, 2001). Esse volume de 
exploração pode aumentar em função das condições de ocorrência e recarga, bem como 
dos meios técnicos e financeiros disponíveis, considerando que a soma das extrações 
com as descargas naturais do aquífero para rios e oceano não pode ser superior à recarga 
natural do aquífero. 


81 FUNÇÕES DOS AQUÍFEROS 


Além de suprir água suficiente para manter os cursos de água superficiais estáveis 
(função de produção), os aquíferos também ajudam a evitar seu transbordamento, 
absorvendo o excesso de água da chuva intensa (função de regularização). Na Ásia 
tropical, onde a estação quente pode durar até 9 meses e as chuvas de monção podem 
ser bastante intensas, esse duplo serviço hidrológico é crucial (SAMPAT, 2001). Segundo o 
mesmo autor, os aquíferos também proporcionam uma forma de armazenar água doce sem 
muita perda pela evaporação - outro serviço particularmente valioso em regiões quentes, 
propensas à seca, onde essas perdas podem ser extremamente altas. 

Na África, por exemplo, em média, um terço da água extraída de reservatórios 
todo ano é perdido pela evaporação. Os pântanos, habitats importantes para aves, peixes 
e outras formas de vida silvestre, nutrem-se normalmente de água subterrânea, onde o 
lençol freático aflora à superfície em ritmo constante. Onde há muita exaustão de água 
subterrânea, o resultado é frequentemente leitos secos de rios e pântanos ressecados. 
(Lovell, 2017; Said e Thompson, 2011) 

Portanto, os aquíferos podem cumprir as seguintes funções (Rebouças et al., 2002): 
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* Função de produção: corresponde à sua função mais tradicional de produção 
de água para o consumo humano, industrial ou irrigação. 


* Função de estocagem e regularização: utilização do aquífero para estocar 
excedentes de água que ocorrem durante as enchentes dos rios, correspon- 
dentes à capacidade máxima das estações de tratamento durante os períodos 
de demanda baixa, ou referentes ao reuso de efluentes domésticos e/ou indus- 
triais. 


* Função de filtro: corresponde à utilização da capacidade filtrante e de depura- 
ção bio-geoquímica do maciço natural permeável. Para isso, são implantados 
poços a distâncias adequadas de rios perenes, lagoas, lagos ou reservatórios, 
para extrair água naturalmente clarificada e purificada, reduzindo substancial- 
mente os custos dos processos convencionais de tratamento. 


* Função ambiental: consiste em sua capacidade de manter a qualidade da 
água subterrânea, protegendo-a de contaminação e poluição. Além disso, os 
aquíferos desempenham um papel importante no equilíbrio do ecossistema lo- 
cal, fornecendo água para pântanos, riachos e rios, que por sua vez sustentam 
a biodiversidade de animais e plantas. 


* Função transporte: o aquífero é utilizado como um sistema de transporte de 
água entre zonas de recarga artificial ou natural e áreas de extração excessiva. 


* Função estratégica: a água contida em um aquífero foi acumulada durante 
muitos anos ou até séculos e é uma reserva estratégica para épocas de pouca 
ou nenhuma chuva. 


91 OCORRÊNCIAS DE AQUÍFEROS NO BRASIL 


A combinação das estruturas geológicas com os fatores geomorfológicos e 
climáticos do Brasil resultou na configuração de 10 províncias hidrogeológicas (Figura 5), 
que são regiões com sistemas aquíferos com condições semelhantes de armazenamento, 
circulação e qualidade da água (MMA, 2003). Essas províncias podem ser divididas em 


subprovíncias. 


Ciências exatas e da terra: teorias e princípios 2 Capítulo 6 


115 


Províncias Hidrogeológicas 


BB Amazonas 


Centro-Oeste 
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1:25.000.000 DD Paraná 
o 195 390 780 1170 1.560 [] Parnaiba 
Km EB são Francisco 
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co mmuAty 
TT 


Figura 5 - Províncias Hidrogeológicas. 


Fonte: Adaptado de Pessoa et. al 1980. 


Conforme Feitosa et al., 2008, uma Província Hidrogeológica é uma região com 
características gerais semelhantes em relação às principais ocorrências de águas 
subterrâneas. Segundo Feitosa et al. (2008), os fatores fundamentais para definir uma 
Província Hidrogeológica são os aspectos geológicos (litologia, estrutura e tectônica) e 
fisiográficos (clima, morfologia, hidrografia, solos e vegetação). É importante ressaltar, no 
entanto, que uma Província Hidrogeológica pode ser modificada ou subdividida ao longo 
do tempo. 

Sendo assim, as águas subterrâneas no Brasil ocupam diferentes tipos de 
reservatórios, desde as zonas fraturadas do embasamento cristalino (escudo) até os 
depósitos sedimentares cenozoicos (bacias sedimentares), sendo classificadas em três 
sistemas aquíferos: porosos, fissurados e cársticos, conforme a Tabela 1 (Leal, 1999). 
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Províncias Domínios Sistema Aquífero Área Volume de Total (%) 
Hidrogeológicas Aquíferos Principal (Km?) Agua (Km?) o 
] Substratos 
Escudo Oriental Aflorantes Zonas Fraturadas 600000 80 0,07 
Escudos Manto Rochas 
Setentrional, Central Pa Fraturadas e/ou | 4000000 10000 8,90 
e Meridional Alteradas 
Bacia ; ; 
: Arenitos Barreiras e 
Amazonas Edi Alter do Chão 1300000 32500 28,94 
mazonas 
Arenitos Itapecuru, 
Bacia Cordas Grajaú, 
Parnaíba Sedimentar Motuca, Poti-Piaufí, 700000 17500 15,58 
Maranhão Cabeças, Serra 
Grande 
Bacia 
E Sedimentar Arenitos São Luis e 
Parnaíba São Luis- Itapecuru 50000 250 0,22 
Barrerinhas 
Es e Arenitos Barreiras, 
Costeira Potiguar Açu- Beberibe, 2300 260 0,20 
9 Calcário Jandaía 
Recife 
Bacia 
. Sedimentar Arenitos Barreiras e 
Costeria Sergipe Nantubá 10000 100 0,9 
-Alagoas 
Bacia 
Sedimentar Arenitos Marizal, 
Costeira Jatobá, Tacaratú, São 5600 840 0,75 
Tucano- Sebastião 
Recôncavo 
Arenitos Baurú- 
Paraná Eai “Aqutjauana, 100000 50400 44,88 
Paraná Guarani, Rio Branco, À 
Bassaltos Serra 
Geral 
Depósitos ES 
Diversos Aluviões e Dunas 77300 411 0,37 


Tabela 1- Províncias Hidrogeológicas e Sístemas Aquíferos (MMA 2023). 


Os escudos são compostos por rochas magmáticas e metamórficas e representam 


os primeiros núcleos de rochas emergidas que surgiram durante a formação inicial da 
crosta terrestre. Por outro lado, as bacias sedimentares são depressões que, ao longo do 
tempo, foram preenchidas por detritos ou sedimentos provenientes de áreas próximas ou 
distantes, geralmente dispostas de forma horizontal (Coelho, 1996). 


As Províncias são classificados em domínios hidrogeológicos, os quais são 


definidos como “representativos de um conjunto de unidades geológicas com similaridades 


hidrogeológicas, baseando-se principalmente nas características litológicas das rochas” 
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(Feitosa et al., 2008). Esse conceito foi desenvolvido durante a elaboração do Mapa de 
Domínios Hidrogeológicos da Bahia, que ainda não foi publicado (Bomfim, 2010). 

O Mapa de Domínios/Subdomínios Hidrogeológicos do Brasil (Figura 6), elaborado 
pela CPRM (Bomfim e Jesus, 2007), utilizando o Sistema de Informações Geográficas 
(SIG), de acordo com a definição mencionada anteriormente, estabelece sete domínios 
hidrogeológicos: 


* | Formações Cenozoicas: aquífero poroso; 

* Bacias Sedimentares: aquífero poroso; 

* | Poroso/Fissural: aquífero misto; 

* | Metassedimentos/Metavulcânicas: aquífero fissural; 
* | Vulcânicas: aquífero fissural; 

* Cristalino: aquífero fissural; 


* | Carbonatos/Metacarbonatos: aquífero fissural. 
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Domínios Hidrogeológicas 


Formações Cenozoicas 
[D) Bacias Sedimentares 
EB Porosorrissura 
CHILE * Metasedimento- Meta Vulcânica 
— NVulcânica 
0) cristalino 


(BB carbonatos-Metacarbonatos 


1:25.000.000 
780 11470 1.560 
Km 


Esri, HERE, Garmin, (e) OpenStreetMap contributors, and the GIS user 
community 


Figura 6- Domínios Hidrogeológicos. 


Fonte : Adaptado Bonfin e Jesus 2007. 


Os domínios hidrogeológicos no Brasil abrangem uma ampla variedade de sistemas 
hidrogeológicos, cada um com suas características distintas. Esses sistemas podem ser 
classificados com base em suas propriedades hidráulicas, geológicas e geográficas, sendo 
essencial compreender sua natureza para uma gestão adequada dos recursos hídricos 
subterrâneos. 

Um dos sistemas hidrogeológicos dominantes no país são os aquíferos porosos, 
formados por rochas sedimentares que ocupam cerca de 42% da área total (3,6 milhões de 
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Km?). Esses aquíferos estão distribuídos em cinco províncias hidrogeológicas: Amazonas, 
Paraná, Parnaíba-Maranhão, Centro-Oeste e Costeira. Sua estrutura geológica, com 
camadas permeáveis e impermeáveis alternadas, confere-lhes condições de artesianismo. 
Destacam-se as Bacias do Paraná, Amazonas, Parnaíba e a Subprovíncia Potiguar-Recife, 
devido à sua extensão e potencialidade. (ABAS, 2003) 

Outro sistema hidrogeológico presente no país é o dos aquíferos fissurados ou 
fraturados, que abrangem cerca de 53,8% do território nacional. Esses aquíferos estão 
nas Províncias Hidrogeológicas dos Escudos Setentrional, Central, Oriental e Meridional. 
As duas primeiras províncias, com rochas fraturadas do embasamento, possuem boas 
possibilidades hídricas devido aos altos índices pluviométricos. A Província Oriental está 
dividida em duas subprovíncias: Nordeste e Sudeste. Já a Província Meridional, localizada 
em Santa Catarina e no Rio Grande do Sul, apresenta substrato alterado. Os rios perenes 
e os aquíferos dessa região são recarregados pelos altos índices pluviométricos locais. 
(MMA,2003). 

As duas primeiras províncias com rochas fraturadas do embasamento apresentam 
razoáveis possibilidades hídricas, devido aos altos índices pluviométricos da área. A 
Província Oriental está dividida em duas subprovíncias (Nordeste e Sudeste). 

A Província Meridional, em Santa Catarina e no Rio Grande do Sul é de substrato 
alterado. Os altos índices pluviométricos da região asseguram a perenização dos rios e 
contribuem para a recarga dos aquíferos, cujas reservas são, em parte, restituídas à rede 
hidrográfica (MMA,2008). 

Esse sistema apresenta reservas de águas subterrâneas da ordem de 10.080 kmº 
(Teixeira et. al. 2020; Bittencourt e Savio, 2020; Leal, 1999). As águas são de boa qualidade 
química, podendo ocorrer localmente teores de ferro acima do permitido. No domínio do 
embasamento cristalino subaflorante, como na Província Hidrogeológica Escudo Oriental 
do Nordeste onde está localizada a região semiárida — há pequena disponibilidade hídrica, 
devido à formação de rochas cristalinas (Lima e Conicelli, 2020). É frequente observar teor 
elevado de sais nas águas dessa região, o que restringe ou impossibilita seu uso (MMA, 
2003). Nesse domínio subaflorante é que nascem os rios temporários. 

Além disso, temos os sistemas hidrogeológicos cársticos, compostos pelo sistema 
cárstico-fissural da Província Hidrogeológica do São Francisco e pela Formação Jandaíra. 
Inclui-se também os domínios do calcário do Grupo Bambuí com mais de 350.000 km, 
distribuído nos Estados da Bahia, Goiás e Minas Gerais, e a Formação Caatinga. As 
profundidades de desenvolvimento cárstico variam, com média em torno de 150 metros. 
O Grupo Bambuí apresenta vazões superiores a 200 metros cúbicos por hora, enquanto 
a Formação Jandaíra possui vazões mais baixas, geralmente inferiores a 3,5 mºjh. Outros 
aquíferos cársticos relevantes são o Pirabas, com profundidade média de 220 metros e 
vazão de 135 mºjh, e a Formação Capiru do Grupo Açungui, com vazão média de 180 mºjh 
e profundidade média de 60 m. (Carmo e Costa, 2021; Lisboa e Simões, 2021; MMA,20083). 
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A classificação e compreensão desses sistemas hidrogeológicos são fundamentais 


para uma gestão eficiente e sustentável dos recursos hídricos subterrâneos. 


Os sistemas aquíferos brasileiros abrigam reservatórios significativos de água 
(Figura 7), que alimentam uma das redes de rios perenes mais extensas do mundo, com 
exceção dos rios temporários que se originam nas regiões semiáridas do Nordeste, onde 
afloram as rochas do embasamento geológico (Rebouças et al., 2002). Além disso, esses 
sistemas desempenham um papel socioeconômico importante devido ao seu potencial 


hídrico (MMA, 2003). 
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Figura 7- Principáis Sistemas Aquíferos do Brasil. 
Fonte: Adaptado MMA. 
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Os sistemas aquíferos no Brasil são abundantes e desempenham um papel crucial 
no abastecimento de água e no equilíbrio hidrológico do país. A seguir, apresento os 
principais aquíferos por regiões: 

Região Norte: Aquífero Alter do Chão: Localizado na Bacia Amazônica, abrange 
uma área de aproximadamente 437.000 km? com uma vazão aproximada de 86 mil kms, 
possuindo grande importância na região, fornecendo água para diversos rios perenes. 
(Sousa et al. 2020) 

Região Nordeste: O aquífero Urucuia Estende-se pelos estados de Piauí, Ceará, 
Bahia e Minas Gerais, sendo um dos principais aquíferos da região Nordeste. Possui grande 
potencial hídrico, contribuindo para o abastecimento de água de várias cidades. Com 
extensão de aproximadamente 142 mil km?, e vazão de 420 mº/s. (Gonçalves et. al 2016) 

O aquífero Jandaíra, localizado no Rio Grande do Norte e em parte da Paraíba e 
do Ceará, é um importante aquífero cárstico-fissural, responsável pelo abastecimento de 
água de diversas comunidades da região. Ocorre numa superfície total de 17.756 km? 
esta área de ocorrência inclui as zonas recobertas pela Formação Barreiras, aluviões e 
dunas, que totalizam cerca de 5.980 km?. Estudos regionais indicam que a espessura da 
Formação Jandaíra pode chegar a 600 m na porção mais profunda da Bacia Potiguar, 
embora as espessuras mais importantes, do ponto de vista hidrogeológico, tenham 
sido registradas entre 50 e 150 m com vazões variando entre 10 a 50 mº h'! conforme 
Vasconscelos et.al 2010. 

Região Centro-Oeste: O aquífero Guarani, considerado um dos maiores sistemas 
aquíferos do mundo, estende-se por partes dos estados do Mato Grosso, Mato Grosso 
do Sul, Goiás e São Paulo, além de abranger áreas do Paraguai, Argentina e Uruguai. É 
uma fonte estratégica de água doce, utilizada para abastecimento público, agricultura e 
indústria. A maior parte (70 por cento ou 840 000 km?) da área ocupada pelo aquífero — 
cerca de 1 200 000 km? — está no subsolo do centro-sudoeste do Brasil. O restante se 
distribui entre o nordeste da Argentina (255 000 km?), noroeste do Uruguai (58 500 km?) e 
sudeste do Paraguai (58 500 km?), nas bacias do rio Paraná e do Chaco-Paraná. A vazão 
é muito variável. Por exemplo: é de mais de 200 000 litros por hora na região do Alto Rio 
Uruguai, no Rio Grande do Sul, mas possui raros poços acima de 5 000 litros por hora na 
região das Missões, no mesmo estado. Em outros estados, já foram registradas vazões 
da ordem de 800 000 litros por hora. (Borges et. al 2013; Borges et. al. 2020) 

Região Sudeste: O Aquífero Bauru, localizado principalmente no estado de São 
Paulo, abrange uma área significativa no Centro-Oeste. É uma importante reserva de água 
subterrânea utilizada para múltiplos fins, como abastecimento urbano e irrigação. Ocupa 
uma área de com área total de 370.000 km?, as camadas de sedimentação intermediária 
e de topo mostram vazão variável entre 3,0 e 20 mº/h, enquanto porções mais restritas na 
base do pacote sedimentar exibem vazões entre 20 e 50 m*/h. De modo geral, a média da 
vazão de produção dos poços é de 8 m* /h (DAEE, 2005; Birsa et. al. 2019). 
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O aquífero Tubarão: Situado no estado de São Paulo, é um aquífero sedimentar que 
desempenha um papel fundamental no abastecimento de água para a região. Seus poços 
normalmente apresentam baixa produtividade, com fluxos inferiores a 2,8 litros/segundo. 
Contudo, em alguns pontos é possível atingir vazões de até 11 litros/segundo, com uma 
área de 20.000 Km? (CETESB,2010; Silva et. al. 2020). 

Região Sul: O aquífero Serra Geral, estende-se pelos estados de São Paulo, 
Paraná e Santa Catarina. É um importante aquífero fissural que fornece água para várias 
cidades e atividades agrícolas na região. Os basaltos fissurados afloram numa extensão 
de cerca de 20.000 km?, localizado em camadas inferiores aos sedimentos do Grupo 
Bauru. Sua espessura varia desde poucos metros, aumentando para Oeste, até 1.000 
metros. Assim sendo, suas características hidrodinâmicas ficam melhor demonstradas 
pelos valores de vazão específica que variam entre 0,08 e 50 m%/h/m., com valor médio 
de 1,0 mº /h/m (Fernandes et al. 2006, Gonçalves et. al. 2020). 

Os sistemas aquíferos no Brasil desempenham um papel fundamental na 
disponibilidade de água doce e no funcionamento dos rios perenes em diversas 
regiões do país. Eles abrigam reservatórios significativos de água e contribuem para o 
desenvolvimento socioeconômico, fornecendo recursos hídricos para atividades agrícolas, 
industriais e domésticas. Além disso, a compreensão dos diferentes sistemas aquíferos e 
suas características, como extensão, vazão e potencialidade, é essencial para o manejo 
sustentável dos recursos hídricos subterrâneos. A preservação e conservação desses 
sistemas são de suma importância para garantir a disponibilidade de água de qualidade no 
presente e para as gerações futuras. 


10 | CONSIDERAÇÕES FINAIS 


Ao considerarmos todas as informações discutidas sobre águas subterrâneas, 
podemos concluir que elas desempenham um papel crucial no abastecimento de água 
doce, tanto para as necessidades humanas quanto para a manutenção dos ecossistemas. 
A ocorrência e o volume das águas subterrâneas são influenciados por fatores geológicos, 
climáticos e hidrológicos, e seu uso sustentável requer uma compreensão detalhada de 
suas características e comportamento. 

A exploração de águas subterrâneas é uma prática comum em muitas regiões, 
proporcionando uma fonte confiável de água para consumo humano, agricultura e indústria. 
No entanto, é essencial que essa exploração seja feita de forma responsável, considerando 
a recarga natural dos aquíferos e evitando o esgotamento excessivo dos recursos hídricos. 

A qualidade das águas subterrâneas é de extrema importância, pois elas podem 
ser vulneráveis à contaminação por atividades humanas e processos naturais. Portanto, 
é fundamental adotar medidas de proteção e monitoramento para preservar a qualidade 
desses recursos vitais. 
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Existem diferentes tipos de aquíferos no Brasil, como aquíferos porosos, fissurados 
e cársticos, cada um com suas características hidrodinâmicas específicas. Além disso, é 
essencial destacar as áreas de reabastecimento e descarga dos aquíferos, que influenciam 
o fluxo e a disponibilidade das águas subterrâneas. 

Os aquíferos desempenham funções essenciais na regulação do ciclo hidrológico, 
armazenando e liberando água conforme necessário. Eles contribuem para a manutenção 
dos rios perenes, atuam como reservatórios estratégicos e têm um papel fundamental na 
sustentabilidade hídrica. 

No Brasil, há ocorrências significativas de aquíferos, como o Aquífero Guarani, o 
Aquífero Alter do Chão, o Aquífero Urucuia, entre outros. Esses aquíferos fornecem água 
para múltiplos usos, contribuindo para o desenvolvimento socioeconômico das regiões 
onde estão localizados. 

Em suma, a compreensão dos aspectos relacionados às águas subterrâneas, 
incluindo sua ocorrência, volume, exploração, qualidade, tipos de aquíferos, áreas de 
reabastecimento e descarga, funções e ocorrências no Brasil, é fundamental para uma 
gestão adequada e sustentável desses recursos preciosos, garantindo seu uso consciente 
e a preservação desses importantes ecossistemas. 
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